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FIELD OF THE INVENTION 

The present invention is directed generally to a system and method for locating people or 

15 objects, and in particular, to a system and method for locating a wireless mobile station using a plurality of 
mobile station location estimators. More generally, the present invention is directed to a computational 
system and method for calibrating the relative performance of multiple models, wherein each such model 
is capable of being activated for generating hypotheses (e.g., estimates and/or predictions) of an 
unknown condition. Additionally, the present invention is directed to a computational system and method 

20 for generating enhanced hypotheses of the unknown condition, wherein the model generated hypotheses 
are used as queries into an archive that associates: (a) historical model generated hypotheses, (b) model 
input data used in generating the model hypotheses, and (c) verified hypotheses to which the model input 
data is known to correspond. 

BACKGROUND 
25 Introduction 

Wireless communications systems are becoming increasingly important worldwide. Wireless 
cellular telecommunications systems are rapidly replacing conventional wire-based telecommunications 
systems in many applications. Cellular radio telephone networks ("CRT"), and specialized mobile radio 
and mobile data radio networks are examples. The general principles of wireless cellular telephony have 
30 been described variously, for example in U. S. Patent 5,295,180 to Vendetti, et al, which is incorporated 
herein by reference. 

There is great interest in using existing infrastructures for wireless communication systems for 
locating people and/or objects in a cost effective manner. Such a capability would be invaluable in a 
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variety of situations, especially in emergency or crime situations. Due to the substantial benefits of such a 
location system, several attempts have been made to design and implement such a system. 

Systems have been proposed that rely upon signal strength and ^lateralization techniques to 
permit location include those disclosed in U.S. Patents 4,818,998 and 4,908,629 to Apsell et al. ("the 
5 Apsell patents") and 4,891 ,650 to Sheffer ("the Sheffer patent*). However, these systems have 

drawbacks that include high expense in that special purpose electronics are required. Furthermore, the 
systems are generally only effective in line-of-sight conditions, such as rural settings. Radio wave surface 
reflections, refractions and ground clutter cause significant distortion, in determining the location of a 
signal source in most geographical areas that are more than sparsely populated. Moreover, these 

10 drawbacks are particularly exacerbated in dense urban canyon (city) areas, where errors and/or conflicts 
in location measurements can result in substantial inaccuracies. 

Another example of a location system using time of arrival and triangulation for location are 
satellite-based systems, such as the military and commercial versions of the Global Positioning Satellite 
system ("GPS"). GPS can provide accurate position determination (i.e., about 100 meters error for the 

1 5 commercial version of GPS) from a time-based signal received simultaneously from at least three 
satellites. A ground-based GPS receiver at or near the object to be located determines the difference 
between the time at which each satellite transmits a time signal and the time at which the signal is 
received and, based on the time differentials, determines the object's location. However, the GPS is 
impractical in many applications. The signal power levels from the satellites are low and the GPS 

20 receiver requires a clear, line-of-sight path to at least three satellites above a horizon of about 60 degrees 
for effective operation. Accordingly, inclement weather conditions, such as clouds, terrain features, such 
as hills and trees, and buildings restrict the ability of the GPS receiver to determine its position. 
Furthermore, the initial GPS signal detection process for a GPS receiver is relatively long (i.e., several 
minutes) for determining the receiver's position. Such delays are unacceptable in many applications such 

25 as, for example, emergency response and vehicle tracking. 

Summary of Factors Affecting RF Propagation 

The physical radio propagation channel perturbs signal strength, frequency (causing rate 
changes, phase delay, signal to noise ratios (e.g., C/l for the analog case, or Emo , RF energy per bit, 
over average noise density ratio for the digital case) and Doppler-shift. Signal strength is usually 
30 characterized by: 

• Free Space Path Loss (Lp) 

• Slow fading loss or margin (Uow) 
■ Fast fading loss or margin (Ust) 

Loss due to slow fading includes shadowing due to clutter blockage (sometimes included in Lp). 
35 Fast fading is composed of multipath reflections which cause: 1 ) delay spread; 2) random phase shift or 
Rayleigh fading; and 3) random frequency modulation due to different Doppler shifts on different paths. 
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Summing the path loss and the two fading margin loss components from the above yields a total 
path loss of: 

Uotal = Lp + Ulow + Lfast 



5 - Referring to Fig. 3, the figure illustrates key components of a typical cellular and PCS power budget 

design process. The cell designer increases the transmitted power P T x by the shadow fading margin Ukw 
which is usually chosen to be within the 1-2 percentile of the slow fading probability density function 
(PDF) to minimize the probability of unsatisfactorily low received power level Prx at the receiver. The Prx 
level must have enough signal to noise energy level (e.g., 10 dB) to overcome the receiver's internal 

io noise level (e.g., -1 18dBm in the case of cellular 0.9 GHz), for a minimum voice quality standard. Thus in 
the example Prx must never be below -108 dBm, in order to maintain the quality standard. 

Additionally the short term fast signal fading due to multipath propagation is taken into account 
by deploying fast fading margin Uast, which is typically also chosen to be a few percentiles of the fast 
fading distribution. The 1 to 2 percentiles compliment other network blockage guidelines. For example the 

15 cell base station traffic loading capacity and network transport facilities are usually designed for a 1-2 
percentile blockage factor as well. However, in the worst-case scenario both fading margins are 
simultaneously exceeded, thus causing a fading margin overload. 

SUMMARY 

OBJECTS OF THE DISCLOSURE 

20 It is an objective of the present disclosure to present a system and method for to wireless 

telecommunication systems for accurately locating people and/or objects in a cost effective manner. 
Additionally, it is an objective of the present disclosure to present such location capabilities using the 
measurements from wireless signals communicated between mobile stations and a network of base 
stations, wherein the same communication standard or protocol is utilized for location as is used by the 

25 network of base stations for providing wireless communications with mobile stations for other purposes 
such as voice communication and/or visual communication (such as text paging, graphical or video 
communications). Related objectives for the present disclosure include presenting a system and method 
that: 

(1.1) can be readily incorporated into existing commercial wireless telephony systems with few, if any, 
30 modifications of a typical telephony wireless infrastructure; 

(1 .2) can use the native electronics of typical commercially available telephony wireless mobile 
stations (e.g., handsets) as location devices; 

(1 .3) can be used for effectively locating people and/or objects wherein there are few (if any) line-of- 
sight wireless receivers for receiving location signals from a mobile station (herein also denoted MS); 

35 (1 .4) can be used not only for decreasing location determining difficulties due to multipath phenomena 
but in fact uses such multipath for providing more accurate location estimates; 
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(1 .5) can be used for integrating a wide variety of location techniques in a straight-forward manner; 

(1 .6) can substantially automatically adapt and/or (re)train and/or (re)calibrate itself according to 

changes in the environment and/or terrain of a geographical area where the novel method and/or system 
of the present disclosure is utilized; and 
5 (1 .7) - can utilize a plurality of wireless location estimators based on different wireless location 
technologies (e.g., GPS location techniques, terrestrial base station signal timing techniques for 
triangulation and/or trilateration, wireless signal angle of arrival location techniques, techniques for 
determining a wireless location within a building, techniques for determining a mobile station location 
using wireless location data collected from the wireless coverage area for, e.g., location techniques using 

10 base station signal coverage areas, signal pattern matching location techniques and/or stochastic 
techniques), wherein each such estimator may be activated independently of one another, whenever 
suitable data is provided thereto and/or certain conditions, e.g., specific to the estimator are met; 
(1 .8) can provide a common interface module from which a plurality of the location estimators can be 
activated and/or provided with input; 

15 (1 .9) provides resulting mobile station location estimates to location requesting applications (e.g., for 
91 1 emergency, the fire or police departments, taxi services, vehicle location, etc.) via an output gateway, 
wherein this gateway: 

(a) routes the mobile station location estimates to the appropriate location application(s) 
via a communications network such as a wireless network, a public switched 

20 telephone network, a short messaging service (SMS), and the Internet, 

(b) determines the location granularity and representation desired by each location 
application requesting a location of a mobile station, and/or 

(c) enhances the received location estimates by, e.g., performing additional processing 
such as "snap to street" functions for mobile stations known to reside in a vehicle. 

25 Yet another objective is to provide a low cost location system and method, adaptable to wireless 

telephony systems, for using simultaneously a plurality of location techniques for synergistically 
increasing MS location accuracy and consistency. 

It is yet another objective that at least some of the following MS location techniques can be 
utilized by various embodiments of the novel method and/or system of the present disclosure: 

30 (2.1 ) time-of-arrival wireless signal processing techniques; 

(2.2) time-difference-of-arrival wireless signal processing techniques; 

(2.3) adaptive wireless signa> processing techniques having, for example, learning capabilities and 
including, for instance, artificial neural net and genetic algorithm processing; 

(2.4) signal processing techniques for matching MS location signals with wireless signal 
35 characteristics of known areas; 

(2.5) conflict resolution techniques for resolving conflicts in hypotheses for MS location estimates; 
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(2.6) enhancement of MS location estimates through the use of both heuristics and historical data 
associating MS wireless signal characteristics with known locations and/or environmental conditions. 

Yet another objective is to provide the capability to activate serially one or more groups of one or 
more location estimating models, wherein each model may generate a location estimate for locating the 
5 same MS. Thus, for example, if a first group of models did not yield an appropriate location estimate, 
then another group may be activated. Accordingly, a first group may include a model based on global 
positioning satellite (GPS) technology, whereas a second group (or other group) may include base station 
triangulation models and pattern recognition models. 

Yet another objective is to provide location estimates in terms of time vectors, which can be used 
10 to establish motion, speed, and an extrapolated next location in cases where the MS signal subsequently 
becomes unavailable. 

DEFINITIONS 

The following definitions are provided for convenience. In general, the definitions here are also 
defined elsewhere in this document as well. 
15 (3.1) The term "wireless" herein is, in general, an abbreviation for "digital wireless", and in particular, 
"wireless" refers to digital radio signaling using one of standard digital protocols such as CDMA, NAMPS, 
AMPS, TDMA and GSM, as one skilled in the art will understand. 

(3.2) As used herein, the term "mobile station" (equivalently, MS) refers to a wireless device that is at 
least a transmitting device, and in most cases is also a wireless receiving device, such as a portable radio 

20 telephony handset. Note that in some contexts herein instead or in addition to MS, the following terms are 
also used: "personal station" (PS), and "location unit" (LU). In general, these terms may be considered 
synonymous. However, the later two terms may be used when referring to reduced functionality 
communication devices in comparison to a typical digital wireless mobile telephone. 

(3.3) The term, "infrastructure", denotes the network of telephony communication services, and more 
25 particularly, that portion of such a network that receives and processes wireless communications with 

wireless mobile stations. In particular, this infrastructure includes telephony wireless base stations (BS) 
such as those for radio mobile communication systems based on CDMA, AMPS, NAMPS, TDMA, and 
GSM wherein the base stations provide a network of cooperative communication channels with an air 
interface with the MS, and a conventional telecommunications interface with a Mobile Switch Center 

30 (MSC). Thus, an MS user within an area serviced by the base stations may be provided with wireless 
communication throughout the area by user transparent communication transfers (i.e., "handoffs") 
between the user's MS and these base stations in order to maintain effective telephony service. The 
mobile switch center (MSC) provides communications and control connectivity among base stations and 
the public telephone network. 

35 (3.4) The phrase, "composite wireless signal characteristic values" denotes the result of aggregating 
and filtering a collection of measurements of wireless signal samples, wherein these samples are 
obtained from the wireless communication between an MS to be located and the base station 
infrastructure (e.g., a plurality of networked base stations). However, other phrases are also used herein 



to denote this collection of derived characteristic values depending on the context and the likely 
orientation of the reader. For example, when viewing these values from a wireless signal processing 
perspective of radio engineering, as in the descriptions of the subsequent Detailed Description sections 
concerned with the aspects of the present disclosure for receiving MS signal measurements from the 
5 base station infrastructure, the phrase typically used is: "RF signal measurements". Alternatively, from a 
data processing perspective, the phrases: "location signature dusted and "location signal data" are used 
to describe signal characteristic values between the MS and the plurality of infrastructure base stations 
substantially simultaneously detecting MS transmissions. Moreover, since the location communications 
between an MS and the base station infrastructure typically include simultaneous communications with 
1 o more than one base station, a related useful notion is that of a "location signature" which is the composite 
wireless signal characteristic values for signal samples between an MS to be located and a single base 
station. Also, in some contexts, the phrases: "signal characteristic values" or "signal characteristic data" 
are used when either or both a location signature(s) and/or a location signature clusters) are intended. 

SUMMARY DISCUSSION 

1 5 In one embodiment, the present disclosure relates to a novel method and system for performing 

wireless mobile station location. In particular, the present disclosure presents a wireless mobile station 
location computing method and system that utilizes multiple wireless location computational estimators 
(these estimators also denoted herein as MS location hypothesizing computational models, "first order 
models" and/or location estimating models"), for providing a plurality of location estimates of a target 

20 mobile station, wherein ambiguities and/or conflicts between the location estimates may be effectively and 
straightforwardly resolved. More particularly, the present disclosure presents a technique for calibrating 
the performance of each of the location estimators so that a confidence value (e.g., a probability) can be 
assigned to each generated location estimate. Additionally, the present disclosure presents provides a 
straightforward technique for using the confidence values (probabilities) for deriving a resulting most likely 

25 location estimate of a target wireless mobile station. 

More generally, the present disclosure relates to a novel computational method and architecture 
for synergistically combining the results of a plurality of computational models in a straightforward way 
that allows the models to be calibrated relative to one another so that differences in results generated by 
the models can be readily resolved. Accordingly, the computational method and architecture of the 

30 present disclosure may be applied to wide range applications where synergies between multiple models 
are expected to be enhance performance. 

For a particular application having a plurality of computational models (each generating an 
estimate of a desired result(s) in a space of hypothesis results), the novel method and/or system 
presented herein may be described, at a high level, as any method or system that performs the following 

35 steps: 

(4.1 .1 ) A step of determining a classification scheme for determining an input class for each 

input data set supplied (e.g., substantially simultaneously) to the plurality of 

computational models (FOMs), wherein for each range, R, of a plurality ranges of 
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desired results in the hypothesis space, there is an input class, and the input data sets 
of this input class are expected to have their corresponding desired result(s) in the 
range R. Some examples will be illustrative. For a wireless location system, the 
present step determines geographical subareas of a wireless network coverage area 
that have "similar* wireless signal characteristics. Such subareas may be relatively 
easy to determine, and there may be no constraint on the size of the subareas. The 
intention is to determine: (a) such a subarea as only a general area where a target MS 
must reside, and (b)the subarea should be relatively homogeneous in its wireless 
signaling characteristics. Accordingly, (a) and (b) are believed to be substantially 
satisfied by grouping together into the same input class the wireless signal data sets 
(i.e., input data sets) from corresponding target MS locations wherein at each of the 
target MS locations: (i) the set of base stations detected by the target MS (at the 
location) is substantially the same, and/or (b) the set of base stations detecting the 
target MS is substantially the same set of base stations. 

Note that there are numerous techniques and commercial packages for 
determining such a classification scheme. In particular, the statistically based system, 
"CART" (an acronym for Classification and Regression Trees) by ANGOSS Software 
International Limited of Toronto, Canada is one such package. Further, note that this 
step is intended to provide reliable but not necessarily highly accurate ranges R for the 
desired results. Also note that in some applications there may be only a single input 
class, thus assuring high reliability (albeit, likely low accuracy). Accordingly, in this 
latter case the present step may be omitted. 
(4.1 .2) A step of calibrating each of the plurality of computational models (FOMs) so that each 
subsequent hypothesis generated by one of the models has a confidence value (e.g., 
probability) associated therewith that is indicative of the likeliness of the hypothesis 
being correct. The calibrating of this step is performed using the input classification 
scheme determined in the above step (4.1 .1). In one embodiment of this step, each 
model is supplied with inputs from a given fixed input class, wherein each of these 
inputs have corresponding known results that constitute a correct hypothesis (i.e., a 
desired result). Subsequently, the performance of each model is determined for the 
input class and a confidence value is assigned to the model for inputs received from the 
input class. Note that this procedure is repeated with each input class available from 
the input classification scheme. In performing this procedure, an application domain 
specific criteria is used to determine whether the hypotheses generated by the models 
identify the desired results in the hypothesis space. Accordingly, for each of the 
models, when supplied with an input data set from a fixed input class, the hypothesis 
generated by the model will be given the confidence value determined for this input 
class as an indication of the likelihood of the generated hypothesis being correct (i.e., 



the desired result). Note that the confidence value for each generated hypothesis may 
be computed as a probability that the hypothesis is correct. 

Note that for a wireless location application, the criteria (in one embodiment) is 
whether a location hypothesis contains the actual location where the MS was when the 
corresponding input data set (wireless signal measurements) were communicated 
between this MS and the wireless network. 

For applications related to the diagnosis of electronic systems, this criteria may 
be whether an hypothesis identifies a proper functional unit such as a circuit board or 
chip. 

For economic forecasting applications, this criteria may be whether an 
hypothesis is within a particular range of the correct hypothesis. For example, if an 
application according to the novel method and/or system disclosed herein predicts the 
U.S. gross national product (GNP) six months into the future according to certain inputs 
(defining input data sets), then hypotheses generated from historical data that has 
associated therewith the actual corresponding GNP (six months later), may be used for 
calibrating each of the plurality of economic forecasting models (FOMs). Thus, the 
application specific criteria for this case may be that a generated hypothesis is within, 
say, 10% of the actual corresponding six month GNP prediction. 

Note that the applications described herein are illustrative, but not 
comprehensive of the scope of the novel method and/or system disclosed herein. 
Further note that this step typically is performed at least once prior to inputting input 
data sets whose resulting hypotheses are to be used to determine the desired or 
correct results. Additionally, once an initial calibration has been performed, this step 
may also be performed: (a) intermittently between the generation of hypotheses, and/or 
(b) substantially continuously and in parallel with the generation of hypotheses by the 
models. 

(4.1 .3 ) A step of providing one or more input data sets to the models (FOMs) for generating a 

plurality of hypotheses, wherein the result(s) desired to be hypothesized are unknown. 
Moreover, note that the generated hypotheses are preferred to have a same data 
structure definition. 

For example, for a wireless location system, the present step provides an input 
data set including the composite signal characteristic values to one or more MS location 
hypothesizing computational models, wherein each such model subsequently 
determines one or more initial estimates (also denoted location hypotheses) of the 
location of the target MS. Note that one or more of these models may be based on, for 
example, the signal processing techniques 2.1 through 2.3 above. 

(4.1 .4) A step of adjusting or modifying the generated hypotheses output by the models, wherein 

for such an hypothesis, adjustments may be performed on one or both of its 
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hypothesized result H.R, and its confidence value for further enhancing the 
performance of the novel method and/or system of the present disclosure. In one 
embodiment of this step, H.R is used as an index to retrieve other results from an 
archival database, wherein this database associates hypothesized results with their 
corresponding desired or correct results. Thus, H.R may be used to identify data from 
other archived hypothesized results that are "nearby" to H.R, and subsequently use the 
nearby data to retrieve the corresponding desired results. Thus, the set of retrieved 
desired results may be used to define a new "adjusted" hypothesis. 

For example, for a wireless location system utilizing the novel method and/or 
system of the present disclosure, each location hypothesis, H, identifies an area for a 
target MS, and H can used to identify additional related locations included in archived 
hypotheses generated by the same FOM as generated H. For instance, such related 
locations may be the area centroids of the archived hypotheses, wherein these 
centroids reside within the area hypothesized by H. Accordingly, such centroids may 
be used to retrieve the corresponding actual verified MS locations (i.e., the 
corresponding desired results), and these retrieved verified locations may be used to 
generate a new adjusted area that is likely to be more accurate than H. In particular, a 
convex hull of the verified locations may be used as a basis for determining a new 
location hypothesis of the target MS. 

For other application domains, the present step requires a first technique to 
determine both "nearby" archived data from previously archived hypotheses, and a 
second technique to determine an "adjusted" hypothesis from the retrieved desired 
results. In general, such techniques can be relatively straightforward to provide when 
the hypothesized results reside in a vector space, and more particularly, in a Cartesian 
product of the real numbers. Accordingly, there are numerous applications that can be 
configured to generate hypothesized results in a vector space (or Cartesian product of 
the real numbers). For instance, economic financial forecasting applications typically 
result in numeric predictions where the first and second techniques can be, e.g., 
substantially identical to the centroid and convex hull techniques for the wireless 
location application.; and 
(4.1 .5 ) A step of subsequently computing a "most likely" estimate is computed, for outputting to 
a requesting application ( e.g., a location requesting application, such as 911 
emergency, the fire or police departments, taxi services, etc.). Note that in computing 
the most likely target MS location estimate a plurality of location hypotheses may be 
taken into account. In fact, it is an important aspect of the novel method and/or system 
disclosed herein that the most likely MS location estimate is determined by 
computationally forming a composite MS location estimate utilizing such a plurality of 



location hypotheses so that, for example, location estimate similarities between location 
hypotheses can be effectively utilized. 
Referring to (4.1 .3) there may be hypotheses for estimating not only desired result(s), but also 
hypotheses may be generated that indicate where the desired result(s) is not. Thus, if the confidence 
5 values are probabilities, an hypothesis may be generated that has a very low (near zero) probability of 
having the desired result. As an aside, note that in general, for each generated hypothesis, H, having a 
probability, P, there is a dual hypothesis H c that may be generated, wherein the H c represents the 
complementary hypothesis that the desired result is in the space of hypothesized results outside of H. 
Thus, the probability that the desired result(s) is outside of the result hypothesized by H is 1-P. 
l o Accordingly, with each location hypothesis having a probability favorably indicating where a desired result 
may be (i.e., P >= 0.5), there is a corresponding probability for the complement hypothesis that indicates 
where the desired result(s) is unlikely to be. Thus, applying this reasoning to a wireless location 
application utilizing the novel method and/or system of the present disclosure, then for an hypothesis H 
indicating that the target MS is in a geographical area A, there is a dual location estimate H c that may be 
1 5 generated, wherein the H c represents the area outside of A and the probability that the target MS is 

outside of A is 1-P. Thus, with each location hypothesis having a probability favorably indicating where a 
target MS may be (i.e., P >= 0.5), there is a corresponding probability for the complement area not 
represented by the location hypothesis that does not favor the target MS being in this complement area. 
Further, note that similar dual hypotheses can be used in other applications using the multiple model 
20 architecture of the present disclosure when probabilities are assigned to hypotheses generated by the 
models of the application. 

Referring to (4.1 .4) as it relates to a wireless location system provided by the novel method 
and/or system of the present disclosure, note that, it is an aspect of the present disclosure to provide 
location hypothesis enhancing and evaluation techniques that can adjust target MS location estimates 
25 according to historical MS location data and/or adjust the confidence values of location hypotheses 

according to how consistent the corresponding target MS location estimate is: (a) with historical MS signal 
characteristic values, (b) with various physical constraints, and (c) with various heuristics. In particular, 
the following capabilities are provided by the novel method and/or system of the present disclosure: 

(5.1) A capability for enhancing the accuracy of an initial location hypothesis, H, generated by a 
30 first order model, FOM H) by using H as, essentially, a query or index into an historical data 

base (denoted herein as the location signature data base), wherein this data base 
includes: (a) a plurality of previously obtained location signature clusters (i.e., composite 
wireless signal characteristic values) such that for each such cluster there is an 
associated actual or verified MS locations where an MS communicated with the base 
35 station infrastructure for locating the MS, and (b) previous MS location hypothesis 

estimates from FOMh derived from each of the location signature clusters stored 
according to (a); 
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(5.2) A capability for analyzing composite signal characteristic values of wireless 
communications between the target MS and the base station infrastructure, wherein such 
values are compared with composite signal characteristics values of known MS locations 
(these latter values being archived in the location signature data base). In one instance, 
the composite signal characteristic values used to generate various location hypotheses for 
the target MS are compared against wireless signal data of known MS locations stored in 
the location signature data base for determining the reliability of the location hypothesizing 
models for particular geographic areas and/or environmental conditions; 

(5.3) A capability for reasoning about the likeliness of a location hypothesis wherein this 
reasoning capability uses heuristics and constraints based on physics and physical 
properties of the location geography; 

(5.4) An hypothesis generating capability for generating new location hypotheses from previous 
hypotheses. 

Note, (5.5) and (5.6) following are alternative descriptions of various inventive aspects of the 
novel method and/or system of the present disclosure: 

(5.5) A method for determining from a plurality of conditions a particular unknown condition, 
wherein for substantially every one of said conditions there is a corresponding set of 
measurements for identifying the condition, comprising: 

classifying or partitioning said plurality of conditions into a collection of classes, wherein 
for each said class, and each of said conditions therein, said corresponding set of 
measurements for identifying the condition is related to said corresponding set of 
measurements for identifying another one of said conditions in said class; 

determining a plurality of estimators for estimating said conditions when supplied with 
said corresponding set of measurements; 

storing a plurality of data item collections, wherein for each of a plurality of said 
conditions, there is one of said data item collections having: 
(a1) a representation of the condition, and 

(a2) a representation of a set of measurements for identifying said condition of 
(al); 

for each of said estimators, and each of class of a plurality of said classes, perform 
(b1)and(b2): 

(b1 ) inputting said representation of (a2) of each of a plurality of said data item 
collections to the estimator for generating a corresponding estimated 
condition, wherein said representation (a2) is contained in the class; 

(b2) comparing, for each of said data item collections, said representation (a1) with 
said corresponding estimate for determining a corresponding performance 
measurement of the estimator for the class; 
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activating a first and a second of said estimators with a particular set of 
measurements for estimating, respectively, a first and a second estimate of said particular 
unknown condition; 

finding a particular one of said classes containing said particular set of 
measurements; 

determining, for said first estimate, a first of said performance measurements for the 
first estimator using an identification of said particular class ; 

determining, for said second estimate, a second of said performance measurements 
for the second estimator using an identification of said particular class; 

obtaining a resulting estimate for the unknown condition using each of said first and 
second estimates, and said first and second performance measurements. 
(5.6) A method for determining from a plurality of conditions a particular unknown condition, 
wherein for substantially every one of said conditions there is a corresponding set of 
measurements for identifying the condition, comprising: 

classifying or partitioning said plurality of conditions into a collection of classes, wherein 
for each said class, and each of said conditions therein, said corresponding set of 
measurements for identifying the condition is related to said corresponding set of 
measurements for identifying another one of said conditions in said class; 

determining a plurality of estimators for estimating said conditions when supplied with 
said corresponding set of measurements; 

storing a plurality of data item collections, wherein for each of said estimators and 
each of a plurality of said conditions, there is one of said data item collections having: 
(a1) a representation of the condition, and 

(a2) a representation of a set of measurements for identifying said condition of 
(a1); and 

(a3) an estimate of said condition generated by said estimator when said 
representation of (a2) is input to said estimator; 
activating a first of said estimators with a particular set of measurements for 
determining a first estimate of said particular unknown condition; 

retrieving one or more of said data items wherein for each of said retrieved data items, 
said estimate (a3) has a desired relationship to said first estimate, and said estimate (a3) 
was generated by said first estimator; 

determining a second estimate of said particular unknown condition using said 
representations of (at) from said retrieved data items. 

As also mentioned above in (2.3), the novel method and/or system of the present disclosure may 
utilize adaptive signal processing techniques. One particularly important utilization of such techniques 

12 



includes the automatic tuning of the novel method and/or system of the present disclosure so that, e.g., 
such tuning can be applied to adjusting the values of location processing system parameters that affect 
the processing performed by the method and/or system of the present disclosure. For example, such 
system parameters as those used for determining the size of a geographical area to be specified when 
5 retrieving location signal data of known MS locations from the historical (location signature) data base can 
substantially affect the location processing. In particular, a system parameter specifying a minimum size 
for such a geographical area may, if too large, cause unnecessary inaccuracies in locating an MS. 
Accordingly, to accomplish a tuning of such system parameters, an adaptation engine is included in the 
novel system disclosed herein for automatically adjusting or tuning parameters used by the system 
10 disclosed herein. Note that in one embodiment, the adaptation engine is based on genetic algorithm 
techniques. 

A novel aspect of the present method and/or system disclosed herein relies on the discovery that 
in many areas where MS location services are desired, the wireless signal measurements obtained from 
communications between the target MS and the base station infrastructure are extensive enough to 

15 provide sufficiently unique or peculiar values so that the pattern of values alone may identify the location 
of the target MS. Further, assuming a sufficient amount of such location identifying pattern information is 
captured in the composite wireless signal characteristic values for a target MS, and that there is a 
technique for matching such wireless signal patterns to geographical locations, then a FOM based on this 
technique may generate a reasonably accurate target MS location estimate. Moreover, if the novel 

20 method and/or system disclosed herein (e.g., the location signature data base) has captured sufficient 
wireless signal data from location communications between MSs and the base station infrastructure 
wherein the locations of the MSs are also verified and captured, then this captured data (e.g., location 
signatures) can be used to train or calibrate such models to associate the location of a target MS with the 
distinctive signal characteristics between the target MS and one or more base stations. Accordingly, in 

25 one embodiment, the novel method and/or system disclosed herein includes one or more FOMs that may 
be generally denoted as classification models wherein such FOMs are trained or calibrated to associate 
particular composite wireless signal characteristic values with a geographical location where a target MS 
could likely generate the wireless signal samples from which the composite wireless signal characteristic 
values are derived. Further, the novel system and/ or method disclosed herein may include the capability 

30 for training and retraining such classification FOMs to automatically maintain the accuracy of these 

models even though substantial changes to the radio coverage area may occur, such as the construction 
of a new high rise building or seasonal variations (due to, for example, foliage variations). As used 
herein, "training" refers to iteratively presenting "training data" to a computational module for changing the 
behavior of the module so that the module may perform progressively better as it learns appropriate 

35 behavioral responses to the training data. Accordingly, training may include, for example, the repeated 
input of training data to an artificial neural network, or repeated statistical regression analyses on different 
and/or enhanced training data (e.g., statistical sample data sets). 
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Note that such classification FOMs, that are trained to identify target MS locations by the 
wireless signal patterns produced, constitute a particularly novel aspect of the present disclosure. It is 
well known in the wireless telephony art that the phenomenon of signal multipath and shadow fading 
renders most analytical location computational techniques such as time-of-arrival (TOA) or time- 
5 difference-of-arrival (TDOA) substantially useless in urban areas and particularly in dense urban areas. 
However, this same multipath phenomenon also may produce substantially distinct or peculiar signal 
measurement patterns, wherein such a pattern coincides with a relatively small geographical area. Thus, 
in one embodiment, the novel method and/or system disclosed herein utilizes multipath as an advantage 
for increasing accuracy where for previous location systems multipath has been a source of substantial 

10 inaccuracies. Moreover, it is worthwhile to note that the utilization of classification FOMs in high multipath 
environments is especially advantageous in that high multipath environments are typically densely 
populated. Thus, since such environments are also capable of yielding a greater density of MS location 
signal data from MSs whose actual locations can be obtained, there can be a substantial amount of 
training or calibration data captured by the novel method and/or system disclosed herein for training or 

1 5 calibrating such classification FOMs and for progressively improving the MS location accuracy of such 
models. 

Moreover, it is also an aspect of the novel method and/or system disclosed herein that 
classification FOMs may be utilized that determine target MS locations by correlating and/or associating 
network anomalous behavior with geographic locations where such behavior occurs. That is, network 

20 behaviors that are problematic for voice and/or data communication may be used advantageously for 
locating a target MS. For example, it is well known that wireless networks typically have within their 
coverage areas persistent subareas where voice quality is problematic due to, e.g., measurements 
related to high total errors, a high error rate, or change in error rate. In particular, such measurements 
may be related to frame error rates, redundancy errors, co-channel interference, excessive handoffe 

25 between base stations, and/or other call quality measurements. Additionally, measurements may be 

used that are related to subareas where wireless communication between the network and a target MS is 
not sufficient to maintain a call (i.e., "deadzones"). Thus, information about such so called problematic 
behaviors may used by, e.g., a location estimator (FOM) to generate a more accurate estimate of a target 
MS. For example, such network behavioral measurements may be provided for training an artificial 

30 neural network and/or for providing to a statistical regression analysis technique and/or statistical 
prediction models (e.g., using principle decomposition, partial least squares, or other regression 
techniques for associating or correlating such measurements with the geographic area for which they 
likely derive. Moreover, note that such network behavioral measurements can also be used to reduce the 
likelihood of a target MS being in an area if such measurements are not what would be expected for the 

35 area. 

It is also a related aspect of the present novel method and/or system disclosed herein to include 
a plurality stationary, low cost, low power location detection base stations" (LBS), each such LBS having 
both restricted range MS detection capabilities, and a built-in MS. Accordingly, a grid of such LBSs can 



be utilized for providing wireless signaling characteristic data (from their built-in MSs) for: (a) (re)training 
such classification FOMs, and (b) calibrating the FOMs so that each generated location hypothesis has a 
reliable confidence value (probability) indicative of the likeliness of the target MS being in an area 
represented by the location hypothesis. 
5 It is a further aspect of the novel method and/or system disclosed herein that the personal 

communication system (PCS) infrastructures currently being developed by telecommunication providers 
offer an appropriate localized infrastructure base upon which to build various personal location systems 
(PLS) employing embodiments of the novel method and/or system disclosed herein, and/or utilizing the 
techniques disclosed herein. In particular, the novel method and/or system is especially suitable for the 

10 location of people and/or objects using code division multiple access (CDMA) wireless infrastructures, 
although other wireless infrastructures, such as, time division multiple access (TDMA) infrastructures and 
GSM are also contemplated. Note that CDMA personal communications systems are described in the 
Telephone Industries Association standard IS-95, for frequencies below 1 GHz, and in the Wideband 
Spread- Spectrum Digital Cellular System Dual-Mode Mobile Station-Base Station Compatibility 

1 5 Standard, for frequencies in the 1 .8-1 .9 GHz frequency bands, both of which are incorporated herein by 
reference. Furthermore, CDMA general principles have also been described, for example, in U. S. Patent 
5,109,390, to Gilhousen, et al, filed November 7, 1989 and CDMA Network Engineering Handbook by 
Qualcomm, Inc., each of which is also incorporated herein by reference. 

As mentioned in the discussion of classification FOMs above, embodiments of the novel method 

20 and/or system disclosed herein can substantially automatically retrain itself to compensate for variations 
in wireless signal characteristics (e.g., multipath) due to environmental and/or topographic changes to a 
geographic area serviced by such an embodiment of the present novel method and/or system. For 
example, in one embodiment, the novel method and/or system optionally includes low cost, low power 
base stations, denoted location base stations (LBS) above, providing, for example, CDMA pilot channels 

25 to a very limited area about each such LBS. The location base stations may provide limited voice traffic 
capabilities, but each is capable of gathering sufficient wireless signal characteristics from an MS within 
the location base station's range to facilitate locating the MS. Thus, by positioning the location base 
stations at known locations in a geographic region such as, for instance, on street lamp poles and road 
signs, additional MS location accuracy can be obtained. That is, due to the low power signal output by 

30 such location base stations, for there to be signaling control communication (e.g., pilot signaling and other 
control signals) between a location base station and a target MS, the MS must be relatively near the 
location base station. Additionally, for each location base station not in communication with the target 
MS, it is likely that the MS is not near to this location base station. Thus, by utilizing information received 
from both location base stations in communication with the target MS and those that are not in 

35 communication with the target MS, the novel method and/or system disclosed herein may substantially 
narrow the possible geographic areas within which the target MS is likely to be. Further, by providing 
each location base station (LBS) with a co-located stationary wireless transceiver (denoted a built-in MS 
above) having similar functionality to an MS, the following advantages are provided: 



(6.1) assuming that the co-located base station capabilities and the stationary transceiver of an LBS 
are such that the base station capabilities and the stationary transceiver communicate with one 
another, the stationary transceiver can be signaled by another component(s) of the novel 
method and/or system disclosed herein to activate or deactivate its associated base station 
capability, thereby conserving power for the LBS that operate on a restricted power such as solar 
electrical power; 

(6.2) the stationary transceiver of an LBS can be used for transferring target MS location information 
obtained by the LBS to a conventional telephony base station; 

(6.3) since the location of each LBS is known and can be used in location processing, the novel 
method and/or system disclosed herein is able to (re)train itself in geographical areas having 
such LBSs. That is, by activating each LBS stationary transceiver so that there is signal 
communication between the stationary transceiver and surrounding base stations within range, 
wireless signal characteristic values for the location of the stationary transceiver are obtained for 
each such base station. Accordingly, such characteristic values can then be associated with the 
known location of the stationary transceiver for training various of the location processing 
modules of the novel method and/or system disclosed herein, such as the classification FOMs 
discussed above. In particular, such training and/or calibrating may include: 

(i) (re)training FOMs; 

(ii) adjusting the confidence value initially assigned to a location hypothesis according to 
how accurate the generating FOM is in estimating the location of the stationary transceiver using 
data obtained from wireless signal characteristics of signals between the stationary transceiver 
and base stations with which the stationary transceiver is capable of communicating; 

(iii) automatically updating the previously mentioned historical data base (i.e., the location 
signature data base), wherein the stored signal characteristic data for each stationary transceiver 
can be used for detecting environmental and/or topographical changes (e.g., a newly built high 
rise or other structures capable of altering the multipath characteristics of a given geographical 
area); and 

(iv) tuning of the location system parameters, wherein the steps of: (a) modifying various 
system parameters and (b) testing the performance of the modified location system on verified 
mobile station location data (including the stationary transceiver signal characteristic data), these 
steps being interleaved and repeatedly performed for obtaining better system location accuracy 
within useful time constraints. 

One embodiment of the novel method and/or system disclosed herein utilizes a mobile (location) 
base station (MBS) that can be, for example, incorporated into a vehicle, such as an ambulance, police 
car, or taxi. Such a vehicle can travel to sites having a transmitting target MS, wherein such sites may be 
randomly located and the signal characteristic data from the transmitting target MS at such a location can 
consequently be archived with a verified location measurement performed at the site by the mobile 
location base station. Moreover, it is important to note that such a mobile location base station as its 
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name implies also includes base station electronics for communicating with mobile stations, though not 
necessarily in the manner of a conventional infrastructure base station. In particular, a mobile location 
base station may only monitor signal characteristics, such as MS signal strength, from a target MS 
without transmitting signals to the target MS. Alternatively, a mobile location base station can periodically 
5 be in bi-directional communication with a target MS for determining a signal time-of-arrival (or time- 
difference-of-arrival) measurement between the mobile location base station and the target MS. 
Additionally, each such mobile location base station includes components for estimating the location of 
the mobile location base station, such mobile location base station location estimates being important 
when the mobile location base station is used for locating a target MS via, for example, time-of-arrival or 
10 time-difference-of-arrival measurements as one skilled in the art will appreciate. In particular, a mobile 
location base station can include: 

(7.1) a mobile station (MS) for both communicating with other components of an embodiment of the 
novel method and/or system disclosed herein (such as a location processing center provided by 
such an embodiment); 

1 5 (7.2) a GPS receiver for determining a location of the mobile location base station; 
(7,3} a gyroscope and other dead reckoning devices; and 
(7.4) devices for operator manual entry of a mobile location base station location. 

Furthermore, a mobile location base station includes modules for integrating or reconciling 
distinct mobile location base station location estimates that, for example, can be obtained using the 
20 components and devices of (7. 1 ) through (7.4) above. That is, location estimates for the mobile location 
base station may be obtained from: GPS satellite data, mobile location base station data provided by the 
location processing center, dead reckoning data obtained from the mobile location base station vehicle 
dead reckoning devices, and location data manually input by an operator of the mobile location base 
station. 

25 The location estimating system embodiment of the novel method and/or system disclosed herein 

offers many advantages over existing location systems. The novel method and/or system , for example, 
is readily adaptable to existing wireless communication systems and can accurately locate people and/or 
objects in a cost effective manner. The embodiments of the novel method and/or system may employ a 
number of distinctly different location estimators which provide a greater degree of accuracy and reliability 

30 than is possible with existing wireless location systems. For instance, the location models provided may 
include not only the radius-radius/TOA and TDOA techniques but also adaptive artificial neural net 
techniques. 

Further, embodiments of the novel method and/or system disclosed herein can adapt to the 
topography of an area in which location service is desired. Such embodiments may also be able to adapt 
35 to environmental changes substantially as frequently as desired. Thus, such embodiments may also be 
able to take into account changes in the location topography over time without extensive manual data 
manipulation. Moreover, embodiments of the novel method and/or system can be utilized with varying 
amounts of signal measurement inputs. Thus, if a location estimate is desired in a very short time interval 
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(e.g., less than approximately one to two seconds), then the present location estimating system can be 
used with only as much signal measurement data as is possible to acquire during an initial portion of this 
time interval. Subsequently, after a greater amount of signal measurement data has been acquired, 
additional more accurate location estimates may be obtained. Note that this capability can be useful in 
the context of 91 1 emergency response in that a first quick coarse wireless mobile station location 
estimate can be used to route a 91 1 call from the mobile station to a 91 1 emergency response center that 
has responsibility for the area containing the mobile station and the 91 1 caller. Subsequently, once the 
91 1 call has been routed according to this first quick location estimate, by continuing to receive additional 
wireless signal measurements, more reliable and accurate location estimates of the mobile station can be 
obtained. 

Moreover, there are numerous additional advantages of embodiments of the novel method 

and/or system when applied to wireless location in CDMA communication systems. A location system 

embodiment of the novel method and/or system readily benefits from the distinct advantages of the 

CDMA spread spectrum scheme. Namely, these advantages include the exploitation of radio frequency 

spectral efficiency and isolation by (a) monitoring voice activity, (b) management of two-way power 

control, (c) provisioning of advanced variable-rate modems and error correcting signal encoding, (d) 

inherent resistance to fading, (e) enhanced privacy, and (f) multiple "rake" digital data receivers and 

searcher receivers for correlation of signal multipaths. 

At a more general level, it is within the scope of the novel method and/or system disclosed 
herein to utilize various novel computational paradigms such as: 

(8.1) providing a multiple hypothesis computational architecture (as illustrated best in Fig. 8 and/or Fig. 
13) wherein the hypotheses may be: 

(8.1.1) generated by modular independent hypothesizing computational models (FOMs), 
wherein the FOMs have been calibrated to thereby output confidence values (probabilities) 
related to the likelihood of correspondingly generated hypotheses being correct; 

(8.1 .2) the FOMs are embedded in the computational architecture in a manner wherein the 
architecture allows for substantial amounts of application specific processing common or generic 
to a plurality of the models to be straightforwardly incorporated into the computational 
architecture; 

(8.1 .3) the computational architecture enhances the hypotheses generated by the models both 
according to past performance of the models and according to application specific constraints 
and heuristics without requiring complex feedback loops for recalibrating one or more of the 
FOMs; 

(8.1 .4) the FOMs are relatively easily integrated into, modified and extracted from the 
computational architecture; 

(8.2) providing a computational paradigm for enhancing an initial estimated solution to a problem by 
using this initial estimated solution as, effectively, a query or index into an historical data base of 
previous solution estimates and corresponding actual solutions for deriving an enhanced solution 
estimate based on past performance of the module that generated the initial estimated solution. 
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Further, note that the present multiple model (FOM) architecture provides additional advantages 
in that large software systems may be developed more easily. For example, if for a given application 
(e.g., wireless location, or another application such as those listed hereinbelow), there are a plurality of 
computational models that may be used, but each model in and of itself is less than completely 
satisfactory, then by integrating two or more such models into the architecture of the novel method and/or 
system disclosed herein, at least the following advantages are provided: 

(a) the models themselves may be integrated into the desired application embodiment as 
substantially a "black box" thus reducing the development effort; 

(b) the models may be easily added and deleted as appropriate. Thus, if an enhanced 
version of a model and/or an entirely new model becomes available, these models can 
be straightforwardly incorporated into the application; 

(c) a generic embodiment of the architecture of the novel method and/or system disclosed 
herein may be provided wherein, e.g., messages (having a common structure 
regardless of the application) may be passed between various components of the 
embodiment, an application development system can be provided that is based on the 
architecture of the novel method and/or system which may be used in a number of 
different technical fields 

(d) the models may be activated in parallel on different computational devices. Thus, the 
architecture of the novel method and/or system lends itself to parallel computing. 

Thus, the multiple hypothesis architecture provided herein is useful in implementing solutions in 
a wide range of applications. In fact, most of the Detailed Description hereinbelow can be immediately 
translated into other application areas, as one skilled in the art of computer application architectures will 
come to appreciate. For example, the following additional applications are within the scope of the method 
and/or system disclosed herein: 

(9.1) document scanning applications for transforming physical documents in to electronic forms of the 
documents. Note that in many cases the scanning of certain documents (books, publications, 
etc.) may have a 20% character recognition error rate. Thus, the novel computation architecture 
of the novel method and/or system can be utilized by (i) providing a plurality of document 
scanning models as the first order models, (ii) building a character recognition data base for 
archiving a correspondence between characteristics of actual printed character variations and 
the intended characters (according to, for example, font types), and additionally archiving a 
correspondence of performance of each of the models on previously encountered actual printed 
character variations (note, this is analogous to the Signature Data Base of the MS location 
application described herein), and (iii) determining any generic constraints and/or heuristics that 
are desirable to be satisfied by a plurality of the models. Accordingly, by comparing outputs from 
the first order document scanning models, a determination can be made as to whether further 
processing is desirable due to, for example, discrepancies between the output of the models. If 
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further processing is desirable, then an embodiment of the multiple hypothesis architecture 
provided herein may be utilized to correct such discrepancies. Note that in comparing outputs 
from the first order document scanning models, these outputs may be compared at various 
granularities; e.g., character, sentence, paragraph or page; 
5 (9.2) diagnosis and monitoring applications such as medical diagnosis/monitoring, communication 

network diagnosis/monitoring. Note that in many cases, the domain wherein a diagnosis is to be 
performed has a canonical hierarchical order among the components within the domain. For 
example, in automobile diagnosis, the components of an auto may be hierarchically ordered 
according to ease of replacement in combination within function. Thus, within an auto's electrical 

10 system (function), there may be a fuse box, and within the fuse box there will be fuses. Thus, 

these components may be ordered as follows (highest to lowest): auto, electrical system, fuse 
box, fuses. Thus, if different diagnostic FOMs provide different hypotheses as to a problem with 
an auto, the confidence values for each component and its subcomponents maybe summed 
together to provide a likelihood value that the problem within the component. Accordingly, the 

15 lowest component having, for example, at least a minimum threshold of summed confidences 

can be selected as the most likely component for either further analysis and/or replacement. 
Note that such summed confidences may be normalized by dividing by the number of 
hypotheses generated from the same input so that the highest summed confidence is one and 
the lowest is zero. Further note that this example is merely representative of a number of 

20 different diagnosis and/or prediction applications to which the novel method and/or system is 

applicable, wherein there are components that have canonical hierarchical decompositions. For 
example, a technique similar to the auto illustration above may be provided for the diagnosis of 
computer systems, networks (LANs, WANs, Internet and telephony networks), medical diagnosis 
from, e.g., x-rays, MRIs, sonograms, etc; 

25 (9.3) robotics applications such as scene and/or object recognition. That is, various FOMs may 

process visual image input differently, and it may be that for expediency, an object is recognized 
if the summed confidence values for the object being recognized is above a certain threshold ; 

(9.4) seismic and/or geologic signal processing applications such as for locating oil and gas deposits; 

(9.5) recognition of terrestrial and/or airborne objects from satellites, wherein there may be various 
30 spectral bands monitored. 

(9.6) Additionally, note that this architecture need not have all modules co-located. In particular, it is 
an additional aspect of embodiments of the novel method and/or system disclosed herein that various 
modules can be remotely located from one another and communicate with one another via 
telecommunication transmissions such as telephony technologies and/or the Internet. Accordingly, the 

35 novel method and/or system is particularly adaptable to such distributed computing environments. For 
example, some number of the first order models may reside in remote locations and communicate their 
generated hypotheses via the Internet. 
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For instance, in weather prediction applications it is not uncommon for computational models to 
require large amounts of computational resources. Thus, such models running at various remote 
computational facilities can transfer weather prediction hypotheses (e.g., the likely path of a hurricane) to 
a site that performs hypothesis adjustments according to: (i) past performance of the each model; (ii) 
5 particular constraints and/or heuristics, and subsequently outputs a most likely estimate for a particular 
weather condition. 

In an alternative embodiment of the novel method and/or system disclosed herein, the 
processing following the generation of location hypotheses (each having an initial location estimate) by 
the first order models may be such that this processing can be provided on Internet user nodes and the 

10 first order models may reside at Internet server sites. In this configuration, an Internet user may request 
hypotheses from such remote first order models and perform the remaining processing at his/her node. 

In other embodiments of the novel method and/or system disclosed herein, a fast, albeit less 
accurate location estimate may be initially performed for very time critical location applications where 
approximate location information may be required. For example, less than 1 second response for a 

1 5 mobile station location embodiment of the novel method and/or system may be desired for 91 1 

emergency response location requests. Subsequently, once a relatively coarse location estimate has 
been provided, a more accurate most likely location estimate can be performed by repeating the location 
estimation processing a second time with, e.g., additional with measurements of wireless signals 
transmitted between a mobile station to be located and a network of base stations with which the mobile 

20 station is communicating, thus providing a second, more accurate location estimate of the mobile station. 

Additionally, note that it is within the scope of the novel method and/or system disclosed herein 
to provide one or more central location development sites that may be networked to, for example, 
geographically dispersed location centers providing location services according to the novel method 
and/or system, wherein the FOMs may be accessed, substituted, enhanced or removed dynamically via 

25 network connections (via, e.g., the Internet) with a central location development site. Thus, a small but 
rapidly growing municipality in substantially flat low density area might initially be provided with access to, 
for example, two or three FOMs for generating location hypotheses in the municipality's relatively 
uncluttered radio signaling environment. However, as the population density increases and the radio 
signaling environment becomes cluttered by, for example, thermal noise and multipath, additional or 

30 alternative FOMs may be transferred via the network to the location center for the municipality. 

Note that in some embodiments of the novel method and/or system disclosed herein, since there 
is a lack of sequencing between the FOMs and subsequent processing of location hypotheses, the FOMs 
can be incorporated into an expert system, if desired. For example, each FOM may be activated from an 
antecedent of an expert system rule. Thus, the antecedent for such a rule can evaluate to TRUE if the 

35 FOM outputs a location hypothesis, and the consequent portion of such a rule may put the output location 
hypothesis on a list of location hypotheses occurring in a particular time window for subsequent 
processing by the location center. Alternatively, activation of the FOMs may be in the consequents of 
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such expert system rules. That is, the antecedent of such an expert system rule may determine if the 

conditions are appropriate for invoking the FOM(s) in the rule's consequent. 

The novel method and/or system may also be configured as a blackboard system with intelligent 

agents (e.g., FOMs). In this embodiment, each of the intelligent agents is calibrated using archived data 
5 so that for each of the input data sets provided either directly to the intelligent agents or to the blackboard, 

each hypothesis generated and placed on the blackboard by the intelligent agents has a corresponding 

confidence value indicative of an expected validity of the hypothesis. 

Of course, other software architectures may also to used in implementing the processing of the 

location center without departing from scope of the novel method and/or system. In particular, object- 
l o oriented architectures are also within the scope of the novel method and/or system. For example, the 

FOMs may be object methods on an MS location estimator object, wherein the estimator object receives 

substantially all target MS location signal data output by the signal filtering subsystem. Alternatively, 

software bus architectures are contemplated by the novel method and/or system, as one skilled in the art 

will understand, wherein the software architecture may be modular and facilitate parallel processing. 
15 Further features and advantages of the novel method and/or system are provided by the figures 

and detailed description accompanying this summary. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates various perspectives of radio propagation opportunities which may be 
considered in addressing correlation with mobile to base station ranging. 
20 Fig. 2 shows aspects of the two-ray radio propagation model and the effects of urban clutter. 

Fig. 3 provides a typical example of how the statistical power budget is calculated in design of a 
Commercial Mobile Radio Service Provider (CMRS) network. 

Fig. 4 illustrates an overall view of a wireless radio location network architecture, based on AIN 
principles. 

25 Fig. 5 is a high level block diagram of an embodiment of the novel method and/or system for 

locating a mobile station (MS) within a radio coverage area for the novel method and/or system. 

Fig. 6 is a high level block diagram of an embodiment of the location center 142. 

Fig. 7 is a high level block diagram of one embodiment of the hypothesis evaluator for the 
location center. 

30 Fig. 8 is a substantially comprehensive high level block diagram illustrating data and control 

flows between the components of an embodiment of the location center 142, as well the functionality of 
the components. 

Figs. 9A and 9B are a high level data structure diagram describing the fields of a location 
hypothesis object generated by the first order models 1224 of the location center. 
35 Fig. 10 is a graphical illustration of the computation performed by the most likelihood estimator 

1344 of the hypothesis evaluator. 

Fig. 11 is a high level block diagram of the mobile base station (MBS). 



Fig. 12 is a high level state transition diagram describing computational states the Mobile Base 
station enters during operation. 

Fig. 13 is a high level diagram illustrating the data structural organization of the Mobile Base 
station capability for autonomously determining a most likely MBS location from a plurality of potentially 
5 conflicting MBS location estimating sources. 

Fig. 14 illustrates the primary components of the signal processing subsystem. 

Fig. 15 illustrates another embodiment of the location engine 139, wherein the context adjuster 
1326 (denoted in this figure as location hypothesis adjuster modules") includes a module (1436) that is 
capable of adjusting location hypotheses for reliability, and another module (1440) that is capable of 
10 adjusting location hypotheses for accuracy. 

DETAILED DESCRIPTION 

Detailed Description Introduction 

Fig. 4 is a high level diagram of a wireless digital radio location intelligent network architecture 
for the novel method and/or system. Accordingly, this figure illustrates the interconnections between the 
15 components, for example, of a typical PCS network configuration and various components that are 

specific to wireless location embodiments of the novel method and/or system. In particular, as one skilled 
in the art will understand, a typical wireless (PCS) network includes: 

(a) a (large) plurality of conventional wireless mobile stations (MSs) 140 for at least one of voice 

related communication, visual (e.g., text) related communication, and (according to the 
20 novel method and/or system ) location related communication; 

(b) a mobile switching center (MSC) 112; 

(c) a plurality of wireless cell sites in a radio coverage area 120, wherein each cell site includes 

an infrastructure base station such as those labeled 1 22 (or variations thereof such as 1 22A^ 
- 122D). In particular, the base stations 122 denote the standard high traffic, fixed location 

25 base stations used for voice and data communication with a plurality of MSs 140, and, 

according to the novel method and/or system disclosed herein, also used for communication 
of information related to locating such MSs 140. Additionally, note that the base stations 
labeled 152 are more directly related to wireless location enablement. For example, as 
described in greater detail hereinbelow, the base stations 152 may be low cost, low 

30 functionality transponders that are used primarily in communicating MS location related 

information to the location center 142 (via base stations 1 22 and the MSC 1 1 2). Note that 
unless stated otherwise, the base stations 152 will be referred to hereinafter as "location 
base station(s) 15Z or simply TBS(s) 152"); 

(d) a public switched telephone network (PSTN) 124 (which may include signaling system links 
35 106 having network control components such as: a service control point (SCP) 104 , one or 

more signaling transfer points (STPs) 110. 
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Added to this wireless network, an embodiment of the novel method and/or system disclosed 
herein provides the following additional components: 

(1 0.1 ) a location center 142 which is required for determining a location of a target MS 140 using 
signal characteristic values for this target MS; 
5 - (10.2) one or more mobile base stations 148 (MBS) which are optional, for physically traveling toward 
the target MS 140 or tracking the target MS; 

(10.3) a plurality of location base stations 152 (LBS) which are optional, distributed within the radio 
coverage areas 120, each LBS 152 having a relatively small MS 140 detection area 154; 

Since location base stations can be located on potentially each floor of a multi-story building, the 

10 wireless location technology described herein can be used to perform location in terms of height as well 
as by latitude and longitude. 

In operation, the MS 140 may utilize one of the wireless technologies, CDMA, TDMA, AMPS, 
NAMPS or GSM techniques for radio communication with: (a) one or more infrastructure base stations 
122, (b) mobile base station(s) 148 , (c) an LBS 152. 

15 Referring to Fig. 4 again, additional detail is provided of typical base station coverage areas, 

sectorization, and high level components within a radio coverage area 120, including the MSC 1 12. 
Although base stations may be placed in any configuration, a typical deployment configuration is 
approximately in a cellular honeycomb pattern, although many practical tradeoffs exist, such as site 
availability, versus the requirement for maximal terrain coverage area. To illustrate, three such exemplary 

20 base stations (BSs) are 1 22A, 1 22B and 1 22C, each of which radiate referencing signals within their area 
of coverage 169 to facilitate mobile station (MS) 140 radio frequency connectivity, and various timing and 
synchronization functions. Note that some base stations may contain no sectors 130 (e.g. 122E), thus 
radiating and receiving signals in a 360 degree omnidirectional coverage area pattern, or the base station 
may contain "smart antennas" which have specialized coverage area patterns. However, the generally 

25 most frequent base stations 1 22 have three sector 1 30 coverage area patterns. For example, base 
station 122A includes sectors 130, additionally labeled a, b and c. Accordingly, each of the sectors 130 
radiate and receive signals in an approximate 120 degree arc, from an overhead view. As one skilled in 
the art will understand, actual base station coverage areas 169 (stylistically represented by hexagons 
about the base stations 122) generally are designed to overlap to some extent, thus ensuring seamless 

30 coverage in a geographical area. Control electronics within each base station 122 are used to 

communicate with a mobile stations 140. Information regarding the coverage area for each sector 130, 
such as its range, area, and "holes" or areas of no coverage (within the radio coverage area 120), may be 
known and used by the location center 142 to facilitate location determination. Further, during 
communication with a mobile station 1 40, the identification of each base station 1 22 communicating with 

35 the MS 140 as well, as any sector identification information, may be known and provided to the location 
center 142. 

In the case of the base station types 122, 148, and 152 communication of location information, a 
base station or mobility controller 1 74 (BSC) controls, processes and provides an interface between 



originating and terminating telephone calls from/to mobile station (MS) 140, and the mobile switch center 
(MSC) 112. The MSC 122, on-the-other-hand, performs various administration functions such as mobile 
station 140 registration, authentication and the relaying of various system parameters, as one skilled in 
the art will understand. 

5 The base stations 1 22 may be coupled by various transport facilities 1 76 such as leased lines, frame 

relay, T-Carrier links, optical fiber links or by microwave communication links. 

When a mobile station 140 (such as a CDMA, AMPS, NAMPS mobile telephone) is powered on and 
in the idle state, it constantly monitors the pilot signal transmissions from each of the base stations 122 
located at nearby cell sites. Since base station/sector coverage areas may often overlap, such 

l o overlapping enables mobile stations 140 to detect, and, in the case of certain wireless technologies, 
communicate simultaneously along both the forward and reverse paths, with multiple base stations 122 
and/or sectors 130. In Fig. 4 the constantly radiating pilot signals from base station sectors 130, such as 
sectors a, b and c of BS 122A, are detectable by mobile stations 140 within the coverage area 169 for BS 
122A. That is, the mobile stations 140 scan for pilot channels, corresponding to a given base 

15 station/sector identifiers (IDs) , for determining which coverage area 169 (i.e., cell) it is contained. This is 
performed by comparing signals strengths of pilot signals transmitted from these particular cell-sites. 

The mobile station 140 then initiates a registration request with the MSC 1 12, via the base 
station controller 174. The MSC 112 determines whether or not the mobile station 140 is allowed to 
proceed with the registration process (except in the case of a 91 1 call, wherein no registration process is 

20 required). At this point calls may be originated from the mobile station 140 or calls or short message 

service messages can be received from the network. The MSC 1 12 communicates as appropriate, with a 
class 4/5 wireline telephony circuit switch or other central offices, connected to the PSTN 124 network. 
Such central offices connect to wireline terminals, such as telephones, or any communication device 
compatible with the line. The PSTN 124 may also provide connections to long distance networks and 

25 other networks. 

The MSC 112 may also utilize IS/41 data circuits or trunks connecting to signal transfer point 
110, which in turn connects to a service control point 104, via Signaling System #7 (SS7) signaling links 
(e.g., trunks) for intelligent call processing, as one skilled in the art will understand. In the case of 
wireless AIN services such links are used for call routing instructions of calls interacting with the MSC 

30 1 1 2 or any switch capable of providing service switching point functions, and the public switched 
telephone network (PSTN) 124, with possible termination back to the wireless network. 

Referring to Fig. 4 again, the location center (LC) 142 interfaces with the MSC 112 either via 
dedicated transport facilities 178, using for example, any number of LAN/WAN technologies, such as 
Ethernet, fast Ethernet, frame relay, virtual private networks, etc., or via the PSTN 124. The LC 142 

35 receives autonomous (e.g., unsolicited) command/response messages regarding, for example: (a) the 
state of the wireless network of each service provider, (b) MS 140 and BS 122 radio frequency (RF) 
measurements, (c) any MBSs 148, (d) location applications requesting MS locations using the location 
center. Conversely, the LC 142 provides data and control information to each of the above components 



in (a) - (d). Additionally, the LC 142 may provide location information to an MS 140, via a BS 122. 
Moreover, in the case of the use of a mobile base station (MBS) 148, several communications paths may 
exist with the LC 142. 

The MBS 148 acts as a low cost, partially-functional, moving base station, and is, in one 
5 embodiment, situated in a vehicle where an operator may engage in MS 140 searching and tracking 

activities. In providing these activities using CDMA, the MBS 148 provides a forward link pilot channel for 
a target MS 140, and subsequently receives unique BS pilot strength measurements from the MS 140. 
The MBS 148 also includes a mobile station for data communication with the LC 142, via a BS 122. In 
particular, such data communication includes telemetering the geographic position of the MBS 148 as 

10 well as various RF measurements related to signals received from the target MS 140. In some 
embodiments, the MBS 148 may also utilize multiple-beam fixed antenna array elements and/or a 
moveable narrow beam antenna , such as a microwave dish 182. The antennas for such embodiments 
may have a known orientation in order to further deduce a radio location of the target MS 140 with 
respect to an estimated current location of the MBS 148. As will be described in more detail herein 

15 below, the MBS 148 may further contain a global positioning system (GPS), distance sensors, dead- 
reckoning electronics, as well as an on-board computing system and display devices for locating both the 
MBS 148 itself as well as tracking and locating the target MS 140. The computing and display provides a 
means for communicating the position of the target MS 140 on a map display to an operator of the MBS 
148. 

20 Each location base station (LBS) 1 52 is a low cost location device. Each such LBS 1 52 

communicates with one or more of the infrastructure base stations 122 using one or more wireless 
technology interface standards. In some embodiments, to provide such LBS's cost effectively, each LBS 
152 only partially or minimally supports the air-interface standards of the one or more wireless 
technologies used in communicating with both the BSs 122 and the MSs 140. Each LBS 152, when put 

25 in service, is placed at a fixed location, such as at a traffic signal, lamp post, etc., and wherein the location 
of the LBS may be determined as accurately as, for example, the accuracy of the locations of the 
infrastructure BSs 122. Assuming the wireless technology CDMA is used, each BS 1 22 uses a time 
offset of the pilot PN sequence to identify a forward CDMA pilot channel. In one embodiment, each LBS 
152 emits a unique, time-offset pilot PN sequence channel in accordance with the CDMA standard in the 

30 RF spectrum designated for BSs 1 22, such that the channel does not interfere with neighboring BSs 122 
cell site channels, nor would it interfere with neighboring LBSs 152. However, as one skilled in the art will 
understand, time offsets, in CDMA chip sizes, may be re-used within a PCS system, thus providing 
efficient use of pilot time offset chips, thereby achieving spectrum efficiency. Each LBS 152 may also 
contain multiple wireless receivers in order to monitor transmissions from a target MS 140. Additionally, 

35 each LBS 1 52 contains mobile station 140 electronics, thereby allowing the LBS to both be controlled by 
the LC 142, and to transmit information to the LC 142, via at least one neighboring BS 122. 

As mentioned above, when the location of a particular target MS 140 is desired, the LC 142 can 
request location information about the target MS 140 from, for instance, one or more activated LBSs 152 
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in a geographical area of interest. Accordingly, whenever the target MS 140 is in such an area, or is 
suspected of being in the area, either upon command from the LC 142, or in a substantially continuous 
fashion, the LBS's pilot channel appears to the target MS 140 as a potential neighboring base station 
channel, and consequently, is placed, for example, in the CDMA neighboring set, or the CDMA remaining 
5 set, of the target MS 1 40 (as one familiar with the CDMA standards will understand). 

During the normal CDMA pilot search sequence of the mobile station initialization state (in the 
target MS), the target MS 140 will, if within range of such an activated LBS 152, detect the LBS pilot 
presence during the CDMA pilot channel acquisition substate. Consequently, the target MS 140 performs 
RF measurements on the signal from each detected LBS 152. Similarly, an activated LBS 152 can 

10 perform RF measurements on the wireless signals from the target MS 140. Accordingly, each LBS 152 
detecting the target MS 140 may subsequently telemeter back to the LC 142 measurement results related 
to signals from/to the target MS 140. Moreover, upon command, the target MS 140 will telemeter back to 
the LC 142 its own measurements of the detected LBSs 152, and consequently, this new location 
information, in conjunction with location related information received from the BSs 122, can be used to 

1 5 locate the target MS 1 40. 

It should be noted that an LBS 152 will normally deny hand-off requests, since typically the LBS 
does not require the added complexity of handling voice or traffic bearer channels, although economics 
and peak traffic load conditions would dictate preference here. GPS timing information, needed by any 
CDMA base station, is either achieved via a the inclusion of a local GPS receiver or via a telemetry 

20 process from a neighboring conventional BS 122, which contains a GPS receiver and timing information. 
Since energy requirements are minimal in such an LBS 152, (rechargeable) batteries or solar cells may 
be used to power the LBS. No expensive terrestrial transport link is typically required since two-way 
communication is provided by the included MS 140 (or an electronic variation thereof). Thus, LBSs 152 
may be placed in numerous locations, such as: 

25 (a) in dense urban canyon areas (e.g., where signal reception may be poor and/or very noisy); 

(b) in remote areas (e.g., hiking, camping and skiing areas); 

(c) along highways (e.g., for emergency as well as monitoring traffic flow), and their rest stations; 

or 

(d) in general, wherever more location precision is required than is obtainable using other 
30 wireless infrastructure network components. 

Location Center - Network Elements API Description 

A location application programming interface, or L-AP1 14 (see Fig. 14, and including L-API- 
Loc.APP 135, L-API-MSC 136, and L-API-SCP 137 shown in Fig. 4), is required between the location 
center 142 (LC) and the mobile switch center (MSC) network element type, in order to send and receive 
35 various control, signals and data messages. The L-AP1 14 should be implemented using a preferably 

high-capacity physical layer communications interface, such as IEEE standard 802.3 (10 baseT Ethernet), 
although other physical layer interfaces could be used, such as fiber optic ATM, frame relay, etc. Two 



forms of API implementation are possible. In the first case the signals control and data messages are 
realized using the MSC 112 vendor's native operations messages inherent in the product offering, without 
any special modifications. In the second case the L-API includes a full suite of commands and 
messaging content specifically optimized for wireless location purposes, which may require some, 
5 although minor development on the part of the MSC vendor. 

Signal Processor Description 

Referring to Fig. 14, the signal processing subsystem 1220 receives control messages and 
signal measurements and transmits appropriate control messages to the wireless network via the location 
applications programming interface referenced earlier, for wireless location purposes. The signal 

10 processing subsystem additionally provides various signal identification, conditioning and pre-processing 
functions, including buffering, signal type classification, signal filtering, message control and routing 
functions to the location estimate modules. 

One adjustment variable and one factor value are required by the signal processing subsystem 
in the CDMA air interface case: 1 .) instantaneous relative power level in dBm (IRPL) of the mobile station 

1 5 transmitter, and 2.) the mobile station Power Class. By adding the IRPL to the RRSSms, a synthetic 

relative signal strength (SRSSms) of the mobile station 140 signal detected at the BS 122 is derived, which 
can be used by location estimate model analysis, as shown below: 

SRSSms = RRSSms + IRPL (in dBm) 

SRSSms, a corrected indication of the effective path loss in the reverse direction (mobile station to BS), is 

20 now comparable with RRSSbs and can be used to provide a correlation with either distance or shadow 
fading because it now accounts for the change of the mobile station transmitter's power level. The two 
signals RRSSbs and SRSSms can now be processed in a variety of ways to achieve a more robust 
correlation with distance or shadow fading. 

The signal processing subsystem 1220 provides the low level processing wireless signal 

25 measurements received from a target MS 140 via the network of base stations 122. In particular, the 
signal processing subsystem 1220 aggregates such measurements for decreasing the volatility of such 
measurements. Moreover, the signal processing subsystem 1220 outputs ensembles of data that are 
suitable for processing by the FOMs1224. A more detailed description of the signal processing 
subsystem 1220 can be found in the publication, "Location Of A Mobile Station Using A Plurality Of 

30 Commercial Wireless Infrastructures," by LeBlanc, Dupray and Karr having International Patent 

Application No. PCT/US97/15933 filed September 8, 1997, and also having a corresponding U.S. national 
Patent Application No. 09/230,109 filed January 22, 1999, now U.S. Patent No. 6,236,365, wherein these 
documents are fully incorporated herein by reference. 

LOCATION CENTER HIGH LEVEL FUNCTIONALITY 

35 At a very high level the location center 142 computes location estimates for a wireless Mobile 

Station 140 (denoted the "target MS" or "MS") by performing the following steps: 
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(23.1 ) receiving signal transmission characteristics of communications communicated between the 
target MS 140 and one or more wireless infrastructure base stations 122; 

(23.2) filtering the received signal transmission characteristics (by the signal processing subsystem 
1220 illustrated in Fig. 5) as needed so that target MS location data can be generated that is uniform and 

5 consistent with location data generated from other target MSs 140. In particular, such uniformity and 
consistency is both in terms of data structures and interpretation of signal characteristic values provided 
by the MS location data; 

(23.3) inputting the generated target MS location data to one or more MS location estimating models 
(denoted First order models or FOMs, and labeled collectively as 1224 in Fig. 5), so that each such model 

10 may use the input target MS location data for generating a "location hypothesis" providing an estimate of 
the location of the target MS 140; 

(23.4) providing the generated location hypotheses to an hypothesis evaluation module (denoted the 
hypothesis evaluator 1 228 in Fig. 5); 

(a) for adjusting at least one of the target MS location estimates of the generated location 
15 hypotheses and related confidence values indicating the confidence given to each location estimate, 

wherein such adjusting uses archival information related to the accuracy of previously generated location 
hypotheses, 

(b) for evaluating the location hypotheses according to various heuristics related to, for 
example, the radio coverage area 120 terrain, the laws of physics, characteristics of likely movement of 

20 the target MS 140; and 

(c) for determining a most likely location area for the target MS 140, wherein the 
measurement of confidence associated with each input MS location area estimate is used for 
determining a "most likely location area"; and 

(23.5) outputting a most likely target MS location estimate to one or more applications 146 (Fig. 5) 
25 requesting an estimate of the location of the target MS 140. 

Location Hypothesis Data Representation 

In order to describe how the steps (23.1 ) through (23.5) are performed in the sections below, 
some introductory remarks related to the data denoted above as location hypotheses will be helpful. 
Additionally, it will also be helpful to provide introductory remarks related to historical location data and 
30 the data base management programs associated therewith. 

For each target MS location estimate generated and utilized by a wireless location embodiment 
of the novel method and/or system disclosed herein, the location estimate is provided in a data structure 
(or object class) denoted as a "location hypothesis" (illustrated in Table LH-1 ). Brief descriptions of the 
data fields for a location hypothesis is provided in the Table LH-1 . 

35 
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Table LH-1 



FOMJD 


First order model ID (providing this Location Hypothesis); note, 
since it is possible for location hypotheses to be generated by other 

than thp FOMq 1994 in npnpral thte fiplrl irlpntifip^ ihp moHnlp that 

generated this location hypothesis. 


MSJD 


The identification of the target MS 140 to this location hypothesis 


pLest 


The most likely location point estimate of the target MS 140. 


valicLpt 


Boolean indicating the validity of u pt_esf. 


area_est 


Location Area Estimate of the target MS 140 provided by the FOM. 
This area estimate will be used whenever "image.area" below is 
NULL 


valid_area 


Boolean indicating the validity of "area_esf (one of "pLest" and 
°area„est B must be valid). 


adjust 


Boolean (true if adjustments to the fields of this location hypothesis 
are to be performed in the Context adjuster Module). 


pt_covering 


Reference to a substantially minimal area (e.g., mesh cell) covering 
of H pt_est". Note, since this MS 140 may be substantially on a cell 
Dounaary, wis covenng may, in some cases, inciuae more man one 
cell. 


image_area 


Reference to a substantially minimal area (e.g., mesh cell) covering 
of "pt_covering" (see detailed description of the function, 
"confidence.adjuster"). Note that if this field is not NULL, then this 
is the target MS location estimate used by the location center 142 
instead of °area_esf . 


extrapolation_area 


Reference to (if non-NULL) an extrapolated MS target estimate 
area provided by the location extrapolator submodule 1432 of the 
hypothesis analyzer 1332. That is, this field, if non-NULL, is an 
extrapolation of the tt image_area B field if it exists, otherwise this field 
is an extrapolation of the "area_est" field. Note other extrapolation 
fields may also be provided depending on the embodiment of the 
novel method and/or system, such as an extrapolation of the 
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"pLcovering". 


confidence 


In one embodiment, this is a probability indicating a likelihood that 
the target MS 140 is in (or out) of a particular area. If "image.area" 
exists, then this is a measure of the likelihood that the target MS 
140 is within the area represented by "image_area B t or if 
u image_area B has not been computed (e.g., "adjust" is FALSE), 




then u area_est" must be valid and this is a measure of the likelihood 
that the target MS 140 is within the area represented by B area_esf . 
Other embodiments, are also within the scope of the novel method 
and/or system that are not probabilities; e.g., translations and/or 
expansions of the [0, 1] probability range as one skilled in the art 
will understand. 


OriginaLTimestamp 


Date and time that the location signature cluster (defined 
hereinbelow) for this location hypothesis was received by the signal 
processing subsystem 1220. 


Active_Timestamp 


Run-time field providing the time to which this location hypothesis 
has had its MS location estimate(s) extrapolated (in the location 
extrapolator 1432 of the hypothesis analyzer 1332). Note that this 
field is initialized with the value from the "OriginaLTimestamp" field. 


Processing Tags and 
environmental categorizations 


For indicating particular types of environmental classifications not 
readily determined by the "OriginaLTimestamp" field (e.g., weather, 
traffic^, and restrictions on location hypothesis processing. 


loc_sig_cluster 


Provides access to the collection of location signature signal 
characteristics derived from communications between the target MS 
140 and the base station(s) detected by this MS (discussed in detail 
hereinbelow); in particular, the location data accessed here is 
provided to the first order models by the signal processing 
subsystem 1220; i.e., access to the "loc sigs" (received at 
"timestamp" regarding the location of the target MS) 


descriptor 


Original descriptor (from the First order model indicating why/how 
the Location Area Estimate and Confidence Value were 
determined). 
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As can be seen in the Table LH-1 , each location hypothesis data structure includes at least one 
measurement, denoted hereinafter as a confidence value (or simply confidence), that is a measurement 
of the perceived likelihood that an MS location estimate in the location hypothesis is an accurate location 
estimate of the target MS 140. Since such confidence values are an important aspect of the novel 
5 method and/or system disclosed herein, the description and use of such confidence values are provided 
further below; however, a brief description is provided here. 

In one embodiment, each confidence value is a probability indicative of a likeliness that the 
target MS 140 resides within an geographic area represented by the hypothesis to which the confidence 
value applies. Accordingly, each such confidence value is in the range [0, 1]. Moreover, for clarity of 

10 discussion, it is assumed that unless stated otherwise that the probabilistic definition provided here is to 
be used when confidence values are discussed. 

Note, however, other definitions of confidence values are within the scope of the novel method 
and/or system disclosed herein, wherein such other definitions may be more general than probabilities, 
and/or that have different ranges other than [0, 1 ]. For example, one such alternative is that each such 

1 5 confidence value is in the range -1 .0 to 1 .0, wherein the larger the value, the greater the perceived 
likelihood that the target MS 140 is in (or at) a corresponding MS location estimate of the location 
hypothesis to which the confidence value applies. As an aside, note that a location hypothesis may have 
more than one MS location estimate (as will be discussed in detail below) and the confidence value will 
typically only correspond or apply to one of the MS location estimates in the location hypothesis. Further, 

20 values for the confidence value field may be interpreted as: (a) -1 .0 may be interpreted to mean that the 
target MS 140 is NOT in such a corresponding MS area estimate of the location hypothesis area, (b) 0 
may be interpreted to mean that it is unknown as to the likelihood of whether the MS 140 in the 
corresponding MS area estimate, and (c) +1 .0 may be interpreted to mean that the MS 140 is perceived 
to positively be in the corresponding MS area estimate. 

25 Additionally, note that it is within the scope of the novel method and/or system disclosed herein 

that the location hypothesis data structure may also include other related "perception" measurements 
related to a likelihood of the target MS 140 being in a particular MS location area estimate. For example, 
it is within the scope of the novel method and/or system to also utilize measurements such as, (a) 
"sufficiency factors" for indicating the likelihood that an MS location estimate of a location hypothesis is 

30 sufficient for locating the target MS 140; (b) "necessity factors" for indicating the necessity that the target 
MS be in an particular area estimate. However, to more easily describe the novel method and/or system 
disclosed herein, a single confidence field is used having the interpretation given above. 

Additionally, in utilizing location hypotheses in, for example, the location evaluator 1228 as in 
(23.4) above, it is important to keep in mind that each location hypothesis confidence value is a relative 

35 measurement. That is, for confidences, cfi and cf 2 , if cfi <= cf2, then for a location hypotheses Hi and H 2 
having cfi and cf2, respectively, the target MS 140 is expected to more likely reside in a target MS 
estimate of H 2 than a target MS estimate of Hi. Moreover, if an area, A, is such that it is included in a 
plurality of location hypothesis target MS estimates, then a confidence score, CS A , can be assigned to A, 



wherein the confidence score for such an area is a function of the confidences for all the location 
hypotheses whose (most pertinent) target MS location estimates contain A. That is, in order to determine 
a most likely target MS location area estimate for outputting from the location center 142, a confidence 
score is determined for areas within the location center service area. More particularly, if a function, T, is 
5 a function of the confidence(s) of location hypotheses, and f is a monotonic function in its parameters and 

f(cfi, cf 2) cf 3 cf N > = CS A for confidences cf, of location hypotheses H, i=1 t 2,...,N, with A contained in 

the area estimate for Hj, then T is denoted a confidence score function. Accordingly, there are many 
embodiments for a confidence score function f that may be utilized in computing confidence scores with 
the novel method and/or system disclosed herein; e.g., 
i o (a) f(cfi , cf 2 cf N ) = £ cfi = CS A ; 

(b) f(cfi, cf 2 cf N ) = I cfi" = CS A , n = 1 , 2, 3, 4, 5, 

(c) f(cfi, cf2, ... , ch) = £ (Kj * cfi) = CSa , wherein K, i = 1 , 2, ... N are positive system (tunable) 
constants (possibly dependent on environmental characteristics such as topography, time, date, traffic, 
weather, and/or the type of base station(s) 122 from which location signatures with the target MS 140 are 

15 being generated, etc.). 

For the present description of the novel method and/or system, the function f as defined in (c) 
immediately above is utilized. However, for obtaining a general understanding of the novel method and/or 
system, the simpler confidence score function of (a) may be more useful. It is important to note, though, 
that it is within the scope of the novel method and/or system to use other functions for the confidence 

20 score function. 

Coverage Area: Area Types And Their Determination 

The notion of "area type" as related to wireless signal transmission characteristics has been 
used in many investigations of radio signal transmission characteristics. Some investigators, when 
investigating such signal characteristics of areas have used somewhat naive area classifications such as 
25 urban, suburban, rural, etc. However, it is desirable for the purposes of the novel method and/or system 
to have a more operational definition of area types that is more closely associated with wireless signal 
transmission behaviors. 

To describe embodiments of the an area type scheme used in the novel method and/or system 
disclosed herein, some introductory remarks are first provided. Note that the wireless signal transmission 
30 behavior for an area depends on at least the following criteria: 

(23.8.1) substantially invariant terrain characteristics (both natural and man-made) of the area; 
e.g., mountains, buildings, lakes, highways, bridges, building density; 

(23.8.2) time varying environmental characteristics (both natural and man-made) of the area; 
e.g., foliage, traffic, weather, special events such as baseball games; 

35 (23.8.3) wireless communication components or infrastructure in the area; e.g., the arrangement 

and signal communication characteristics of the base stations 122 in the area (e.g., base 
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station antenna downtilt). Further, the antenna characteristics at the base stations 122 may 
be important criteria. 

Accordingly, a description of wireless signal characteristics for determining area types could 
potentially include a characterization of wireless signaling attributes as they relate to each of the above 
5 criteria. Thus, an area type might be: hilly, treed, suburban, having no buildings above 50 feet, with base 
stations spaced apart by two miles. However, a categorization of area types is desired that is both more 
closely tied to the wireless signaling characteristics of the area, and is capable of being computed 
substantially automatically and repeatedly over time. Moreover, for a wireless location system, the 
primary wireless signaling characteristics for categorizing areas into at least minimally similar area types 

10 are: thermal noise and, more importantly, multipath characteristics (e.g., multipath fade and time delay). 

Focusing for the moment on the multipath characteristics, it is believed that (23.8.1) and (23.8.3) 
immediately above are, in general, more important criteria for accurately locating an MS 140 than 
(23.8.2). That is, regarding (23.8.1), multipath tends to increase as the density of nearby vertical area 
changes increases. For example, multipath is particularly problematic where there is a high density of 

15 high rise buildings and/or where there are closely spaced geographic undulations. In both cases, the 
amount of change in vertical area per unit of area in a horizontal plane (for some horizontal reference 
plane) may be high. Regarding (23.8.3), the greater the density of base stations 122, the less 
problematic multipath may become in locating an MS 140. Moreover, the arrangement of the base 
stations 122 in the radio coverage area 120 in Fig. 4 may affect the amount and severity of multipath. 

20 Accordingly, it would be desirable to have a method and system for straightforwardly determining 

area type classifications related to multipath, and in particular, multipath due to (23.8.1) and (23.8.3). The 
novel method and/or system disclosed herein provides such a determination by utilizing a novel notion of 
area type, hereinafter denoted "transmission area type" (or, "area type" when both a generic area type 
classification scheme and the transmission area type discussed hereinafter are intended) for classifying 

25 "similar" areas, wherein each transmission area type class or category is intended to describe an area 
having at least minimally similar wireless signal transmission characteristics. That is, the novel 
transmission area type scheme of the novel method and/or system is based on: (a) the terrain area 
classifications; e.g., the terrain of an area surrounding a target MS 140, (b) the configuration of base 
stations 122 in the radio coverage area 120, and (c) characterizations of the wireless signal transmission 

30 paths between a target MS 140 location and the base stations 1 22. 

In one embodiment of a method and system for determining such (transmission) area type 
approximations, a partition (denoted hereinafter as Po) is imposed upon the radio coverage area 120 for 
partitioning for radio coverage area into subareas, wherein each subarea is an estimate of an area having 
included MS 140 locations that are likely to have is at least a minimal amount of similarity in their wireless 

35 signaling characteristics. To obtain the partition P 0 of the radio coverage area 1 20, the following steps are 
performed: 

(23.8.4.1) Partition the radio coverage area 120 into subareas, wherein in each subarea is: (a) 
connected, (b) the subarea is not too oblong, e.g., the variations in the lengths of chords 

34 



sectioning the subarea through the centroid of the subarea are below a predetermined 
threshold, (c) the size of the subarea is below a predetermined value, and (d) for most 
locations (e.g., within a first or second standard deviation) within the subarea whose 
wireless signaling characteristics have been verified, it is likely (e.g., within a first or 
second standard deviation ) that an MS 140 at one of these locations will detect (forward 
transmission path) and/or will be detected (reverse transmission path) by a same 
collection of base stations 122, For example, in a CDMA context, a first such collection 
may be (for the forward transmission path) the active set of base stations 122, or, the 
union of the active and candidate sets, or, the union of the active, candidate and/or 
remaining sets of base stations 122 detected by "most" MSs 140 in the subarea. 
Additionally (or alternatively), a second such collection may be the base stations 122 that 
are expected to detect MSs 140 at locations within the subarea. Of course, the union or 
intersection of the first and second collections is also within the scope of the method 
and/or system for partitioning the radio coverage area 120 according to (d) above. It is 
worth noting that it is believed that base station 122 power levels will be substantially 
constant. However, even if this is not the case, one or more collections for (d) above may 
be determined empirically and/or by computationally simulating the power output of each 
base station 122 at a predetermined level. Moreover, it is also worth mentioning that this 
step is relatively straightforward to implement using the data stored in the location 
signature data base 1320 (i.e., the verified location signature clusters discussed in detail 
hereinbelow). Denote the resulting partition here as Pi. 

(23.8.4.2) Partition the radio coverage area 120 into subareas, wherein each subarea appears to 
have substantially homogeneous terrain characteristics. Note, this may be performed 
periodically substantially automatically by scanning radio coverage area images obtained 
from aerial or satellite imaging. For example, EarthWatch Inc. of Longmont, CO can 
provide geographic with 3 meter resolution from satellite imaging data. Denote the 
resulting partition here as P2. 

(23.8.4.3) Overlay both of the above partitions, Pi and P2 of the radio coverage area 120 to 
obtain new subareas that are intersections of the subareas from each of the above 
partitions. This new partition is P 0 (i.e., Po = Pi intersect P2), and the subareas of it are 
denoted as "Po subareas". 

Now assuming Po has been obtained, the subareas of Po are provided with a first classification or 
categorization as follows: 

(23.8.4.4) Determine an area type categorization scheme for the subareas of Pi. For example, a 
subarea, A, of Pi, may be categorized or labeled according to the number of base stations 
122 in each of the collections used in (23.8.4. 1)(d) above for determining subareas of Pi. 
Thus, in one such categorization scheme, each category may correspond to a single 
number x (such as 3), wherein for a subarea, A, of this category, there is a group of x 
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(e.g., three) base stations 122 that are expected to be detected by a most target MSs 140 
in the area A. Other embodiments are also possible, such as a categorization scheme 
wherein each category may correspond to a triple: of numbers such as (5, 2, 1), wherein 
for a subarea A of this category, there is a common group of 5 base stations 1 22 with two- 
way signal detection expected with most locations (e.g., within a first or second deviation) 
within A, there are 2 base stations that are expected to be detected by a target MS 140 in 
A but these base stations can not detect the target MS, and there is one base station 122 
that is expected to be able to detect a target MS in A but not be detected. 

(23.8.4.5) Determine an area type categorization scheme for the subareas of P2. Note that the 
subareas of P2 may be categorized according to their similarities. In one embodiment, 
such categories may be somewhat similar to the naive area types mentioned above (e.g., 
dense urban, urban, suburban, rural, mountain, etc.). However, it is also an aspect of 
wireless location embodiments of the novel method and/or system disclosed herein that 
more precise categorizations may be used, such as a category for all areas having 
between 20,000 and 30,000 square feet of vertical area change per 1 1 ,000 square feet of 
horizontal area and also having a high traffic volume (such a category likely 
corresponding to a "moderately dense urban" area type). 

(23.8.4.6) Categorize subareas of P 0 with a categorization scheme denoted the tt P 0 
categorization," wherein for each Po subarea, A, a a Po area type" is determined for A 
according to the following substep(s): 

(a) Categorize A by the two categories from (23.8.4.4) and (23.8.5) with which it is 
identified. Thus, A is categorized (in a corresponding Po area type) both 
according to its terrain and the base station infrastructure configuration in the 
radio coverage area 120. 

(23.8.4.7) For each Po subarea, A, of Po perform the following step(s): 

(a) Determine a centroid, C(A), for A; 

(b) Determine an approximation to a wireless transmission path between C(A) and 
each base station 122 of a predetermined group of base stations expected to 
be in (one and/or two-way) signal communication with most target MS 140 
locations in A. For example, one such approximation is a straight line between 
C(A) and each of the base stations 1 22 in the group. However, other such 
approximations are within the scope of wireless location embodiments of the 
novel method and/or system, such as, a generally triangular shaped area as 
the transmission path, wherein a first vertex of this area is at the 
corresponding base station for the transmission path, and the sides of the 
generally triangular shaped defining the first vertex have a smallest angle 
between them that allows A to be completely between these sides. 
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5 (c) For each base station 1 22, BSs, in the group mentioned in (b) above, create an 

empty list, BSHist, and put on this list at least the P 0 area types for the 
"significant" Po subareas crossed by the transmission path between C(A) and 
BS. Note that "significant" Po subareas may be defined as, for example, the 
- 5 Po subareas through which at least a minimal length of the transmission path 

traverses. Alternatively, such "significant" P 0 subareas may be defined as 
those Po subareas that additionally are known or expected to generate 
substantial multipath. 
(d) Assign as the transmission area type for A as the collection of BSHists. Thus, 
l o any other Po subarea having the same (or substantially similar) collection of 

lists of Po area types will be viewed as having approximately the same radio 
transmission characteristics. 
Note that other transmission signal characteristics may be incorporated into the transmission 
area types. For example, thermal noise characteristics may be included by providing a third radio 
1 5 coverage area 120 partition, P3, in addition to the partitions of Pi and P2 generated in (23.8.4.1) and 
(23.8.4.2) respectively. Moreover, the time varying characteristics of (23.8.2) may be incorporated in the 
transmission area type frame work by generating multiple versions of the transmission area types such 
that the transmission area type for a given subarea of Po may change depending on the combination of 
time varying environmental characteristics to be considered in the transmission area types. For instance, 
20 to account for seasonality, four versions of the partitions Pi and P2 may be generated, one for each of the 
seasons, and subsequently generate a (potentially) different partition Po for each season. Further, the 
type and/or characteristics of base station 122 antennas may also be included in an embodiment of the 
transmission area type. 

Other embodiments of area types are also within the scope of the novel method and/or system 
25 disclosed herein. As mentioned above, each of the first order models 1224 have default confidence 
values associated therewith, and these confidence values may be probabilities. More precisely, such 
probability confidence values can be determined as follows. Assume there is a partition of the coverage 
area into subareas, each subarea being denoted a "partition area." For each partition area, activate each 
first order model 1224 with historical location data in the Location Signature Data Base 1320 (Fig. 6), 
30 wherein the historical location data has been obtained from corresponding known mobile station locations 
in the partition area. For each first order model, determine a probability of the first order model 
generating a location hypothesis whose location estimate contains the corresponding known mobile 
[ station location. To accomplish this, assume the coverage area is partitioned into partition areas A, 

wherein each partition area A is specified as the collection of coverage area locations such that for each 
35 location, the detected wireless transmissions between the network base stations and a target mobile 
station at the location can be straightforwardly equated with other locations of area A. For example, one 
such partition, Po, can be defined wherein each partition area A is specified in terms of three sets of base 
* station identifiers, namely, (a) the base station identifiers of the base stations that can be both detected at 
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each location of A and can detect a target mobile station at each location, (b) the identifiers for base 
stations that can detect a target mobile station at each location of A, but can not be detected by the target 
mobile station, and (c) the identifiers for base stations that can be detected by a target mobile station at 
each location of A, but these base stations can not detect the target mobile station. That is, two locations, 
5 h and l 2 . are identified as being in A if and only if the three sets of (a), (b), and (c) for h are, respectively, 
identical to the three sets of (a), (b), and (c) for 

Accordingly, assuming the partition P 0 as described immediately above is used, a description 
can be given as to how probabilities may be assigned as the confidence values of location hypotheses 
generated by the first order models 1224. For each partition area A, a first order model 1224 is supplied 

10 with wireless measurements of archived location data in the Location Signature Data Base associated 
with corresponding verified mobile station locations. Thus, a probability can be determined as to how 
likely the first order model is to generate a location hypothesis having a location estimate containing the 
corresponding verified mobile station location. Accordingly, a table of partition area probabilities can be 
determined for each first order model 1224. Thus, when a location hypothesis is generated and identified 

1 5 as belonging to one of the partition areas, the corresponding probability for that partition area may be 

assigned as the confidence value for the location hypothesis. The advantages to using actual probabilities 
here is that, as will be discussed below, the most likelihood estimator 1344 can compute a straightforward 
probability for each distinct intersection of the multiple location hypotheses generated by the multiple first 
order models, such that each such probability indicates a likelihood that the target mobile station is in the 

20 corresponding intersection. 

Location Information Data Bases And Data 

Location Data Bases Introduction 

It is an aspect of wireless location embodiments of the novel method and/or system disclosed 
herein that MS location processing performed by the location center 142 should become increasingly better at 

25 locating a target MS 1 40 both by (a) building an increasingly more detailed model of the signal characteristics 
of locations in the service area for such a wireless location embodiment , and also (b) by providing capabilities 
for the location center processing to adapt to environmental changes. 

One way the functionality of (a) and (b) immediately above is realized is by providing one or more 
data base management systems and data bases for 

30 (a) storing and associating wireless MS signal characteristics with known locations of MSs 140 

used in providing the signal characteristics. Such stored associations may not only provide an increasingly 
better model of the signal characteristics of the geography of the service area, but also provide an 
increasingly better model of more changeable signal characteristic affecting environmental factors such as 
weather, seasons, and/or traffic patterns; 

35 (b) adaptively updating the signal characteristic data stored so that it reflects changes in the 

environment of the service area such as, for example, a new high rise building or a new highway. 
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Referring again to Fig. 5 of the collective representation of these data bases is the location 
information data bases 1232. Included among these data bases is a data base for providing training and/or 
calibration data to one or more trainable/calibratable FOMs 1224, as well as an archival data base for 
archiving historical MS location information related to the performance of the FOMs. These data bases will be 
5 discussed as necessary hereinbelow. However, a further brief introduction to the archival data base is 

provided here. Accordingly, the term, "location signature data base" is used hereinafter to denote the archival 
data base and/or data base management system depending on the context of the discussion. The location 
signature data base (shown in, for example, Fig. 6 and labeled 1 320) is a repository for wireless signal 
characteristic data derived from wireless signal communications between an MS 140 and one or more base 

l o stations 1 22, wherein the corresponding location of the MS 140 is known and also stored in the location 
signature data base 1320. More particularly, the location signature data base 1320 associates each such 
known MS location with the wireless signal characteristic data derived from wireless signal communications 
between the MS 140 and one or more base stations 1 22 at this MS location. Accordingly, it is an aspect of 
wireless location embodiments of the novel method and/or system disclosed herein to utilize such historical 

1 5 MS signal location data for enhancing the correctness and/or confidence of certain location hypotheses as will 
be described in detail in other sections below. 

Data Representations for the Location Signature Data Base 

There are four fundamental entity types (or object classes in an object oriented programming 
paradigm) utilized in the location signature data base 1 320. Briefly, these data entities are described in the 
20 items (24.1 ) through (24.4) that follow: 

(24.1 ) (verified) location signatures: Each such (verified) location signature describes the wireless signal 
characteristic measurements between a given base station (e.g., BS 122 or LBS 152) and an MS 140 at a 
(verified or known) location associated with the (verified) location signature. That is, a verified location 
signature corresponds to a location whose coordinates such as latitude-longitude coordinates are known, 

25 while simply a location signature may have a known or unknown location corresponding with it. Note that the 
term (verified) location signature is also denoted by the abbreviation, "(verified) loc sig" hereinbelow; 

(24.2) (verified) location signature clusters: Each such (verified) location signature cluster includes a 
collection of (verified) location signatures corresponding to all the location signatures between a target MS 
140 at a (possibly verified) presumed substantially stationary location and each BS (e.g., 122 or 152) from 

30 which the target MS 140 can detect the BS's pilot channel regardless of the classification of the BS in the 
target MS (i.e., for CDMA, regardless of whether a BS is in the MS's active, candidate or remaining base 
station sets, as one skilled in the art will understand). Note that for simplicity here, it is presumed that each 
location signature cluster has a single fixed primary base station to which the target MS 140 synchronizes or 
obtains its timing; 

35 (24.3) "composite location objects (or entities)": Each such entity is a more general entity than the verified 
location signature cluster. An object of this type is a collection of (verified) location signatures that are 
associated with the same MS 140 at substantially the same location at the same time and each such loc sig is 
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associated with a different base station. However, there is no requirement that a loc sig from each BS 1 22 
for which the MS 140 can detect the BS's pilot channel is included in the "composite location object (or 
entity)"; and 

(24.4) MS location estimation data that includes MS location estimates output by one or more MS location 
. 5 estimating first order models 1224, such MS location estimate data is described in detail hereinbelow. 

It is important to note that a loc sig is, in one embodiment, an instance of the data structure 
containing the signal characteristic measurements output by the signal filtering and normalizing 
subsystem also denoted as the signal processing subsystem 1220 describing the signals between: (i) a 
specific base station 122 (BS) and (ii) a mobile station 140 (MS), wherein the BS's location is known and 
10 the MS's location is assumed to be substantially constant (during a 2 to 5 second interval in one 

embodiment of the novel method and/or system), during communication with the MS 140 for obtaining a 
single instance of loc sig data, although the MS location may or may not be known. Further, for notational 
purposes, the BS 122 and the MS 140 for a loc sig hereinafter will be denoted the "BS associated with the 
loc sig", and the "MS associated with the loc sig" respectively. Moreover, the location of the MS 140 at 
1 5 the time the loc sig data is obtained will be denoted the "location associated with the loc sig" (this location 
possibly being unknown). 

Loc sigs have the following functions or object methods associated therewith: 
(26.1) A "normalization" method for normalizing loc sig data according to the associated MS 140 and/or BS 

20 1 22 signal processing and generating characteristics. That is, the signal processing subsystem 1 220, 

one embodiment being described in the PCT patent application PCT/US97/15933, titled, "Wireless 
Location Using A Plurality of Commercial Network Infrastructures," now U.S. Patent No. 6,236,365, by 
F. W. LeBlanc and the present inventor, provides (methods for loc sig objects) for "normalizing" each 
loc sig so that variations in signal characteristics resulting from variations in (for example) MS signal 

25 processing and generating characteristics of different types of MS's may be reduced. In particular, 

since wireless network designers are typically designing networks for effective use of hand set MS's 
140 having a substantially common minimum set of performance characteristics, the normalization 
methods provided here transform the loc sig data so that it appears as though the loc sig was 
provided by a common hand set MS 140. However, other methods may also be provided to 

30 "normalize" a loc sig so that it may be compared with loc sigs obtained from other types of MS's as 

well. Note that such normalization techniques include, for example, interpolating and extrapolating 
according to power levels so that loc sigs may be normalized to the same power level for, e.g., 
comparison purposes. 

Normalization for the BS 1 22 associated with a loc sig is similar to the normalization for MS signal 
35 processing and generating characteristics. Just as with the MS normalization, the signal processing 

subsystem 1220 provides a loc sig method for "normalizing" loc sigs according to base station signal 
processing and generating characteristics. 



40 



Note, however, loc sigs stored in the location signature data base 1320 are NOT 
"normalized" according to either MS or BS signal processing and generating characteristics. That is, 
"raw" values of the wireless signal characteristics are stored with each loc sig in the location signature 
data base 1320. 

5 (26.2) A method for determining the "area type" corresponding to the signal transmission characteristics of the 
area(s) between the associated BS 122 and the associated MS 140 location for the loc sig. Note, 
such an area type may be designated by, for example, the techniques for determining transmission 
area types as described hereinabove. 
(26.3) Other methods are contemplated for determining additional environmental characteristics of the 
l o geographical area between the associated BS 1 22 and the associated MS 1 40 location for the loc sig; 

e.g., a noise value indicating the amount of noise likely in such an area. 

Referring now to the composite location objects and verified location signature clusters of (24.3) and 
(24.2) respectively, the following information is contained in these aggregation objects: 
15 (27.1 .1) an identification of the BS 122 designated as the primary base station for communicating with the 
target MS 140; 

(27.1 .2) a reference to each loc sig in the location signature data base 1320 that is for the same MS location 
at substantially the same time with the primary BS as identified in (27.1); 

(27.1 .3) an identification of each base station (e.g., 122 and 152) that can be detected by the MS 140 at the 
20 time the location signal measurements are obtained. Note that in one embodiment, each composite 

location object includes a bit string having a corresponding bit for each base station, wherein a "1 " for 
such a bit indicates that the corresponding base station was identified by the MS, and a "0" indicates 
that the base station was not identified. In an alternative embodiment, additional location signal 
measurements may also be included from other non-primary base stations. For example, the target 

25 MS 140 may communicate with other base stations than ifs primary base station. However, since 

the timing for the MS 140 is typically derived from ifs primary base station and since timing 
synchronization between base stations is not exact (e.g., in the case of CDMA, timing variations may 
be plus or minus 1 microsecond)at least some of the location signal measurements may be less 
reliable that the measurements from the primary base station, unless a forced hand-off technique is 

30 used to eliminate system timing errors among relevant base stations; 

(27.1 .4) a completeness designation that indicates whether any loc sigs for the composite location object 
have been removed from (or invalidated in) the location signature data base 1320. 

Note, a verified composite location object is designated as Incomplete" if a loc sig initially referenced 
35 by the verified composite location object is deleted from the location signature data base 1 320 (e.g., because 
of a confidence that is too tow). Further note that if all loc sigs for a composite location object are deleted, then 
the composite object is also deleted from the location signature data base 1320. Also note that common 
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fields between loc sigs referenced by the same composite location object may be provided in the composite 
location object only (e.g., timestamp, etc.). 

Accordingly, a composite location object that is complete (i.e., not incomplete) is a verified location 
signature cluster as described in (24.2). 

5 Location Center Architecture 

Overview of Location Center Functional Components 

Fig. 5 presents a high level diagram of the location center 142 and the location engine 1 39 in the 
context of the infrastructure for the entire location system according to the novel method and/or system 
disclosed herein. 

10 It is important to note that the architecture for the location center 142 and the location engine 1 39 

disclosed herein is designed for extensibility and flexibility so that MS 140 location accuracy and reliability may 
be enhanced as further location data become available and as enhanced MS location techniques become 
available. In addressing the design goals of extensibility and flexibility, the high level architecture for 
generating and processing MS location estimates may be considered as divided into the following high level. 

1 5 functional groups described hereinbelow. 

Low Level Wireless Signal Processing Subsystem for Receiving and Conditioning Wireless Signal 
Measurements 

A first functional group of location engine 139 modules is for performing signal processing and 
filtering of MS location signal data received from a conventional wireless (e.g., CDMA) infrastructure, as 
20 discussed in the steps (23.1) and (23.2) above. This group is denoted the signal processing subsystem 1 220 
herein. One embodiment of such a subsystem is described in the PCT patent application PCT/US97/1 5933, 
titled, "Wireless Location Using A Plurality of Commercial Network Infrastructures," now U.S. Patent No. 
6,236,365 by F. W. LeBlanc, Karr, and the present inventor filed July 8, 1999. 

Initial Location Estimators: First Order Models 

25 A second functional group of location engine 1 39 modules is for generating various target MS 140 

location initial estimates, as described in step (23.3). Accordingly, the modules here use input provided by the 
signal processing subsystem 1220. This second functional group includes one or more signal analysis 
modules or models, each hereinafter denoted as a first order model 1 224 (FOM), for generating location 
hypotheses for a target MS 140 to be located. Note that it is intended that each such FOM 1 224 use a 

30 different technique for determining a location area estimate for the target MS 1 40. A brief description of some 
types of first order models is provided immediately below. Note that Fig. 8 illustrates another, more detailed 
view of the location system according to the novel method and/or system disclosed herein. In particular, this 
figure illustrates some of the FOMs 1224 contemplated by a wireless location embodiment of the novel 
method and/or system, and additionally illustrates the primary communications with other modules of the 

35 present wireless location system. However, it is important to note that the embodiments are not limited to the 
FOMs 1 224 shown and discussed herein. That is, it is a primary aspect of wireless location embodiments of 



the novel method and/or system to easily incorporate FOMs using other signal processing and/or 

computational location estimating techniques than those presented herein. Further, note that each FOM type 

may have a plurality of its models incorporated into an embodiment of a wireless location system according to 

the novel method and/or system disclosed herein. 
5 For example, (as will be described in further detail below), one such type of model or FOM 1 224 

(hereinafter models of this type are referred to as "distance models") may be based on a range or distance 

computation and/or on a base station signal reception angle determination between the target MS 140 from 

each of one or more base stations. Basically, such distance models 1224 determine a location estimate of 

the target MS 1 40 by determining a distance offset from each of one or more base stations 1 22, possibly in a 

l o particular direction from each (some of) the base stations, so that an intersection of each area locus defined 

by the base station offsets may provide an estimate of the location of the target MS 

In one embodiment, such a distance model may perform the following steps: 

(a) Determines a minimum distance between the target MS and each BS using TOA, TDOA, 

signal strength on both forward and reverse paths; 

15 (b) Generates an estimated error; 

(c) Outputs a location hypothesis for estimating a location of a MS: each such hypothesis 

having: (a) one or more (nested) location area estimates for the MS, each location estimate 

having a confidence value (e.g., provided using the estimated error) indicating a perceived 

accuracy, and (b) a reason for both the location estimate (e.g., substantial multipath, etc) 

20 and the confidence. 

Another type of FOM 1 224 is a statistically based first order model 1 224, wherein a statistical 

technique, such as regression techniques (e.g., least squares, partial least squares, principle decomposition), 

or e.g., Bollenger Bands (e.g., for computing minimum and maximum base station offsets). In general, 

models of this type output location hypotheses that are determined by performing one or more statistical 

25 techniques or comparisons between the verified location signatures in location signature data base 1 320, and 

the wireless signal measurements from a target MS. Models of this type are also referred to hereinafter as a 

"stochastic signal (first order) moder or a "stochastic FOM" or a "statistical model." 

In one embodiment, such a stochastic signal model may output location hypotheses determined by 

one or more statistical comparisons with loc sigs in the Location Signature database 1320 (e.g., 

30 comparing MS location signals with verified signal characteristics for predetermined geographical areas). 

Still another type of FOM 1 224 is an adaptive learning model, such as an artificial neural net or a 

genetic algorithm, wherein the FOM may be trained to recognize or associate each of a plurality of locations 

with a corresponding set of signal characteristics for communications between the target MS 140 (at the 

location) and the base stations 1 22. Moreover, typically such a FOM is expected to accurately 

3 5 interpolate/extrapolate target MS 1 40 location estimates from a set of signal characteristics from an unknown 

target MS 140 location. Models of this type are also referred to hereinafter variously as "artificial neural net 

models" or "neural net models" or "trainable models" or "learning models." Note that a related type of FOM 



1224 is based on pattern recognition. These FOMs can recognize patterns in the signal characteristics of 
communications between the target MS 140 (at the location) and the base stations 122 and thereby estimate 
a location area of the target MS. However, such FOMs may not be trainable. 

In one embodiment, an adaptive learning model such as a model based on an artificial neural 
5 network may determine an MS 140 location estimate using base station IDs, data on signal-to-noise, other 
signal data (e.g., a number of signal characteristics including, e.g., all CDMA fingers). Moreover, the 
output from such a model may include: a latitude and longitude for a center of a circle having radius R (R 
may be an input to such an artificial neural network), and is in the output format of the distance model(s). 
Yet another type of FOM 1224 can be based on global positioning satellite (GPS) signals. For 
10 example, in some embodiments of the MS 140, there may be electronics that at least are capable of receiving 
GPS signals and using such signals for determining the position of the MS 1 40. In general, if an appropriate 
number of GPS signals from different satellites are received by the MS 140 (e.g., signals from two, three or 
four satellites depending on the GPS technology used and the location resolution desired), such signals may 
provide highly accurate location estimates. However, such signals may not be available in, e.g., high rise 
15 buildings, parking garages, and dense urban canyons. Moreover, such signals can in some circumstances be 
subject to the multipath phenomenon which could cause gross errors in a location estimate of the target MS 
140, or no location estimate being generated. 

Yet another type of FOM 1 224 can be based on a collection of dispersed low power, low cost fixed 
location wireless transceivers (also denoted "location base stations 152" hereinabove) that are provided for 
20 detecting a target MS 140 in areas where, e.g., there is insufficient base station 1 22 infrastructure coverage 
for providing a desired level of MS 140 location accuracy. 

In one embodiment, such a location base station model may perform the following steps: 
(a) If an input is received then the target MS 140 is detected by a location base station 152 (i.e., a 
LBS being a unit having a reduced power BS and a MS). 
25 (b) If an input is obtained, then the output is a hypothesis data structure having a small area of the 

highest confidence. 

(c) If no input is received from a LBS then a hypothesis having an area with highest negative 
confidence is output. 

Yet another type of FOM 1224 can be based on input from a mobile base station 148, wherein 
30 location hypotheses may be generated from target MS 140 location data received from the mobile base 
station 1 48. In one embodiment, such a mobile base station model may provide output similar to the 
distance FOM 1224 described hereinabove. 

Still other types of FOM 1224 can be based on various techniques for recognizing wireless signal 
measurement patterns and associating particular patterns with locations in the coverage area 120. For 
35 example, artificial neural networks or other learning models can used as the basis for various FOMs. 
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It is important to keep in mind that a novel aspect of the novel method and/or system disclosed 
herein is the use or activation, simultaneously or serially, of a potentially large number of such first order 
models 1224, wherein such FOMs are not limited to those described herein. Thus, a wireless location system 
according to the novel method and/or system provides a framework for incorporating MS location estimators 
5 to be subsequently provided as new FOMs in a straightforward manner. For example, a FOM 1 224 based on 
wireless signal time delay measurements from a distributed antenna system 168 for wireless communication 
may be incorporated into a wireless location system according to the novel method and/or system for 
locating a target MS 140 in an enclosed area serviced by the distributed antenna system (such a FOM is 
more fully described in the U.S. Patent 6,236,365 filed July 8, 1999). Accordingly, by using such a distributed 
l o antenna FOM, such a wireless location system may determine the floor of a multi-story building from which a 
target MS is transmitting. Thus, MSs 140 can be located in three dimensions using such a distributed 
antenna FOM. 

It is important to note the following aspects of the novel method and/or system relating to FOMs 

1224: 

1 5 (28.1 ) Each such first order model 1 224 may be relatively easily incorporated into and/or removed from an 
embodiment of the novel method and/or system disclosed herein. For example, assuming that the 
signal processing subsystem 1220 provides uniform input to the FOMs, and there is a uniform FOM 
output interface, it is believed that a large majority (if not substantially all) viable MS location estimation 
strategies may be accommodated. Thus, it is straightforward to add or delete such FOMs 1224. 

20 (28.2) Each such first order model 1224 may be relatively simple and still provide significant MS 140 locating 
functionality and predictability. For example, much of what is believed to be common or generic MS 
location processing has been coalesced into, for example: a location hypothesis evaluation subsystem, 
denoted the hypotheses evaluator 1228 and described immediately below. Thus, a wireless location 
system according to the novel method and/or system disclosed herein may be modular and 

25 extensible such that, for example, (and importantly) different first order models 1 224 may be utilized 

depending on the signal transmission characteristics of the geographic region serviced by the wireless 
location system. Thus, a simple configuration of the wireless location system may have a small 
number of FOMs 1224 for a simple wireless signal environment (e.g., flat terrain, no urban canyons 
and low population density). Alternatively, for complex wireless signal environments such as in cities 

30 like San Francisco, Tokyo or New York, a large number of FOMs 1 224 may be simultaneously or 

serially utilized for generating MS location hypotheses. 
(28.3) A wireless location system according to the novel method and/or system disclosed herein provides 
the capability to activate serially one or more groups of one or more models 1224, wherein each 
model may generate a location estimate for locating the same target MS 140. Thus, for example, if 

35 a first group of models do not yield an appropriate location estimate (e.g., having a sufficiently high 

confidence value), then another group of models may be activated. Accordingly, a first group may 
include a model based on global positioning satellite (GPS) technology, whereas a second group 
may include base station triangulation models and pattern recognition models. Thus, for a target 
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MS 140 having the capability to receive GPS signals, if measurements of such signals can be 
received from the MS, then in many situations a location estimate derived therefrom will be 
sufficient for locating the MS for such applications as embodiments of the novel method and/or 
system is intended (e.g., E-91 1 , auto tracking, selective wireless advertising based on location, 
5 railroad box car tracking, tracking of trucks between cities or states, penal restrictions on the 

location areas to which an individual has access). However, if such GPS signals are not available, 
due to environmental conditions and/or due to the target MS 140 not having the electronics for 
receiving GPS signals, then a second group of models may be activated. This capability for 
serially activating groups of models may be embodied using a data driven approach wherein the 

1 o location center control subsystem 1 350 activates, e.g., all models that are capable of being 

activated from the types of wireless measurement data received from a target MS 140 according to 
a particular priority. Thus, assuming GPS data is received, possibly along with additional signal 
time delay measurements, then a GPS location estimating model may be invoked as a first level of 
location processing. Thus, if each confidence value for the (one or more) GPS model generated 

1 5 location hypotheses have a sufficiently high probability (and substantially overlap if more than one), 

then no further processing may be required, thereby conserving computational resources of the 
location center 142. However, in the case that such confidence values are not sufficient (or that 
the estimates do not substantially overlap), then an additional second group of models 1224 may 
be activated that are based on other computational techniques for obtaining estimates of the target 

20 MS 140. Thus, the second group may include the models 1 224 that perform various forms of 

pattern recognition. Additionally, if the combination of the first and second groups do not provide 
an estimate of sufficiently high confidence, then an third group of, e.g., triangulation based models 
1224 may be invoked. 

An Introduction to an Evaluator for Location Hypotheses: Hypothesis Evaluator 

25 A third functional group of location engine 1 39 modules evaluates location hypotheses output by the 

first order models 1 224 and thereby provides a "most likely" target MS location estimate. The modules for this 
functional group are collectively denoted the hypothesis evaluator 1228. 

Hypothesis Evaluator Introduction 

A primary purpose of the hypothesis evaluator 1228 is to mitigate conflicts and ambiguities related to 
30 location hypotheses output by the first order models 1 224 and thereby output a "most likely" estimate of an 
MS for which there is a request for it to be located. In providing this capability, there are various related 
embodiments of the hypothesis evaluator that are within the scope of the novel method and/or system 
disclosed herein. Since each location hypothesis includes both an MS location area estimate and a 
corresponding confidence value indicating a perceived confidence or likelihood of the target MS being within 
35 the corresponding location area estimate, there is a monotonic relationship between MS location area 
estimates and confidence values. That is, by increasing an MS location area estimate, the corresponding 
confidence value may also be increased (in an extreme case, the location area estimate could be the entire 
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coverage area 120 and thus the confidence value may likely correspond to the highest level of certainty; i.e., 
+1 .0). Accordingly, given a target MS location area estimate (of a location hypothesis), an adjustment to its 
accuracy may be performed by adjusting the MS location area estimate and/or the corresponding confidence 
value. Thus, if the confidence value is, for example, excessively low then the area estimate may be increased 
5 as a technique for increasing the confidence value. Alternatively, if the estimated area is excessively large, 
and there is flexibility in the corresponding confidence value, then the estimated area may be decreased and 
the confidence value also decreased. Thus, if at some point in the processing of a location hypothesis, if the 
location hypothesis is judged to be more (less) accurate than initially determined, then (i) the confidence value 
of the location hypothesis may be increased (decreased), and/or (ii) the MS location area estimate can be 

10 decreased (increased). Moreover, note that when the confidence values are probabilities, such adjustments 
are may require the reactivation of one or more FOMs 1224 with requests to generate location hypotheses 
having location estimates of different sizes. Alternatively, adjuster modules 1436 and/or 1440 (Fig. 15 
discussed hereinbelow) may be invoked for generating location hypotheses having area estimates of different 
sizes. Moreover, the confidence value on such an adjusted location hypothesis (actually a new location 

1 5 hypothesis corresponding to the originally generated hypothesis) may also be a probability in that 

combinations of FOMs 1 224 and adjuster modules 1 436 and 1440 can also be calibrated for thereby yielding 
probabilities as confidence values to the resulting location hypotheses. 

In a first class of embodiments (typically wherein the confidence values are not maintained as 
probabilities), the hypothesis evaluator 1228 evaluates location hypotheses and adjusts or modifies only their 

20 confidence values for MS location area estimates and subsequently uses these MS location estimates with 
the adjusted confidence values for determining a "most likely" MS location estimate for outputting. 
Alternatively, in a second class of embodiments for the hypothesis evaluator 1 228 (also typically wherein the 
confidence values are not maintained as probabilities), MS location area estimates can be adjusted while 
confidence values remain substantially fixed. However, in one preferred embodiment of the present 

25 embodiment, both location hypothesis area estimates and confidence values are modified. 

The hypothesis evaluator 1228 may perform any or most of the following tasks depending on the 
embodiment of the hypothesis evaluator. That is, 

(30.1) it may enhance the accuracy of an initial location hypothesis generated by an FOM by using the 

initial location hypothesis as, essentially, a query or index into the location signature data base 
30 1320 for obtaining one or more corresponding enhanced location hypotheses, wherein the 

enhanced location hypotheses have both an adjusted target MS location area estimates and an 
adjusted confidences based on past performance of the FOM in the location service surrounding 
the target MS location estimate of the initial location hypothesis; 
Additionally, for embodiments of the hypothesis evaluator 1228 wherein the confidence values for location 
35 hypotheses are not maintained as probabilities, the following additional tasks (30.2) through (30.7) may be 
performed: 

(30.2) the hypothesis evaluator 1228 may utilize environmental information to improve and reconcile location 

hypotheses supplied by the first order models 1224. A basic premise in this context is that the 
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accuracy of the individual first order models may be affected by various environmental factors such 
as, for example, the season of the year, the time of day, the weather conditions, the presence of 
buildings, base station failures, etc.; 

(30.3) the hypothesis evaluator 1 228 may determine how well the associated signal characteristics used for 
locating a target MS compare with particular verified loc sigs stored in the location signature data base 
1320 (see the location signature data base section for further discussion regarding this aspect of the 
novel method and/or system). That is, for a given location hypothesis, verified loc sigs (which were 
previously obtained from one or more verified locations of one or more MS's) are retrieved for an area 
corresponding to the location area estimate of the location hypothesis, and the signal characteristics of 
these verified loc sigs are compared with the signal characteristics used to generate the location 
hypothesis for determining their similarities and subsequently an adjustment to the confidence of the 
location hypothesis (and/or the size of the location area estimate); 

(30.4) the hypothesis evaluator 1228 may determine if (or how well) such location hypotheses are consistent 
with well known physical constraints such as the laws of physics. For example, if the difference 
between a previous (most likely) location estimate of a target MS and a location estimate by a current 
location hypothesis requires the MS to: 

(a1 ) move at an unreasonably high rate of speed (e.g., 200 mph), or 
(b1 ) move at an unreasonably high rate of speed for an area (e.g., 80 mph in a com patch), 
or 

(d) make unreasonably sharp velocity changes (e.g., from 60 mph in one direction to 60 
mph in the opposite direction in 4 sec), then the confidence in the current Location 
Hypothesis is likely to be reduced. 
Alternatively, if for example, the difference between a previous location estimate of a target MS 
and a current location hypothesis indicates that the MS is: 

(a2) moving at an appropriate velocity for the area being traversed, or 
(b2) moving along an established path (e.g., a freeway), 
then the confidence in the current location hypothesis may be increased. 

(30.5) the hypothesis evaluator 1228 may determine consistencies and inconsistencies between location 
hypotheses obtained from different first order models. For example, if two such location hypotheses, 
for substantially the same timestamp, have estimated location areas where the target MS is likely to be 
and these areas substantially overlap, then the confidence in both such location hypotheses may be 
increased. Additionally, note that a velocity of an MS may be determined (via deltas of successive 
location hypotheses from one or more first order models) even when there is low confidence in the 
location estimates for the MS, since such deltas may, in some cases, be more reliable than the actual 
target MS location estimates; 

(30.6) the hypothesis evaluator 1 228 determines new (more accurate) location hypotheses from other 
location hypotheses. For example, this module may generate new hypotheses from currently active 
ones by decomposing a location hypothesis having a target MS location estimate intersecting two 
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radically different area types. Additionally, this module may generate location hypotheses 
indicating areas of poor reception; and 
(307) the hypothesis evaluator 1 228 determines and outputs a most likely location hypothesis for a target 
MS. 

5 

Note that the hypothesis evaluator may accomplish the above tasks, (30.1) - (30.7), by employing 
various data processing tools including, but not limited to, fuzzy mathematics, genetic algorithms, neural 
networks, expert systems and/or blackboard systems. 

Note that in the embodiments where confidence values are not maintained as probabilities, such as 

l o those of Figs. 6 and 7, the hypothesis evaluator 1 228 includes the following four high level modules for 

processing output location hypotheses from the first order models 1224: a context adjuster 1326, a hypothesis 
analyzer 1 332, an MS status repository 1 338 and a most likelihood estimator 1 334. These four modules are 
briefly described hereinbelow. Alternatively, Fig. 1 5 shows an embodiment of the hypothesis evaluator 1 228 
where confidence values are maintained as probabilities. 

1 5 As mentioned above, each of the first order models 1 224 have default confidence values 

associated therewith, and these confidence values may be probabilities. More precisely, such probability 
confidence values can be determined as follows. Assume there is a partition of the coverage area into 
subareas, each subarea being denoted a "partition area." For each partition area, activate each first order 
model 1224 with historical location data in the Location Signature Data Base 1320 (Fig. 6), wherein the 

20 historical location data has been obtained from corresponding known mobile station locations in the 
partition area. For each first order model, determine a probability of the first order model generating a 
location hypothesis whose location estimate contains the corresponding known mobile station location. 
To accomplish this, assume the coverage area is partitioned into partition areas A, wherein each partition 
area A is specified as the collection of coverage area locations such that for each location, the detected 

25 wireless transmissions between the network base stations and a target mobile station at the location can 
be straightforwardly equated with other locations of area A. For example, one such partition, Po, can be 
defined wherein each partition area A is specified in terms of three sets of base station identifiers, namely, 
(a) the base station identifiers of the base stations that can be both detected at each location of A and can 
detect a target mobile station at each location, (b) the identifiers for base stations that can detect a target 

30 mobile station at each location of A, but can not be detected by the target mobile station, and (c) the 

identifiers for base stations that can be detected by a target mobile station at each location of A, but these 
base stations can not detect the target mobile station. That is, two locations, U and b. are identified as 
being in A if and only if the three sets of (a), (b), and (c) for h are, respectively, identical to the three sets 
of (a), (b), and (c) for l 2 . Other embodiments of Po may be provided by the area types described in the 

35 "Coverage Area: Area Types And Their Determination" section. 

Accordingly, assuming such a partition Po is used, a description can be given as to how 
probabilities may be assigned as the confidence values of location hypotheses generated by the first 
order models 1224. For each partition area A, a first order model 1224 is supplied with wireless 
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measurements of archived location data in the Location Signature Data Base associated with 
corresponding verified mobile station locations. Thus, a probability can be determined as to how likely the 
first order model is to generate a location hypothesis having a location estimate containing the 
corresponding verified mobile station location. Accordingly, a table of partition area probabilities can be 
5 determined for each first order model 1224. Thus, when a location hypothesis is generated and identified 
as belonging to one of the partition areas, the corresponding probability for that partition area may be 
assigned as the confidence value for the location hypothesis. The advantages to using actual probabilities 
here is that, as will be discussed below, the most likelihood estimator 1344 can compute a straightforward 
probability for each distinct intersection of the multiple location hypotheses generated by the multiple first 
10 order models, such that each such probability indicates a likelihood that the target mobile station is in the 
corresponding intersection. 

Reliability and Accuracy of a Multiple Model Approach 

To illustrate how both reliability and accuracy can both be enhanced by this multiple model 
architecture of the novel method and/or system disclosed herein, examples believed to illustrate the 
1 5 plausibility of such enhancements immediately follow. However, note that the general architecture is 
within the scope of the novel method and/or system independently of these examples. 

(i) Theoretical Computations Related To Location Estimation Reliability 

Regarding reliability, let the reliability for each first order model in a particular area, A, be defined 
as the probability of a target MS being in a corresponding location estimate of a location hypothesis 

20 generated by the first order model, wherein the probability is determined from an analysis of past 

performance of the model. Thus, a first order model whose past location estimates were shown to have 
included the corresponding target MS 6 out of 10 times has a reliability of 0.6. Now suppose there are 
five first order models, Fi, F 2 , F 3 , F 4 , and F 5 . Further, suppose that historically these models have the 
following reliabilities: Fi is 0.7, F 2 is 0.6, F 3 is 0.55, F 4 is 0.8, and F 5 is 0.75. Additionally, assume that 

25 these models Fi are substantially "independent" in a statistical sense. For a particular target MS, also 
suppose that each of the models Fj generates a single corresponding location hypothesis, Hi, 1 <= i <= 5 
for estimating a location of the particular target MS. Accordingly, one "system reliability" measurement of 
the reliability of the combination of all the models is the probability that the particular target MS resides in 
a subcombination of the location hypotheses Hi, H2, H3, H4, and H5 simultaneously, wherein, for example, 

30 the subcombination always includes at least three of the location hypothesis Hi. Accordingly, this 
probability can be determined using the reliabilities of the first order models Fj. Accordingly, this 
probability is approximately 0.92. Note that this resulting reliability is substantially greater than the 
reliability of each of the first order models. 

Moreover, suppose a sixth and seventh first order model, F6 and F 7 , are introduced, and further 

35 suppose that each of these models have only 0.5 reliability in area A (i.e., a target MS is as likely to be in 

a location estimate generated by this model as it is not). Additionally, suppose that the criteria for 

measuring reliability is made more rigorous in that it is required that the target MS must be in four location 
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estimates simultaneously. Under these conditions, the reliability is approximately 0.85. That is, even with 
a more rigorous reliability definition, and with a majority of the first order models providing low reliability 
(i.e., less than or equal to 0.60), the synergies of the simultaneous location estimates produces a 
substantial magnification in the reliability of the entire location system. 
5 It is important, however, to keep in mind that, as mentioned above, the analysis performed here 

assumes that the models are substantially independent. It is believed that seven such independent first 
order models can be obtained that substantially satisfy this constraint.. Also note that in many cases 
multiple instantiations of such models can be activated simultaneously with wireless signal measurements 
from different collections of base stations. Moreover, since the forward transmissions from the base 

1 o station network to a target MS is in a different bandwidth from the reverse transmissions from the target 
MS to the base station network, the forward and reverse signals can behave differently. Thus, for a first 
order model based on forward signal measurements, there is likely to be a corresponding first order 
model based on reverse signal measurements that has a similar computational technique. It is also 
important to note that the analysis performed here does not take into account the location hypothesis 

15 adjustment step mentioned above (and described further below). Accordingly, it is believed that even 
better system reliabilities can be obtained by performing the adjusting step. 

(ii) Theoretical Computations Related To Location Estimation Accuracy 

Turning now to the accuracy aspect of wireless location, accuracy as used herein is intended to 
substantially correspond to the area sizes of MS location estimates, that is, assuming the shapes of the 

20 location estimates are not, in general, grossly oblong. Accordingly, circles, squares, hexagons and other 
- similar polygonal area shapes are considered appropriate. Moreover, for simplicity it is assumed that 
there are no "holes" (i.e., excluded areas encircled by a location estimate) in the location estimates. 
However, even with such assumptions in place, to compute accuracy by computing expected area sizes 
for overlaps in location estimates is beyond the scope of this document. Accordingly, in order to provide a 

25 quantitative indication as to accuracy enhancements possible by using multiple simultaneous independent 
MS location estimates, squares of the same size will be assumed for the location estimates generated by 
the first order models. Thus, to estimate accuracy, an average or expected area size overlap between 
two of the generated square location estimates will be first computed. 

To compute the average area size overlap, a grid is superimposed on each such generated 

30 location estimate square so that a cell mesh of smaller cell squares is obtained. In particular, the mesh 
provided is an N by N array of cells, where N >=3. Additionally, to estimate an average size of an overlap 
(i.e., intersection) between two generated location estimate squares, having N by N cell meshes, also 
assume that each intersection must be on mesh boundaries of the intersecting generated location 
estimates. The following diagram illustrates two generated location estimate squares overlapping in the 

35 two pattern filled smaller cell mesh squares (for a 4 by 4 cell mesh). 
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Further, assume that all possible overlaps have the same probability of occurring. 

Accordingly, the area size of an average intersection between two such generated location 
estimate squares can be computed by determining the number of different intersections and the areas of 

5 each intersection. Note that for two location estimate squares having N by N cell meshes, the number of 
different intersections is: [N * (N - 1) * 4] +1. Thus, using the average intersection size, a ratio of the area 
of the average intersection to the total area of one of the generated location estimate squares can be 
computed. Unfortunately, this ratio varies with the number N. However, the ratio appears to 
asymptotically converge as TABLE A below illustrates (i.e., the deltas of the third column). Thus, it is 

10 believed that this ratio is generally less than 0.3 of the original square area as the second column of the 
following table indicates. 
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Thus, using 0.3 as a conservative estimate of the ratio of the area of the average intersection to 
the area of one of the generated location estimate squares, this implies that an average intersection size 
will be only 0.3 of the area of one of the generated location estimates. Additionally, note that in all 
5 examples of TABLE A, by comparing the average number of mesh cells in an intersection (column four) 
with the corresponding average difference between the length and the width of an intersection (column 
five), it can be seen that a typical intersection is very close to also being a square. 

Accordingly, to illustrate the accuracy enhancements possible simply by using multiple 
simultaneous independent MS location estimates, let h.2 be a non-empty intersection of generated 

10 location estimate squares Li and L2, and, let h,4 be a non-empty intersection of generated location 

estimate squares U and U Then as TABLE A illustrates, both intersections are likely to be very close to 
being squares. However, h i2 and b,4 may be of very different sizes. Assume that h,2 is of a size less than 
or equal to that of b.4, then if cell rows and/or cell columns were added to li ( 2 so that it has a size 
comparable to h.4 and additionally is square (or substantially so), then it is expected that the intersection 

15 of h,2 and l 3 ,4 would be 0.3 (or less) of 13,4. Thus, if the generated location estimate squares Li, L 2 , L3, and 
U are generated in a manner that is substantially independent of the way the others are generated, then it 
is expected that the reduction in area size over the size of the original generated location estimate 
squares is on the order of 0.09. 

Referring back to the seven first order models, Fj, discussed above, an example can now be 

20 given tying the above accuracy analysis with the reliability analysis preceding it. If each of the seven first 
order models R simultaneously generates a location hypothesis having a corresponding location estimate 
square of equal size, then it is expected that approximately 85% of the time a simultaneous intersection of 
four of the location estimates contains the target MS, and further that the intersection of these location 
estimates yields an area that is likely nearly square with an area on the order of 0.09 of one of the original 

25 location estimate squares. Thus, for illustration, if each of the original generated location estimate 

squares is approximately 1,460 feet on a side, then such squares can be circumscribed by a circle with a 
diameter of approximately 2,065 feet (i.e., about 0.4 miles). Further, the simultaneous intersection of four 
of these squares is likely to provide an area of 191,844 sq. ft. in a generally square area. Thus, this 
intersection area should likely be close to a square of 438 feet per side. Accordingly, such a square can 

30 be contained within circular area, C, having a diameter of 61 9 feet. Accordingly, approximately 85% of 
the time the target MS is within 31 0 feet of the center of C. 

Thus, in one simple embodiment of the novel method and/or system disclosed herein, each of 
the first order models can be provided with a corresponding set of reliability measurements, wherein for 
each set, each measurement therein is a measurement of the reliability of the model (corresponding to 

35 the set) within a particular area of the wireless coverage area. In particular, the reliability measurements 
of each set can be determined by activating each model on archived historical measurements of wireless 
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signals from target mobile stations whose locations have been both verified and also archived. 
Accordingly, when a particular first order model outputs a location hypothesis having an MS location area 
estimate, E, the reliability measurement for the model in an area containing A can be used as a 
confidence value indicating the likelihood of the target MS being in E. Since the confidence values are 
5 probabilities, calculations similar to (or identical to) those used to determine reliability above can be 
performed for determining a measurement of the likelihood that the target MS is in each of the (any) 
intersections of the location hypotheses output by the models. 

It is also worthwhile to mention that further enhancements in reliability can be obtained by 
increasing the reliability of the models R. The above illustration is believed to be conservative in that four 

10 of the first order models F, have relatively poor reliability (i.e., 0.60 or less). Said other way, reliabilities 
exceeding those of at least some the models R are believed to be attainable. For example, since a 
wireless location system according to the novel method and/or system disclosed herein may utilize multiple 
trainable location generating models such as artificial neural networks (ANNs), and since such ANNs are 
commonly trained to perform at levels of 90% or better in many problem domains, it is believed that 

1 5 increases in reliability and/or accuracy of the generated location hypotheses is achievable. Moreover, 
such increases are likely in areas of high population density since in such areas there is likely to be a 
greater number of network base stations detected by and/or detecting a target MS simultaneously. That 
is, there is the capability for increasing the amount of independent wireless signal measurement data that 
can be utilized by such trainable first order models in locating the target MS. 

20 Context Adjuster Introduction. 

The context adjuster (alternatively denoted location adjuster modules) 1 326 module enhances both 
the comparability and predictability of the location hypotheses output by the first order models 1 224. In one 
embodiment (typically where confidence values of location hypotheses are not maintained as probabilities), 
this module modifies location hypotheses received from the FOMs 1224 so that the resulting location 

25 hypotheses output by the context adjuster 1 326 may be further processed uniformly and substantially without 
concern as to differences in accuracy between the first order models from which location hypotheses 
originate.. Further, embodiments of the context adjuster may determine those factors that are perceived to 
impact the perceived accuracy (e.g., confidence) of the location hypotheses:. For instance, environmental 
characteristics may be taken into account here, such as time of day, season, month, weather, geographical 

30 area categorizations (e.g., dense urban, urban, suburban, rural, mountain, etc.), area subcategorizations (e.g., 
heavily treed, hilly, high traffic area, etc.). 

It is also an important and novel aspect of the context adjuster 1326 that the methods for 
adjusting location hypotheses provided in this module may be generalized and thereby also utilized with 
multiple hypothesis computational architectures related to various applications wherein a terrain, surface, 

35 volume or other "geometric" interpretation (e.g., a metric space of statistical samples) may be placed on 
a large body of stored application data for relating hypothesized data to verified data. Moreover, it is 
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important to note that various techniques for "visualizing data" may provide such a geometric 
interpretation. Thus, the methods herein may be utilized in applications such as: 

(a) sonar, radar, x-ray or infrared identification of objects such as occurs in robotic navigation, 
medical image analysis, geological, and radar imaging. 

5 More generally, the novel computational paradigm of the context adjuster 1 326 may be utilized in a 

number of applications wherein there is a large body of archived information providing verified or actual 
application process data related to the past performance of the application process. 

It is worth mentioning that the computational paradigm used in the context adjuster 1 326 is a hybrid 
of a hypothesis adjuster and a data base query mechanism. For example, the context adjuster 1 326 uses an 

l o input (location) hypothesis both as an hypothesis and as a data base query or index into the location 
signature data base 1 320 for constructing a related but more accurate location hypothesis. 

The context adjuster 1 326 tends to create hypotheses that are more accurate than the hypotheses 
generated by the initial hypotheses generators. That is, for each hypothesis, H, provided by one of the initial 
hypothesis generators, G (e.g., a FOM 1224), a corresponding enhanced hypothesis, provided by the context 

1 5 adjuster 1 326, is generated by mapping the past performance of G into the archived verified application data 
(as will be discussed in detail hereinbelow). In particular, the context adjuster hypothesis generation is based 
on the archived verified (or known) performance application data that is related to both G and H. For 
example, in the present wireless location application, if a FOM 1224, G, substantially consistently generates, 
in a particular geographical area, location hypotheses that are biased approximately 1000 feet north of the 

20 actual verified MS 140 location, then the context adjuster 1326 can generate corresponding hypotheses 
without this bias. Thus, the context adjuster 1326 tends to filter out inaccuracies in the initially generated 
hypotheses. 

Therefore in a multiple hypothesis architecture where typically the generated hypotheses may be 
evaluated and/or combined for providing a "most likely" result, it is believed that a plurality of relatively simple 
25 (and possibly inexact) initial hypothesis generators may be used in conjunction with the hybrid computational 
paradigm represented by the context adjuster 1 326 for providing enhanced hypotheses with substantially 
greater accuracy. 

Additionally, note that this hybrid paradigm applies to other domains that are not geographically 
based. For instance, this hybrid paradigm applies to many prediction and/or diagnostic applications for which: 
30 (a) the application data and the application are dependent on a number of parameters whose 

values characterize the range of outputs for the application. That is, there is a set of parameters, pi, p& P3, ... 
, pn from which a parameter space pi x P2 x p3 x ... x pn is derived whose points characterize the actual and 
estimated (or predicted) outcomes. As examples, in the MS location system, pi = latitude and p2 = longitude; 

(b) there is historical data from which points for the parameter space, pi x p 2 x p3 x ... x Pn can 
35 be obtained, wherein this data relates to (or indicates) the performance of the application, and the points 

obtained from this data are relatively dense in the space (at least around the likely future actual outcomes that 
the application is expected to predict or diagnose). For example, such historical data may associate the 
predicted outcomes of the application with corresponding actual outcomes; 
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(c) there is a metric or distance-like evaluation function that can be applied to the parameter 
space for indicating relative closeness or accuracy of points in the parameter space, wherein the evaluation 
function provides a measurement of closeness that is related to the actual performance of the application. 
Note that there are numerous applications for which the above criteria are applicable. For 
5 instance, computer aided control of chemical processing plants are likely to satisfy the above criteria. 
Certain robotic applications may also satisfy this criteria. In fact, it is believed that a wide range of signal 
processing applications satisfy this criteria. 

A second embodiment of the context adjuster will now be described. Referring now to Fig. 19, 
each location hypothesis generated by the first order models 1 224 can be provided to one or more 

l o adjuster modules (Fig. 1 5) for adjusting the generated location hypothesis according to the performance of 
the first order model in the specific localized area "nearby" the newly generated location estimate of the 
location hypothesis. Thus, for partition areas where there is substantial archived verified location data in 
the Location Signature Data Base 1320, such an adjuster module can generate an additional derived 
location estimate. That is, for archived mobile station location estimates previously generated by the first 

15 order model wherein these archived location estimates are also "nearby" or adjacent to the centroid or 
point location estimate of a newly generated location estimate, these archived location estimates are used 
to retrieve their corresponding verified or actual mobile station locations. Subsequently, the retrieved 
verified mobile station locations are used to derive an area that becomes another location estimate of the 
target mobile station, this derived location estimate being denoted herein as an "adjusted image". In 

20 particular, the retrieved locations may be used to generate a convex hull in the coverage area, and 

subsequently cover the convex hull with a minimal covering of cells from a cell mesh that also partitions 
the coverage area (each cell of this mesh being, for example, a 300 ft by 300 ft square). Accordingly, it is 
believed that such adjusted location hypotheses can, in many circumstances, provide enhanced mobile 
station location estimating performance. Moreover, a probability can also be computed for each adjusted 

25 location hypothesis indicating the likelihood that the target mobile station is contained in the location 
estimate of the adjusted location hypothesis. Note that this probability can be computed in a manner 
similar to the computation of the probabilities discussed above for location hypotheses generated directly 
from the first order models. That is, by considering the first order model 1224 and the adjuster module as 
a single location hypothesis generating computational unit, and providing this unit, iteratively, with 

30 archived measurements of wireless transmissions between the base station network and target mobile 
stations at various corresponding archived verified locations within a partition area A, a probability for the 
partition area can be determined as to how likely this computational unit is to generate a location 
hypothesis having a location estimate containing a target mobile station location for the partition area A. 
Accordingly, a table of location partition area probabilities can be determined for each such computational 

35 unit. Thus, when an adjusted location hypothesis is derived for one of the partition areas, the 
corresponding probability may be assigned as the confidence value. 

In Fig. 15, two such adjuster modules are shown, namely, an adjuster for enhancing reliability 
1436 and an adjuster for enhancing accuracy 1440. Both of these adjusters perform their location 
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hypothesis adjustments in the manner described above. The difference between these two adjuster 
modules 1436 and 1440 is primarily the size of the localized area "nearby" the newly generated location 
estimate. In particular, since it is believed that the larger (smaller) the localized nearby area is, the more 
likely (less likely) the corresponding adjusted image is to contain the target mobile station location, the 
5 adjuster for enhancing reliability 1436 may determine its localized areas "nearby" a newly generated 
location estimate as, for example, having a 40% larger diameter (alternatively, area) than the location 
area estimate generated by a first order model 1224. Alternatively, the adjuster for enhancing accuracy 
1444 may determine its localized areas "nearby" a newly generated location estimate as, for example, 
having a 30% smaller diameter (alternatively, area) than the location area estimate generated by a first 
10 order model 1224. Thus, each newly generated location hypothesis can potentially be used to derive at 
least two additional adjusted location hypotheses with some of these adjusted location hypotheses being 
more reliable and some being more accurate than the location hypotheses generated directly from the 
first order models 1224. 

MS Status Repository Introduction 

1 5 The MS status repository 1 338 is a run-time storage manager for storing location hypotheses from 

previous activations of the location engine 1 39 (as well as for storing the output "most likely" target MS 
location estimate(s)) so that a target MS 140 may be tracked using target MS location hypotheses from 
previous location engine 139 activations to determine, for example, a movement of the target MS 140 
between evaluations of the target MS location. 

20 Control and Output Gating Modules 

A fourth functional group of location engine 139 modules is the control and output gating modules 
which includes the location center control subsystem 1350, and the output gateway 1356. The location 
control subsystem 1350 provides the highest level of control and monitoring of the data processing performed 
by the location center 142. In particular, this subsystem performs the following functions: 
25 (a) controls and monitors location estimating processing for each target MS 140. Note that this 

includes high level exception or error handling functions; 

(b) receives and routes external information as necessary. For instance, this subsystem may 

receive (via, e.g., the public telephone switching network 124 and Internet 468) such 
environmental information as increased signal noise in a particular service area due to increase 
30 traffic, a change in weather conditions, a base station 1 22 (or other infrastructure provisioning), 

change in operation status (e.g., operational to inactive); 

(c) receives and directs location processing requests from other location centers 142 (via, e.g., the 

Internet 468); 

(d) performs accounting and billing procedures; 

35 (e) interacts with location center operators by, for example, receiving operator commands and 

providing output indicative of processing resources being utilized and malfunctions; 
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(f) provides access to output requirements for various applications requesting location estimates. 
For example, an Internet location request from a trucking company in Los Angeles to a location 
center 142 in Denver may only want to know if a particular truck or driver is within the Denver 
area. Alternatively, a local medical rescue unit is likely to request as precise a location estimate 
5 as possible. 

Note that in Fig. 6 (a) - (d) above are, at least at a high level, performed by utilizing the operator 
interface 1374. 

Referring now to the output gateway 1356, this module routes target MS 140 location estimates to 
the appropriate location application(s). For instance, upon receiving a location estimate from the most 
i o likelihood estimator 1344, the output gateway 1 356 may determine that the location estimate is for an 
automobile being tracked by the police and therefore must be provided according to a particular protocol. 

System Tuning and Adaptation: The Adaptation Engine 

A fifth functional group of location engine 1 39 modules provides the ability to enhance the MS 
locating reliability and/or accuracy of a wireless location system according to the novel method and/or 

1 5 system disclosed herein by providing it with the capability to adapt to particular operating configurations, 
operating conditions and wireless signaling environments without performing intensive manual analysis of the 
performance of various embodiments of the location engine 139. That is, this functional group automatically 
enhances the performance of the location engine for locating MSs 140 within a particular coverage area 120 
using at least one wireless network infrastructure therein. More precisely, this functional group allows the 

20 wireless location system to adapt by tuning or optimizing certain system parameters according to location 
engine 139 location estimate accuracy and reliability. 

Implementations of First Order Models 

Further descriptions of various first order models 1 224 are provided in this section. 
Distance First Order Models (TOA/TDOA) 

25 As discussed in the Location Center Architecture Overview section herein above, distance 

models determine a presumed direction and/or distance that a target MS 140 is from one or more base 
stations 122. In some embodiments of distance models, the target MS location estimate(s) generated are 
obtained using radio signal analysis techniques that are quite general and therefore are not capable of 
taking into account the peculiarities of the topography of a particular radio coverage area. For example, 

30 substantially all radio signal analysis techniques using conventional procedures (or formulas) are based 
on "signal characteristic measurements" such as: 

(a) signal timing measurements (e.g., TOA and TDOA), 

(b) signal strength measurements, and/or 

(c) signal angle of arrival measurements. 
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Furthermore, such signal analysis techniques are likely predicated on certain very general assumptions 
that can not fully account for signal attenuation and multipath due to a particular radio coverage area 
topography. 

Taking a CDMA or TDMA base station network as an example, each base station (BS) 122 is 
5 required to emit a constant signal-strength pilot channel pseudo-noise (PN) sequence on the forward link 
channel identified uniquely in the network by a pilot sequence offset and frequency assignment. It is 
possible to use the pilot channels of the active, candidate, neighboring and remaining sets, maintained in 
the target MS, for obtaining signal characteristic measurements (e.g., TOA and/or TDOA measurements) 
between the target MS 140 and the base stations in one or more of these sets. 

l o Based on such signal characteristic measurements and the speed of signal propagation, signal 

characteristic ranges or range differences related to the location of the target MS 140 can be calculated. 
Using TOA and/or TDOA ranges as exemplary, these ranges can then be input to either the radius-radius 
multilateration or the time difference multilateration algorithms along with the known positions of the 
corresponding base stations 122 to thereby obtain one or more location estimates of the target MS 140. 

1 5 For example, if there are, four base stations 1 22 in the active set, the target MS 1 40 may cooperate with 
each of the base stations in this set to provide signal arrival time measurements. Accordingly, each of the 
resulting four sets of three of these base stations 122 may be used to provide an estimate of the target 
MS 140 as one skilled in the art will understand. Thus, potentially (assuming the measurements for each 
set of three base stations yields a feasible location solution) there are four estimates for the location of the 

20 target MS 140. Further, since such measurements and BS 122 positions can be sent either to the network 
or the target MS 140, location can be determined in either entity. 

Since many of the signal measurements utilized by embodiments of distance models are subject 
to signal attenuation and multipath due to a particular area topography. Many of the sets of base stations 
from which target MS location estimates are desired may result in either no location estimate, or an 

25 inaccurate location estimate. 

Accordingly, some embodiments of distance FOMs may attempt to mitigate such ambiguity or 
inaccuracies by, e.g., identifying discrepancies (or consistencies) between arrival time measurements and 
other measurements (e.g., signal strength), these discrepancies (or consistencies) may be used to filter 
out at least those signal measurements and/or generated location estimates that appear less accurate. In 

30 particular, such identifying by filtering can be performed by, for example, an expert system residing in the 
distance FOM. 

A second approach for mitigating such ambiguity or conflicting MS location estimates is 
particularly novel in that each of the target MS location estimates is used to generate a location 
hypothesis regardless of its apparent accuracy. Accordingly, these location hypotheses are input to an 
35 alternative embodiment of the context adjuster 1 326 that is substantially (but not identical to) the context 
adjuster as described in detail in Appendix D of U.S. Patent Application Serial No. 09/194,367 (from 
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which the present patent application claims priority), so that each location hypothesis may be adjusted to 
enhance its accuracy. In contradistinction to the embodiment of the context adjuster 1326 of APPENDIX 
D of U.S. Patent Application Serial No. 09/194,367, where each location hypothesis is adjusted according 
to past performance of its generating FOM 1224 in an area of the initial location estimate of the location 
5 hypothesis (the area, e.g., determined as a function of distance from this initial location estimate), this 
alternative embodiment adjusts each of the location hypotheses generated by a distance first order model 
according to a past performance of the model as applied to signal characteristic measurements from the 
same set of base stations 122 as were used in generating the location hypothesis. That is, instead of 
only using only an identification of the distance model (i.e., its FOMJD) to, for example, retrieve archived 

l o location estimates generated by the model in an area of the location hypothesis' estimate (when 

determining the model's past performance), the retrieval retrieves only the archived location estimates 
that are, in addition, derived from the signal characteristics measurement obtained from the same 
collection of base stations 1 22 as was used in generating the location hypothesis. Thus, the adjustment 
performed by this embodiment of the context adjuster 1 326 adjusts according to the past performance of 

15 the distance model and the collection of base stations 122 used. 

Stochastic First Order Model 

The stochastic first order models may use statistical prediction techniques such as principle 
decomposition, partial least squares, or other regression techniques for predicting, for example, expected 
minimum and maximum distances of the target MS from one or more base stations 1 22, e.g., Bollenger 
20 Bands. Additionally, some embodiments may use Markov processes and Random Walks (predicted 
incremental MS movement) for determining an expected area within which the target MS 140 is likely to. 
be. That is, such a process measures the incremental time differences of each pilot as the MS moves for 
predicting a size of a location area estimate using past MS estimates such as the verified location 
signatures in the location signature data base 1320. 

25 Pattern Recognition and Adaptive First Order Models 

It is a particularly important aspect of a wireless location system according to the novel method 
and/or system disclosed herein to provide: 

(a) one or more FOMs 1224 that generate target MS 140 location estimates by using pattern 
recognition or associativity techniques, and/or 
30 (b) one or more FOMs 1 224 that are adaptive or trainable so that such FOMs may generate 

increasingly more accurate target MS location estimates from additional training. 

Statistically Based Pattern Recognition First Order Models 

Regarding FOMs 1224 using pattern recognition or associativity techniques, there are many 
such techniques available. For example, there are statistically based systems such as "CART" (an 
35 acronym for Classification and Regression Trees) by ANGOSS Software International Limited of Toronto, 
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Canada that may be used for automatically detecting or recognizing patterns in data that were unprovided 
(and likely previously unknown). Accordingly, by imposing a relatively fine mesh or grid of cells on the 
radio coverage area, wherein each cell is entirely within a particular area type categorization such as the 
transmission area types (discussed in the section, "Coverage Area: Area Types And Their Determination" 
5 above), the verified location signature clusters within the cells of each area type may be analyzed for 
signal characteristic patterns. If such patterns are found, then they can be used to identify at least a likely 
area type in which a target MS is likely to be located. That is, one or more location hypotheses may be 
generated having target MS 140 location estimates that cover an area having the likely area type wherein 
the target MS 140 is located. Further note that such statistically based pattern recognition systems as 

10 "CART" include software code generators for generating expert system software embodiments for 
recognizing the patterns detected within a training set (e.g., the verified location signature clusters). 

Accordingly, although an embodiment of a FOM as described here may not be exceedingly 
accurate, it may be very reliable. Thus, since a fundamental aspect of the novel method and/or system 
disclosed herein is to use a plurality MS location techniques for generating location estimates and to 

15 analyze the generated estimates (likely after being adjusted) to detect patterns of convergence or 

clustering among the estimates, even large MS location area estimates are useful. For example, it can be 
the case that four different and relatively large MS location estimates, each having very high reliability, 
have an area of intersection that is acceptably precise and inherits the very high reliability from each of 
the large MS location estimates from which the intersection area was derived. 

20 Adaptive/Trainable First Order Models 

Adaptive/Trainable First Order Models 

The term adaptive is used to describe a data processing component that can modify its data 
processing behavior in response to certain inputs that are used to change how subsequent inputs are 
processed by the component. Accordingly, a data processing component may be "explicitly adaptive" by 

25 modifying its behavior according to the input of explicit instructions or control data that is input for 

changing the component's subsequent behavior in ways that are predictable and expected. That is, the 
input encodes explicit instructions that are known by a user of the component. Alternatively, a data 
processing component may be "implicitly adaptive" in that its behavior is modified by other than 
instructions or control data whose meaning is known by a user of the component. For example, such 

30 implicitly adaptive data processors may learn by training on examples, by substantially unguided 

exploration of a solution space, or other data driven adaptive strategies such as statistically generated 
decision trees. Accordingly, it is an aspect of the novel method and/or system disclosed herein to utilize 
not only explicitly adaptive MS location estimators within FOMs 1224, but also implicitly adaptive MS 
location estimators. In particular, artificial neural networks (also denoted neural nets and ANNs herein) 

35 are used in some embodiments as implicitly adaptive MS location estimators within FOMs. Thus, in the 
sections below, neural net architectures and their application to locating an MS is described. 
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Artificial Neural Networks For MS Location 

Artificial neural networks may be particularly useful in developing one or more first order models 
1 224 for locating an MS 140, since, for example, ANNs can be trained for classifying and/or associatively 
pattern matching of various RF signal measurements such as the location signatures. That is, by training 
5 one or more artificial neural nets using RF signal measurements from verified locations so that RF signal 
transmissions characteristics indicative of particular locations are associated with their corresponding 
locations, such trained artificial neural nets can be used to provide additional target MS 140 location 
hypotheses. Moreover, it is an aspect of the novel method and/or system disclosed herein that the 
training of such artificial neural net based FOMs (ANN FOMs) is provided without manual intervention as 
l o will be discussed hereinbelow. 

Detailed Description of the Hypothesis Evaluator Modules 

Context Adjuster Embodiments 

The context adjuster 1326 performs the first set of potentially many adjustments to at least the 
confidences of location hypotheses, and in some important embodiments, both the confidences and the target 

1 5 MS location estimates provided by FOMs 1 224 may be adjusted according to previous performances of the 
FOMs. More particularly, as mentioned above, the context adjuster adjusts confidences so that, assuming 
there is a sufficient density verified location signature clusters captured in the location signature data base 
1320, the resulting location hypotheses output by the context adjuster 1326 may be further processed 
uniformly and substantially without concern as to differences in accuracy between the first order models from 

20 which location hypotheses originate. Accordingly, the context adjuster adjusts location hypotheses both to 
environmental factors (e.g., terrain, traffic, time of day, etc., as described in 30.1 above), and to how 
predictable or consistent each first order model (FOM) has been at locating previous target MS's whose 
locations were subsequently verified. 

Of particular importance is the novel computational paradigm utilized herein. That is, if there is a 

25 sufficient density of previous verified MS location data stored in the location signature data base 1 320, then 
the FOM location hypotheses are used as an "index" into this data base (i.e., the location signature data base) 
for constructing new target MS 1 40 location estimates. A more detailed discussion of this aspect of the novel 
method and/or system disclosed herein is given hereinbelow. Accordingly, only a brief overview is provided 
here. Thus, since the location signature data base 1320 stores previously captured MS location data 

30 including: 

(a) clusters of MS location signature signals (see the location signature data base section for a 
discussion of these signals) and 

(b) a corresponding verified MS location, for each such cluster, from where the MS signals 
originated, 

35 the context adjuster 1 326 uses newly created target MS location hypotheses output by the FOMs as indexes 

or pointers into the location signature data base for identifying other geographical areas where the target MS 

140 is likely to be located based on the verified MS location data in the location signature data base. 
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In particular, at least the following two criteria are addressed by the context adjuster 1326: 

(32.1 ) Confidence values for location hypotheses are to be comparable regardless of first order 

models from which the location hypotheses originate. That is, the context adjuster 
moderates or dampens confidence value assignment distinctions or variations between first 
order models so that the higher the confidence of a location hypothesis, the more likely (or 
unlikely, if the location hypothesis indicates an area estimate where the target MS is NOT) 
the target MS is perceived to be in the estimated area of the location hypothesis regardless 
of the First Order Model from which the location hypothesis was output; 

(32.2) Confidence values for location hypotheses may be adjusted to account for current 

environmental characteristics such as month, day (weekday or weekend), time of day, 
area type (urban, rural, etc.), traffic and/or weather when comparing how accurate the 
first order models have previously been in determining an MS location according to such 
environmental characteristics. For example, in one embodiment, such environmental 
characteristics are accounted for by utilizing a transmission area type scheme (as 
discussed in section 5.9 above) when adjusting confidence values of location 
hypotheses. 

Note that in satisfying the above two criteria, the context adjuster 1326, at least in one embodiment, 
may use heuristic (fuzzy logic) rules to adjust the confidence values of location hypotheses from the first 
order models. Additionally, the context adjuster may also satisfy the following criteria: 

(33.1) The context adjuster may adjust location hypothesis confidences due to BS failure(s), 

(33.2) Additionally in one embodiment, the context adjuster may have a calibration mode for at 

least one of: 

(a) calibrating the confidence values assigned by first order models to their location 
hypotheses outputs; 

(b) calibrating itself. 

Additional details of embodiments of the context adjuster 1326 can be found in the publication, 
"Wireless Location System," by Dupray and Karr having International Patent Application No. 
PCT/US97/15892 filed September 8, 1997, and also having a corresponding U.S. national Patent 
Application No. 09/194,367 filed November 24, 1998, wherein both of these documents are fully 
incorporated herein by reference. 

Location Hypothesis Analyzer Embodiment 

The location hypothesis analyzer is primarily for embodiments of a wireless location system 
according to the novel method and/or system disclosed herein where the location hypothesis confidence 
values are not probabilities, and may therefore be modified easily. Referring now to Fig. 7, an embodiment of 
the Hypothesis Analyzer is illustrated. The control component is denoted the control module 1400 . Thus, 
this control module manages or controls access to the run time location hypothesis storage area 1410. The 
control module 1400 and the run time location hypothesis storage area 1410 may be implemented as a 
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blackboard system and/or an expert system. Accordingly, in the blackboard embodiment, the control module 
1400 determines when new location hypotheses may be entered onto the blackboard from other processes 
such as the context adjuster 1326 as well as when location hypotheses may be output to the most likelihood 
estimator 1344. 

The following is a brief description of each submodule included in the location hypothesis analyzer 1 332. 

(35.1 ) A control module 1400 for managing or controlling further processing of location hypotheses received 
from the context adjuster. This module controls all location hypothesis processing within the location 
hypothesis analyzer as well as providing the input interface with the context adjuster. There are 
numerous embodiments that may be utilized for this module, including, but not limited to, expert 
systems and blackboard managers. 

(35.2) A run-time location hypothesis storage area 1410 for retaining location hypotheses during their 
processing by the location hypotheses analyzer. This can be, for example, an expert system fact base 
or a blackboard. Note that in some of the discussion hereinbelow, for simplicity, this module is referred 
to as a "blackboard". However, it is not intended that such notation be a limitation; i.e., the term 
"blackboard" hereinafter will denote a run-time data repository for a data processing paradigm wherein 
the flow of control is substantially data-driven. 

(35.3) An analytical reasoner module 1416 for determining if (or how well) location hypotheses are consistent 
with well known physical or heuristic constraints as, e.g., mentioned in (30.4) above. Note that this 
module may be a daemon or expert system rule base. 

(35.4) An historical location reasoner module 1424 for adjusting location hypotheses' confidences according 
to how well the location signature characteristics (i.e., loc sigs) associated with a location hypothesis 
compare with "nearby" loc sigs in the location signature data base as indicated in (30.3) above. Note 
that this module may also be a daemon or expert system rule base. 

(35.5) A location extrapolator module 1432 for use in updating previous location estimates for a target MS 
when a more recent location hypothesis is provided to the location hypothesis analyzer 1 332. That is, 
assume that the control module 1 400 receives a new location hypothesis for a target MS for which 
there are also one or more previous location hypotheses that either have been recently processed (i.e., 
they reside in the MS status repository 1338, as shown best in Fig. 6), or are currently being processed 
(i.e., they reside in the run-time location hypothesis storage area 1410). Accordingly, if the 
activejimestamp (see Figs. 9A, B regarding location hypothesis data fields) of the newly received 
location hypothesis is sufficiently more recent than the activejimestamp of one of these previous 
location hypotheses, then an extrapolation may be performed by the location extrapolator module 1432 
on such previous location hypotheses so that all target MS location hypotheses being concurrently 
analyzed are presumed to include target MS location estimates for substantially the same point in time. 
Thus, initial location estimates generated by the FOMs using different wireless signal measurements, 
from different signal transmission time intervals, may have their corresponding dependent location 
hypotheses utilized simultaneously for determining a most likely target MS location estimate. Note that 
this module may also be daemon or expert system rule base. 
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(35.6) hypothesis generating module 1428 for generating additional location hypotheses according to, for 
example, MS location information not adequately utilized or modeled. Note, location hypotheses may 
also be decomposed here if, for example it is determined that a location hypothesis includes an MS 
area estimate that has subareas with radically different characteristics such as an MS area estimate 
that includes an uninhabited area and a densely populated area. Additionally, the hypothesis 
generating module 1428 may generate "poor reception" location hypotheses that specify MS location 
areas of known poor reception that are "near" or intersect currently active location hypotheses. Note, 
that these poor reception location hypotheses may be specially tagged (e.g., with a distinctive FOM JD 
value or specific tag field) so that regardless of substantially any other location hypothesis confidence 
value overlapping such a poor reception area, such an area will maintain a confidence value of 
"unknown" (i.e., zero). Note that substantially the only exception to this constraint is location 
hypotheses generated from mobile base stations 148. Note that this module may also be daemon or 
expert system rule base. 

In the blackboard system embodiment of the location hypothesis analyzer, a blackboard system is 
the mechanism by which the last adjustments are performed on location hypotheses and by which additional: 
location hypotheses may be generated.. Briefly, a blackboard system can be described as a particular class 
of software that typically includes at least three basic components. That is: 

(36.1 ) a data base called the "blackboard," whose stored information is commonly available to a 
collection of programming elements known as "daemons", wherein, in a wireless location 
system according to the novel method and/or system disclosed herein, the blackboard 
includes information concerning the current status of the location hypotheses being 
evaluated to determine a "most likely" MS location estimate. Note that this data base is 
provided by the run time location hypothesis storage area 1410; 

(36.2) one or more active (and typically opportunistic) knowledge sources, denoted conventionally 
as "daemons," that create and modify the contents of the blackboard. The blackboard 
system employed requires only that the daemons have application knowledge specific to 
the MS location problem addressed by a wireless location system according to the novel 
method and/or system disclosed herein. As shown in Fig. 7, the knowledge sources or 
daemons in the hypothesis analyzer include the analytical reasoner module 1 41 6, the 
hypothesis generating module 1428, and the historical location reasoner module 1416; 

(36.3) a control module that enables the realization of the behavior in a serial computing 
environment. The control element orchestrates the flow of control between the various 
daemons. This control module is provided by the control module 1400. 

Note that this blackboard system may be commercially available, however, the knowledge sources, i.e., 
daemons, may be developed specifically for an embodiment of the novel method and/or system disclosed 
herein. For further information regarding such blackboard systems, the following references are incorporated 
herein by reference: (a) Jagannathan, V,, Dodhiawala, R, & Baum, L. S. (1989). Blackboard architectures 
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and applications. Boston, MA: Harcourt Brace Jovanovich Publishers; (b) Engelmore, R., & Morgan, T. 
(1988). Blackboard systems. Reading, MA: Addison-Wesley Publishing Company. 

Alternatively, the control module 1400 and the run-time location hypothesis storage area 1410 may be 
implemented as an expert system or as a fuzzy rule inferencing system, wherein the control module 1400 

5 activates or 'fires" rules related to the knowledge domain (in the present case, rules relating to the accuracy of 
MS location hypothesis estimates), and wherein the rules provide a computational embodiment of, for 
example, constraints and heuristics related to the accuracy of MS location estimates. Thus, the control 
module 1400 for the present embodiment is also used for orchestrating, coordinating and controlling the 
activity of the individual rule bases of the location hypothesis analyzer (e.g. as shown in Fig. 7, the analytical 

10 reasoner module 1416, the hypothesis generating module 1428 , the historical location reasoner module 

1424, and the location extrapolator module 1432). For further information regarding such expert systems, the 
following reference is incorporated herein by reference: Waterman, D. A. (1970). A guide to expert systems. 
Reading, MA: Addison-Wesley Publishing Company. 

MS Status Repository Embodiment 

1 5 The MS status repository 1 338 is a run-time storage manager for storing location hypotheses from 

previous activations of the location engine 139 (as well as the output target MS location estimate(s)) so 
that a target MS may be tracked using target MS location hypotheses from previous location engine 139 
activations to determine, for example, a movement of the target MS between evaluations of the target MS 
location. Thus, by retaining a moving window of previous location hypotheses used in evaluating 

20 positions of a target MS, measurements of the target MS's velocity, acceleration, and likely next position 
may be determined by the location hypothesis analyzer 1332. Further, by providing accessibility to recent 
MS location hypotheses, these hypotheses may be used to resolve conflicts between hypotheses in a 
current activation for locating the target MS; e.g., MS paths may be stored here for use in extrapolating a 
new location 

25 Most Likelihood Estimator Embodiment 

The most likelihood estimator 1 344 is a module for determining a "most likely" location estimate for a 
target MS 140 being located (e.g., as in (30.7) above). In one embodiment, the most likelihood estimator 
performs an integration or summing of all location hypothesis confidence values for any geographic region(s) 
of interest having at least one location hypothesis that has been provided to the most likelihood estimator, and 

30 wherein the location hypothesis has a relatively (or sufficiently) high confidence. That is, the most likelihood 
estimator 1 344 determines the area(s) within each such region having high confidences (or confidences 
above a threshold) as the most likely target MS 140 location estimates. 

In one embodiment of the most likelihood estimator 1344, this module utilizes an area mesh, M, over 
which to integrate, wherein the mesh cells of M are preferably smaller than the greatest location accuracy 

35 desired. That is, each cell, c, of M is assigned a confidence score indicating a likelihood that the target MS 
140 is located in c, wherein the confidence score for c is determined by the confidence values of the target 
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MS location estimates provided to the most likelihood estimator 1344. Thus, to obtain the most likely location 
determination(s) the following steps are performed: 

(a) For each of the active location hypotheses output by, e.g., the hypothesis analyzer 1332 

(alternatively, the context adjuster 1326), each corresponding MS location area estimate, LAE, 
is provided with a smallest covering, Clea, of cells c from M. 

(b) Subsequently, assuming confidence values are probabilities, adjust the confidence score of each 

of the cells of Clea by adding to it the confidence value for LAE divided by the number of cells in 
Clea. Also, for each cell outside of Clea and within some reasonable larger area A (e.g., the 
intersection of the coverage areas of each base station detected by the target MS 1 40), add: 
(1 - (the confidence value of Clea)) divided by (the number of cells in A - Clea) 

(c) Given that the confidence values are within a particular range (e.g., 0 to 1 when the confidence 
values are probabilities), and that this range has been partitioned into intervals, Int, having 
lengths of, e.g., 0.05, for each interval, Int, perform a cluster analysis function for clustering cells 
with confidences that are in Int. Thus, a topographical-type map may be constructed from the 
resulting cell clusters, wherein higher confidence areas are analogous to representations of 
areas having higher elevations. 

(d) Output a representation of the resulting clusters for each Int to the output gateway 1356 for 

determining the location granularity and representation desired by each location application 146 
requesting the location of the target MS 140. 

Of course, variations in the above algorithm are also within the scope of the novel method and/or 
system disclosed herein. For example, some embodiments of the most likelihood estimator 1344 may: 

(e) Perform special processing for areas designated as "poor reception" areas. For example, the 

most likelihood estimator 1 344 may be able to impose a confidence value of zero (i.e., meaning 
it is unknown as to whether the target MS is in the area) on each such poor reception area 
regardless of the location estimate confidence values unless there is a location hypothesis from 
a reliable and unanticipated source. That is, the mesh cells of a poor reception area may have 
their confidences set to zero unless, e.g., there is a location hypothesis derived from target MS 
location data provided by a mobile base station 148 that: (a) is near the poor reception area, (b) 
able to detect that the target MS 140 is in the poor reception area, and (c) can relay target MS 
location data to the location center 1 42. In such a case, the confidence of the target MS 
location estimate from the MBS location hypothesis may take precedence. 

(f) Additionally, in some embodiments of the most likelihood estimator 1344, cells c of M that are 

"near" or adjacent to a covering Clea may also have their confidences adjusted according to 
how near the cells c are to the covering. That is, the assigning of confidences to cell meshes 
may be "fuzzified" in the terms of fuzzy logic so that the confidence value of each location 
hypothesis utilized by the most likelihood estimator 1344 is provided with a weighting factor 
depending on its proxity to the target MS location estimate of the location hypothesis. More 
precisely, it is believed that "nearness," in the present context, should be monotonic with the 
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"wideness" of the covering; i.e., as the extent of the covering increases (decreases) in a 
particular direction, the cells c affected beyond the covering also increases (decreases). 
Furthermore, in some embodiments of the most likelihood estimator 1344, the greater (lesser) 
the confidence in the LEA, the more (fewer) cells c beyond the covering have their confidences 
affected. To describe this technique in further detail, reference is made to Fig. 1 0, wherein an 
area A is assumed to be a covering Clea having a confidence denoted "conf. Accordingly, to 
determine a confidence adjustment to add to a cell c not in A (and additionally, the centroid of A 
not being substantially identical with the centroid of c which could occur if A were donut 
shaped), the following steps may be performed: 

(i) Determine the centroid of A, denoted Cent(A). 

(ii) Determine the centroid of the cell c ( denoted Q. 

(iii) Determine the extent of A along the line between Cent(A) and Q, denoted L 

(iv) For a given type of probability density function, P(x), such as a Gaussian function, 
let T be the beginning portion of the function that lives on the x-axis interval [0, t], 
wherein P(t) = ABS(conf) = the absolute value of the confidence of Clea. 

(v) Stretch T along the x-axis so that the stretched function, denoted sT(x), has an x- 
axis support of [0, UiUeW^**®'^)], where a is in range of 3.0 to 10.0; e.g., 5.0. 
Note that sT(x) is the function, 

P(x * (1 +e-ta(ABS(oonf y on this stretched extent. Further note that for 
confidences of +1 and -1 , the support of sT(x) is [0, L] and for confidences at (or 
near) zero this support. Further, the term, 

L/(1 +e {a(ABS{cortO-1)]) 

is monotonically increasing with L and ABS(conf). 

(vi) Determine D = the minimum distance that Q is outside of A along the line 
between Cent(A) and Q, 

(vii) Determine the absolute value of the change in the confidence of c as sT(D). 

(viii) Provide the value sT(D) with the same sign as conf, and provide the potentially 
sign changed value sT(D) as the confidence of the cell c. 

Additionally, in some embodiments, the most likelihood estimator 1344, upon receiving one or more 
location hypotheses from the hypothesis analyzer 1 332, also performs some or all of the following tasks: 

(37.1) Determines the area of interest over which to perform the integration. In one embodiment, 
this area is a convex hull including each of the MS area estimates from the received location 
hypotheses (wherein such location hypotheses have not been removed from consideration 
by the filtering process of (37.1)); 

(37.2) Determines, once the integration is performed, one or more collections of contiguous area 

mesh cells that may be deemed a "most likely" MS location estimate, wherein each such 
collection includes one or more area mesh cells having a high confidence value. 
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Hypothesis Generating Module 

The hypothesis generating module 1428 (Fig. 7) is used for generating additional location hypotheses 
according to, for example, MS location information not adequately utilized or modeled. Note, location 
hypotheses may also be decomposed here if, for example it is determined that a location hypothesis includes 
5 an MS area estimate that has subareas with radically different characteristics such as an area that includes an 
uninhabited area and a densely populated area. Additionally, the hypothesis generating module 1428 may 
generate "poor reception" location hypotheses that specify MS location areas of known poor reception that are 
"near" or intersect currently active location hypotheses. Note, that these poor reception location hypotheses 
may be specially tagged (e.g., with a distinctive FOMJD value or specific tag field) so that regardless of 
i o substantially any other location hypothesis confidence value overlapping such a poor reception area, such an 
area will mantain a confidence value of "unknown" (i.e., zero). Note that substantially the only exception to 
this constraint is location hypotheses generated from mobile base stations 148. 

Mobile Base Station Location Subsystem Description 

Mobile Base Station Subsystem Introduction 

1 5 Any collection of mobile electronics (denoted mobile location unit) that is able to both estimate a 

location of a target MS 140 and communicate with the base station network may be utilized by a wireless 
location system according to the novel method and/or system disclosed herein to more accurately locate 
the target MS. Such mobile location units may provide greater target MS location accuracy by, for 
example, homing in on the target MS and by transmitting additional MS location information to the 

20 location center 142. There are a number of embodiments for such a mobile location unit that are within 
the scope of a wireless location system according to the novel method and/or system. For example, in a 
minimal version, such the electronics of the mobile location unit may be little more than an onboard MS 
140, a sectored/directional antenna and a controller for communicating between them. Thus, the onboard 
MS is used to communicate with the location center 142 and possibly the target MS 140, while the 

25 antenna monitors signals for homing in on the target MS 140. In an enhanced version of the mobile 

location unit, a GPS receiver may also be incorporated so that the location of the mobile location unit may 
be determined and consequently an estimate of the location of the target MS may also be determined. 
However, such a mobile location unit is unlikely to be able to determine substantially more than a 
direction of the target MS 140 via the sectored/directional antenna without further base station 

30 infrastructure cooperation in, for example, determining the transmission power level of the target MS or 
varying this power level. Thus, if the target MS or the mobile location unit leaves the coverage area 120 
or resides in a poor communication area, it may be difficult to accurately determine where the target MS is 
located. None-the-less, such mobile location units may be sufficient for many situations, and in fact the 
wireless location system disclosed herein contemplates their use. However, in cases where direct 

35 communication with the target MS is desired without constant contact with the base station infrastructure, 
the wireless location system may include a mobile location unit that is also a scaled down version of a 
base station 122. Thus, given that such a mobile base station or MBS 148 includes at least an onboard 
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MS 140, a sectored/directional antenna, a GPS receiver, a scaled down base station 122 and sufficient 
components (including a controller) for integrating the capabilities of these devices, an enhanced 
autonomous MS mobile location system can be provided that can be effectively used in, for example, 
emergency vehicles, air planes and boats. Accordingly, the description that follows below describes an 
5 embodiment of an MBS 148 having the above mentioned components and capabilities for use in a 
vehicle. 

As a consequence of the MBS 148 being mobile, there are fundamental differences in the 
operation of an MBS in comparison to other types of BS's 1 22 (1 52). In particular, other types of base 
stations have fixed locations that are precisely determined and known by the location center, whereas a 

10 location of an MBS 148 may be known only approximately and thus may require repeated and frequent 
re-estimating. Secondly, other types of base stations have substantially fixed and stable communication 
with the location center (via possibly other BS's in the case of LBSs 152) and therefore although these 
BS's may be more reliable in their in their ability to communicate information related to the location of a 
target MS with the location center, accuracy can be problematic in poor reception areas. Thus, MBS's 

1 5 may be used in areas (such as wilderness areas) where there may be no other means for reliably and 
cost effectively locating a target MS 140 (i.e., there may be insufficient fixed location BS's coverage in an 
area). 

Fig. 1 1 provides a high level block diagram architecture of one embodiment of the MBS location 
subsystem 1508. Accordingly, an MBS may include components for communicating with the fixed 

20 location BS network infrastructure and the location center 142 via an on-board transceiver 1512 that is 
effectively an MS 140 integrated into the location subsystem 1508. Thus, if the MBS 148 travels through 
an area having poor infrastructure signal coverage, then the MBS may not be able to communicate 
reliably with the location center 142 (e.g., in rural or mountainous areas having reduced wireless 
telephony coverage). So it is desirable that the MBS 148 must be capable of functioning substantially 

25 autonomously from the location center. In one embodiment, this implies that each MBS 148 must be 
capable of estimating both its own location as well as the location of a target MS 140. 

Additionally, many commercial wireless telephony technologies require all BS's in a network to 
be very accurately time synchronized both for transmitting MS voice communication as well as for other 
services such as MS location. Accordingly, the MBS 148 will also require such time synchronization. 

30 However, since an MBS 148 may not be in constant communication with the fixed location BS network 
(and indeed may be off-line for substantial periods of time), on-board highly accurate timing device may 
be necessary. In one embodiment, such a device may be a commercially available ribidium oscillator 
1520 as shown in Fig. 11. 

Since the MBS 148 , includes a scaled down version of a BS 122 (denoted 1522 in Fig. 1 1), it is 

35 capable of performing most typical BS 1 22 tasks, albeit on a reduced scale. In particular, the base station 
portion of the MBS 148 can: 

(a) raise/lower its pilot channel signal strength, 

(b) be in a state of soft hand-off with an MS 140, and/or 



(c) be the primary BS 1 22 for an MS 140, and consequently be in voice communication with the 
target MS (via the MBS operator telephony interface 1524) if the MS supports voice 
communication. 

Further, the MBS 148 can, if it becomes the primary base station communicating with the MS 140, 
5 request the MS to raise/lower its power or, more generally, control the communication with the MS (via 
the base station components 1522). However, since the MBS 148 will likely have substantially reduced 
telephony traffic capacity in comparison to a standard infrastructure base station 122, note that the pilot 
channel for the MBS is preferably a nonstandard pilot channel in that it should not be identified as a 
conventional telephony traffic bearing BS 122 by MS's seeking normal telephony communication. Thus, a 

l o target MS 140 requesting to be located may, depending on its capabilities, either automatically configure 
itself to scan for certain predetermined MBS pilot channels, or be instructed via the fixed location base 
station network (equivalently BS infrastructure) to scan for a certain predetermined MBS pilot channel. 

Moreover, the MBS 148 has an additional advantage in that it can substantially increase the reliability 
of communication with a target MS 140 in comparison to the base station infrastructure by being able to 

1 5 move toward or track the target MS 140 even if this MS is in (or moves into) a reduced infrastructure base 
station network coverage area. Furthermore, an MBS 148 may preferably use a directional or smart 
antenna 1526 to more accurately locate a direction of signals from a target MS 140. Thus, the sweeping 
of such a smart antenna 1 526 (physically or electronically) provides directional information regarding 
signals received from the target MS 140. That is, such directional information is determined by the signal 

20 propagation delay of signals from the target MS 140 to the angular sectors of one of more directional 
antennas 1526 on-board the MBS 148. 

Before proceeding to further details of the MBS location subsystem 1508, an example of the 
operation of an MBS 148 in the context of responding to a 91 1 emergency call is given. In particular, this 
example describes the high level computational states through which the MBS 148 transitions, these 

25 states also being illustrated in the state transition diagram of Fig. 12. Note that this figure illustrates the 
primary state transitions between these MBS 148 states, wherein the solid state transitions are indicative 
of a typical "ideal" progression when locating or tracking a target MS 140, and the dashed state transitions 
are the primary state reversions due, for example, to difficulties in locating the target MS 140. 

Accordingly, initially the MBS 148 may be in an inactive state 1700, wherein the MBS location 

30 subsystem 1508 is effectively available for voice or data communication with the fixed location base 
station network but the MS 140 locating capabilities of the MBS are not active. From the inactive state 
1700 the MBS (e.g., a police or rescue vehicle) may enter an active state 1704 once an MBS operator 
has logged onto the MBS location subsystem of the MBS, such logging being for authentication, 
verification and journaling of MBS 148 events. In the active state 1704, the MBS may be listed by a 91 1 

35 emergency center and/or the location center 142 as eligible for service in responding to a 91 1 request. 
From this state, the MBS 148 may transition to a ready state 1708 signifying that the MBS is ready for use 
in locating and/or intercepting a target MS 140. That is, the MBS 148 may transition to the ready state 
1708 by performing the following steps: 



(1 a) Synchronizing the timing of the location subsystem 1508 with that of the base station network 
infrastructure. In one embodiment, when requesting such time synchronization from the base 
station infrastructure, the MBS 148 will be at a predetermined or well known location so that the 
MBS time synchronization may adjust for a known amount of signal propagation delay in the 
5 synchronization signal. 

(1b) Establishing the location of the MBS 148. In one embodiment, this may be accomplished by, for 
example, an MBS operator identifying the predetermined or well known location at which the 
MBS 148 is located. 

(1c) Communicating with, for example, the 91 1 emergency center via the fixed location base station 
10 infrastructure to identify the MBS 148 as in the ready state. 

Thus, while in the ready state 1708, as the MBS 148 moves, it has its location repeatedly (re)- 
estimated via, for example, GPS signals, location center 142 location estimates from the base stations 
122 (and 152), and an on-board deadreckoning subsystem 1527 having an MBS location estimator 
according to the programs described hereinbelow. However, note that the accuracy of the base station 
1 5 time synchronization (via the ribidium oscillator 1 520) and the accuracy of the MBS 1 48 location may 
need to both be periodically recalibrated according to (1a) and (1b) above. 

Assuming a 911 signal is transmitted by a target MS 140, this signal is transmitted, via the fixed 
location base station infrastructure, to the 91 1 emergency center and the location center 142, and 
assuming the MBS 148 is in the ready state 1708, if a corresponding 91 1 emergency request is 
20 transmitted to the MBS (via the base station infrastructure) from the 91 1 emergency center or the location 
center, then the MBS may transition to a seek state 1712 by performing the following steps: 

(2a) Communicating with, for example, the 91 1 emergency response center via the fixed location 
base station network to receive the PN code for the target MS to be located (wherein this 
communication is performed using the MS-like transceiver 1 51 2 and/or the MBS operator 
25 telephony interface 1524). 

(2b) Obtaining a most recent target MS location estimate from either the 91 1 emergency center or 

the location center 142. 
(2c) Inputting by the MBS operator an acknowledgment of the target MS to be located, and 
transmitting this acknowledgment to the 91 1 emergency response center via the transceiver 
30 1512. 



Subsequently, when the MBS 148 is in the seek state 1712, the MBS may commence toward the 
target MS location estimate provided. Note that it is likely that the MBS is not initially in direct signal 
contact with the target MS. Accordingly, in the seek state 1712 the following steps may be, for example, 
35 performed: 

(3a) The location center 142 or the 91 1 emergency response center may inform the target MS, via 
the fixed location base station network, to lower its threshold for soft hand-off and at least 
periodically boost its location signal strength. Additionally, the target MS may be informed to 
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scan for the pilot channel of the MBS 148. (Note the actions here are not, actions performed by 
the MBS 148 in the "seek state"; however, these actions are given here for clarity and 
completeness.) 

(3b) Repeatedly, as sufficient new MS location information is available, the location center 142 
5 provides new MS location estimates to the MBS 148 via the fixed location base station network. 

(3c) The MBS repeatedly provides the MBS operator with new target MS location estimates provided 

substantially by the location center via the fixed location base station network. 
(3d) The MBS 148 repeatedly attempts to detect a signal from the target MS using the PN code for 
the target MS. 

10 (3e) The MBS 148 repeatedly estimates its own location (as in other states as well), and receives 

MBS location estimates from the location center. 

Assuming that the MBS 148 and target MS 140 detect one another (which typically occurs when 
the two units are within .25 to 3 miles of one another), the MBS enters a contact state 1 71 6 when the 
target MS 140 enters a soft hand-off state with the MBS. Accordingly, in the contact state 1716, the 
15 following steps are, for example, performed: 

(4a) The MBS 148 repeatedly estimates its own location. 

(4b) Repeatedly, the location center 142 provides new target MS 140 and MBS location estimates to 

the MBS 148 via the fixed location base infrastructure network. 
(4c) Since the MBS 148 is at least in soft hand-off with the target MS 140, the MBS can estimate the 
20 direction and distance of the target MS itself using, for example, detected target MS signal 

strength and TOA as well as using any recent location center target MS location estimates. 
(4d) The MBS 148 repeatedly provides the MBS operator with new target MS location estimates 

provided using MS location estimates provided by the MBS itself and by the location center via 

the fixed location base station network. 

25 When the target MS 140 detects that the MBS pilot channel is sufficiently strong, the target MS 

may switch to using the MBS 148 as its primary base station. When this occurs, the MBS enters a control 
state 1 720, wherein the following steps are, for example, performed: 
(5a) The MBS 148 repeatedly estimates its own location. 

(5b) Repeatedly, the location center 142 provides new target MS and MBS location estimates to the 
30 MBS 1 48 via the network of base stations 1 22 (1 52). 

(5c) The MBS 148 estimates the direction and distance of the target MS 140 itself using, for 
example, detected target MS signal strength and TOA as well as using any recent location 
center target MS location estimates. 
(5d) The MBS 148 repeatedly provides the MBS operator with new target MS location estimates 
35 provided using MS location estimates provided by the MBS itself and by the location center 142 

via the fixed location base station network. 
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(5e) The MBS 148 becomes the primary base station for the target MS 140 and therefore controls at 

least the signal strength output by the target MS. 
Note, there can be more than one MBS 148 tracking or locating an MS 140. There can also be more 
than one target MS 140 to be tracked concurrently and each target MS being tracked may be stationary 
5 or moving. 

MBS Subsystem Architecture 

An MBS 148 uses MS signal characteristic data for locating the MS 140. The MBS 148 may use 
such signal characteristic data to facilitate determining whether a given signal from the MS is a "direct 
shot" or an multipath signal. That is, in one embodiment, the MBS 148 attempts to determine or detect 

10 whether an MS signal transmission is received directly, or whether the transmission has been reflected or 
deflected. For example, the MBS may determine whether the expected signal strength, and TOA agree in 
distance estimates for the MS signal transmissions. Note, other signal characteristics may also be used, 
if there are sufficient electronics and processing available to the MBS 148; i.e., determining signal phase 
and/or polarity as other indications of receiving a "direct shof from an MS 140. 

1 5 In one embodiment, the MBS 148 (Fig. 1 1 ) includes an MBS controller 1 533 for controlling the 

location capabilities of the MBS 148. In particular, APPENDIX A hereinbelow includes a description of the 
program "mobile.base.station.controller" referenced in Fig.1 1(1). Note, that the description of this 
program describes some of the functionality of the location controller 1535 for activating, e.g., the 
baseline location estimators 1540 (described hereinbelow). In particular, the MBS controller 1533 

20 initiates and controls the MBS state changes as described in Fig. 12. Additionally, the MBS controller 
1533 also communicates with the location controller 1535, wherein this latter controller controls MBS 
activities related to MBS location and target MS location. The location controller 1535 receives data 
input from an event generator 1537 that generates event records to be provided to the location controller 
1535. For example, records may be generated from data input received from: (a) the vehicle movement 

25 detector 1539 indicating that the MBS 148 has moved at least a predetermined amount and/or has 

changed direction by at least a predetermined angle, or (b) the MBS signal processing subsystem 1541 
indicating that the additional signal measurement data has been received from either the location center 
142 or the target MS 140. Note that the MBS signal processing subsystem 1541, in one embodiment, is 
similar to the signal processing subsystem 1 220 of the location center 1 42. Moreover, also note that 

30 there may be multiple command schedulers. In particular, a scheduler 1 528 for commands related to 
communicating with the location center 142, a scheduler 1530 for commands related to GPS 
communication (via GPS receiver 1531), a scheduler 1529 for commands related to the frequency and 
granularity of the reporting of MBS changes in direction and/or position via the MBS deadreckoning 
subsystem 1527 (note that this scheduler is potentially optional and that such commands may be 

35 provided directly to the deadreckoning estimator 1 544), and a scheduler 1 532 for communicating with the 
target MS(s) 140 being located. Further, it is assumed that there is sufficient hardware and/or software to 
perform commands in different schedulers substantially concurrently. 
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In order to display an MBS computed location of a target MS 140, a location of the MBS must be 
known or determined. Accordingly, each MBS 148 has a plurality of MBS location estimators (or 
hereinafter also simply referred to as location estimators) for determining the location of the MBS. Each 
such location estimator computes MBS location information such as MBS location estimates, changes to 
5 MBS location estimates, or, an MBS location estimator may be an interface for buffering and/or translating 
a previously computed MBS location estimate into an appropriate format. In particular, the MBS location 
module 1536, which determines the location of the MBS, may include the following MBS location 
estimators 1540 (also denoted baseline location estimators): 

(a) a GPS location estimator 1540a (not individually shown) for computing an MBS location 
io estimate using GPS signals, 

(b) a location center location estimator 1540b (not individually shown) for buffering and/or 
translating an MBS estimate received from the location center 142, 

(c) an MBS operator location estimator 1540c (not individually shown) for buffering and/or 
translating manual MBS location entries received from an MBS location operator, and 

15 (d) in some MBS embodiments, an LBS location estimator 1540d (not individually shown) for the 

activating and deactivating of LBS's 152. Note that, in high multipath areas and/or stationary 
base station marginal coverage areas, such low cost location base stations 152 (LBS) may be 
provided whose locations are fixed and accurately predetermined and whose signals are 
substantially only receivable within a relatively small range (e.g., 2000 feet), the range 
20 potentially being variable. Thus, by communicating with the LBS's 152 directly, the MBS 148 

may be able to quickly use the location information relating to the location base stations for 
determining its location by using signal characteristics obtained from the LBSs 152. 
Note that each of the MBS baseline location estimators 1540, such as those above, provide an actual 
MBS location rather than, for example, a change in an MBS location. Further note that it is an aspect of a 
25 wireless location system according to the novel method and/or system disclosed herein that additional MBS 
baseline location estimators 1540 may be easily integrated into the MBS location subsystem 1508 as 
such baseline location estimators become available. For example, a baseline location estimator that 
receives MBS location estimates from reflective codes provided, for example, on streets or street signs 
can be straightforwardly incorporated into the MBS location subsystem 1508. 
30 Additionally, note that a plurality of MBS location technologies and their corresponding MBS location 

estimators are utilized due to the fact that there is currently no single location technology available that is 
both sufficiently fast, accurate and accessible in substantially all terrains to meet the location needs of an 
MBS 148. For example, in many terrains GPS technologies may be sufficiently accurate; however, GPS 
technologies: (a) may require a relatively long time to provide an initial location estimate (e.g., greater 
35 than 2 minutes); (b) when GPS communication is disturbed, it may require an equally long time to provide 
a new location estimate; (c) clouds, buildings and/or mountains can prevent location estimates from being 
obtained; (d) in some cases signal reflections can substantially skew a location estimate. As another 
example, an MBS 148 may be able to use triangulation or ^lateralization technologies to obtain a location 
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estimate; however, this assumes that there is sufficient (fixed location) infrastructure BS coverage in the 
area the MBS is located. Further, it is well known that the multipath phenomenon can substantially distort 
such location estimates. Thus, for an MBS 148 to be highly effective in varied terrains, an MBS is 
provided with a plurality of location technologies, each supplying an MBS location estimate. 
5 In fact, much of the architecture of the location engine 139 could be incorporated into an MBS 148. 

For example, in some embodiments of the MBS 148, the following FOMs 1224 may have similar location 
models incorporated into the MBS: 

(a) a variation of the distance FOM 1224 wherein TOA signals from communicating fixed 
location BS's are received (via the MBS transceiver 1512) by the MBS and used for 

l o providing a location estimate; 

(b) a variation of the artificial neural net based FOMs 1224 (or more generally a location 
learning or a classification model) may be used to provide MBS location estimates via, for 
example, learned associations between fixed location BS signal characteristics and 
geographic locations; 

1 5 (c) an LBS location FOM 1 224 for providing an MBS with the ability to activate and deactivate 

LBS's to provide (positive) MBS location estimates as well as negative MBS location 
regions (i.e., regions where the MBS is unlikely to be since one or more LBS's are not 
detected by the MBS transceiver); 
(d) one or more MBS location reasoning agents and/or a location estimate heuristic agents for 
20 resolving MBS location estimate conflicts and providing greater MBS location estimate 

accuracy. For example, modules similar to the analytical reasoner module 1416 and the 
historical location reasoner module 1424. 
However, for those MBS location models requiring communication with the base station 
infrastructure, an alternative embodiment is to rely on the location center 142 to perform the computations 
25 for at least some of these MBS FOM models. That is, since each of the MBS location models mentioned 
immediately above require communication with the network of fixed location BS's 122 (152), it may be 
advantageous to transmit MBS location estimating data to the location center 142 as if the MBS were 
another MS 140 for the location center to locate, and thereby rely on the location estimation capabilities at 
the location center rather than duplicate such models in the MBS 148. The advantages of this approach 
30 are that: 

(a) an MBS is likely to be able to use less expensive processing power and software than that 
of the location center; 

(b) an MBS is likely to require substantially less memory, particularly for data bases, than that 
of the location center. 

35 As will be discussed further below, in one embodiment of the MBS 148, there are confidence values 

assigned to the locations output by the various location estimators 1540. Thus, the confidence for a 
manual entry of location data by an MBS operator may be rated the highest and followed by the 
confidence for (any) GPS location data, followed by the confidence for (any) location center location 142 
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estimates, followed by the confidence for (any) location estimates using signal characteristic data from 
LBSs. However, such prioritization may vary depending on, for instance, the radio coverage area 120. In 
an one embodiment of a wireless location system according to the novel method and/or system disclosed 
herein, it is an aspect that for MBS location data received from the GPS and location center, their 
5 confidences may vary according to the area in which the MBS 148 resides. That is, if it is known that for 
a given area, there is a reasonable probability that a GPS signal may suffer multipath distortions and that 
the location center has in the past provided reliable location estimates, then the confidences for these two 
location sources may be reversed. 

In one embodiment of a location system according to the novel method and/or system disclosed 
10 herein, MBS operators may be requested to occasionally manually enter the location of the MBS 148 
when the MBS is stationary for determining and/or calibrating the accuracy of various MBS location 
estimators. 

There is an additional important source of location information for the MBS 148 that is 
incorporated into an MBS vehicle (such as a police vehicle) that has no comparable functionality in the 

15 network of fixed location BS's. That is, the MBS 148 may use deadreckoning information provided by a 
deadreckoning MBS location estimator 1544 whereby the MBS may obtain MBS deadreckoning location 
change estimates. Accordingly, the deadreckoning MBS location estimator 1544 may use, for example, 
an on-board gyroscope 1550, a wheel rotation measurement device (e.g., odometer) 1554, and optionally 
an accelerometer (not shown). Thus, such a deadreckoning MBS location estimator 1544 periodically 

20 provides at least MBS distance and directional data related to MBS movements from a most recent MBS 
location estimate. More precisely, in the absence of any other new MBS location information, the 
deadreckoning MBS location estimator 1544 outputs a series of measurements, wherein each such 
measurement is an estimated change (or delta) in the position of the MBS 148 between a request input 
timestamp and a closest time prior to the timestamp, wherein a previous deadreckoning terminated. 

25 Thus, each deadreckoning location change estimate includes the following fields: 

(a) an "earliest timestamp" field for designating the start time when the deadreckoning location 
change estimate commences measuring a change in the location of the MBS; 

(b) a "latest timestamp" field for designating the end time when the deadreckoning location 
change estimate stops measuring a change in the location of the MBS; and 

30 (c) an MBS location change vector. 

That is, the "latest timestamp" is the timestamp input with a request for deadreckoning location data, and 
the "earliest timestamp" is the timestamp of the closest time, T, prior to the latest timestamp, wherein a 
previous deadreckoning output has its a timestamp at a time equal to T. 

Further, the frequency of such measurements provided by the deadreckoning subsystem 1527 may 

35 be adaptively provided depending on the velocity of the MBS 148 and/or the elapsed time since the most 
recent MBS location update. Accordingly, the architecture of at least some embodiments of the MBS 
location subsystem 1508 must be such that it can utilize such deadreckoning information for estimating 
the location of the MBS 148. 



In one embodiment of the MBS location subsystem 1508 described in further detail hereinbelow, the 
outputs from the deadreckoning MBS location estimator 1544 are used to synchronize MBS location 
estimates from different MBS baseline location estimators. That is, since such a deadreckoning output 
may be requested for substantially any time from the deadreckoning MBS location estimator, such an 
5 output can be requested for substantially the same point in time as the occurrence of the signals from 
which a new MBS baseline location estimate is derived. Accordingly, such a deadreckoning output can 
be used to update other MBS location estimates not using the new MBS baseline location estimate. 

It is assumed that the error with dead reckoning increases with deadreckoning distance. Accordingly, 
it is an aspect of the embodiment of the MBS location subsystem 1 508 that when incrementally updating 

10 the location of the MBS 148 using deadreckoning and applying deadreckoning location change estimates 
to a "most likely area" in which the MBS 148 is believed to be, this area is incrementally enlarged as well 
as shifted. The enlargement of the area is used to account for the inaccuracy in the deadreckoning 
capability. Note, however, that the deadreckoning MBS location estimator is periodically reset so that the 
error accumulation in its outputs can be decreased. In particular, such resetting occurs when there is a 

15 high probability that the location of the MBS is known. For example, the deadreckoning MBS location 
estimator may be reset when an MBS operator manually enters an MBS location or verifies an MBS 
location, or a computed MBS location has sufficiently high confidence. 

Thus, due to the MBS 148 having less accurate location information (both about itself and a target 
MS 140), and further that deadreckoning information must be utilized in maintaining MBS location 

20 estimates, a first embodiment of the MBS location subsystem architecture is somewhat different from the 
location engine 1 39 architecture. That is, the architecture of this first embodiment is simpler than that of 
the architecture of the location engine 1 39. However, it important to note that, at a high level, the 
architecture of the location engine 139 may also be applied for providing a second embodiment of the 
MBS location subsystem 1508, as one skilled in the art will appreciate after reflecting on the architectures 

25 and processing provided at an MBS 148. For example, an MBS location subsystem 1508 architecture 
may be provided that has one or more first order models 1224 whose output is supplied to, for example, a 
blackboard or expert system for resolving MBS location estimate conflicts, such an architecture being 
analogous to one embodiment of the location engine 139 architecture. 

Furthermore, it is also an important aspect of a wireless location system according to the novel method 

30 and/or system disclosed herein that, at a high level, the MBS location subsystem architecture may also 
be applied as an alternative architecture for the location engine 1 39. For example, in one embodiment of 
the location engine 139, each of the first order models 1224 may provide its MS location hypothesis 
outputs to a corresponding location track," analogous to the MBS location tracks described hereinbelow, 
and subsequently, a most likely MS current location estimate may be developed in a "current location 

35 track" (also described hereinbelow) using the most recent location estimates in other location tracks. 

Further, note that the ideas and methods discussed here relating to MBS location estimators 1540 
and MBS location tracks, and, the related programs hereinbelow are sufficiently general so that these 
ideas and methods may be applied in a number of contexts related to determining the location of a device 
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capable of movement and wherein the location of the device must be maintained in real time. For 
example, the present ideas and methods may be used by a robot in a very cluttered environment (e.g., a 
warehouse), wherein the robot has access: (a) to a plurality of "robot location estimators" that may 
provide the robot with sporadic location information, and (b) to a deadreckoning location estimator. 
5 Each MBS 148, additionally, has a location display (denoted the MBS operator visual user interface 

1 558 in Fig. 1 1 ) where area maps that may be displayed together with location data. In particular, MS 
location data may be displayed on this display as a nested collection of areas, each smaller nested area 
being the most likely area within (any) encompassing area for locating a target MS 140. Note that the 
MBS controller algorithm below may be adapted to receive location center 142 data for displaying the 

l o locations of other MBSs 1 48 as well as target MSs 1 40. 

Further, the MBS 148 may constrain any location estimates to streets on a street map using the MBS 
location snap to street module 1 562. For example, an estimated MBS location not on a street may be 
"snapped to" a nearest street location. Note that a nearest street location determiner may use "normal" 
orientations of vehicles on streets as a constraint on the nearest street location. Particularly, if an MBS 

15 148 is moving at typical rates of speed and acceleration, and without abrupt changes direction. For 
example, if the deadreckoning MBS location estimator 1544 indicates that the MBS 148 is moving in a 
northerly direction, then the street snapped to should be a north-south running street. Moreover, the MBS 
location snap to street module 1562 may also be used to enhance target MS location estimates when, for 
example, it is known or suspected that the target MS 140 is in a vehicle and the vehicle is moving at 

20 typical rates of speed. Furthermore, the snap to street location module 1 562 may also be used in 
enhancing the location of a target MS 140 by either the MBS 148 or by the location engine 139. In 
particular, the location estimator 1 344 or an additional module between the location estimator 1 344 and 
the output gateway 1 356 may utilize an embodiment of the snap to street location module 1 562 to 
enhance the accuracy of target MS 140 location estimates that are known to be in vehicles. Note that this 

25 may be especially useful in locating stolen vehicles that have embedded wireless location transceivers 
(MSs 140), wherein appropriate wireless signal measurements can be provided to the location center 
142. 

MBS Data Structure Remarks 

Assuming the existence of at least some of the location estimators 1 540 that were mentioned above, 
30 the discussion here refers substantially to the data structures and their organization as illustrated in Fig. 
13. 

The location estimates (or hypotheses) for an MBS 148 determining its own location each have an 
error or range estimate associated with the MBS location estimate. That is, each such MBS location 
estimate includes a "most likely MBS point location" within a "most likely area". The "most likely MBS 
35 point location" is assumed herein to be the centroid of the "most likely area." In one embodiment of the 
MBS location subsystem 1508, a nested series of "most likely areas" may be provided about a most likely 
MBS point location. However, to simplify the discussion herein each MBS location estimate is assumed 
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to have a single "most likely area". One skilled in the art will understand how to provide such nested 
"most likely areas" from the description herein. Additionally, it is assumed that such "most likely areas" are 
not grossly oblong; i.e., area cross sectioning lines through the centroid of the area do not have large 
differences in their lengths. For example, for any such "most likely area", A, no two such cross sectioning 
5 lines of A may have lengths that vary by more than a factor of two. 

Each MBS location estimate also has a confidence associated therewith providing a measurement of 
the perceived accuracy of the MBS being in the "most likely area" of the location estimate. 

A (MBS) location track" is an data structure (or object) having a queue of a predetermined length for 
maintaining a temporal (timestamp) ordering of "location track entries" such as the location track entries 
10 1770a, 1770b, 1774a, 1774b, 1778a, 1778b, 1782a, 1782b, and 1786a (Fig. 13), wherein each such MBS 
location track entry is an estimate of the location of the MBS at a particular corresponding time. 

There is an MBS location track for storing MBS location entries obtained from MBS location 
estimation information from each of the MBS baseline location estimators described above (i.e., a GPS 
location track 1750 for storing MBS location estimations obtained from the GPS location estimator 1540, a 
15 location center location track 1754 for storing MBS location estimations obtained from the location 
estimator 1540 deriving its MBS location estimates from the location center 142, an LBS location track 
1758 for storing MBS location estimations obtained from the location estimator 1540 deriving its MBS 
location estimates from base stations 122 and/or 152, and a manual location track 1762 for MBS operator 
entered MBS locations). Additionally, there is one further location track, denoted the "current location 
20 track" 1 766 whose location track entries may be derived from the entries in the other location tracks 

(described further hereinbelow). Further, for each location track, there is a location track head that is the 
head of the queue for the location track. The location track head is the most recent (and presumably the 
most accurate) MBS location estimate residing in the location track. Thus, for the GPS location track 
1750 has location track head 1770; the location center location track 1754 has location track head 1774; 
25 the LBS location track 1758 has location track head 1778; the manual location track 1762 has location 
track head 1782; and the current location track 1766 has location track head 1786. Additionally, for 
notational convenience, for each location track, the time series of previous MBS location estimations (i.e., 
location track entries) in the location track will herein be denoted the "path for the location track." Such 
paths are typically the length of the location track queue containing the path. Note that the length of each 
30 such queue may be determined using at least the following considerations: 

(i) In certain circumstances (described hereinbelow), the location track entries are removed 
from the head of the location track queues so that location adjustments may be made. In 
such a case, it may be advantageous for the length of such queues to be greater than the 
number of entries that are expected to be removed; 
35 (ii) In determining an MBS location estimate, it may be desirable in some embodiments to 

provide new location estimates based on paths associated with previous MBS location 
estimates provided in the corresponding location track queue. 
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Also note that it is within the scope of a wireless location system according to the novel method and/or 
system disclosed herein that the location track queue lengths may be a length of one. 
Regarding location track entries, each location track entry includes: 

(a) a "derived location estimate" for the MBS that is derived using at least one of: 

(i) at least a most recent previous output from an MBS baseline location estimator 1540 
(i.e., the output being an MBS location estimate); 

(ii) deadreckoning output information from the deadreckoning subsystem 1527. 
Further note that each output from an MBS location estimator has a "type" field that is used 
for identifying the MBS location estimator of the output. 

(b) an "earliest timestamp" providing the time/date when the earliest MBS location information 
upon which the derived location estimate for the MBS depends. Note this will typically be 
the timestamp of the earliest MBS location estimate (from an MBS baseline location 
estimator) that supplied MBS location information used in deriving the derived location 
estimate for the MBS 148. 

(c) a "latest timestamp" providing the time/date when the latest MBS location information upon 
which the derived location estimate for the MBS depends. Note that earliest timestamp = 
latest timestamp only for so called "baseline entries" as defined hereinbelow. Further note 
that this attribute is the one used for maintaining the "temporal (timestamp) ordering" of 
location track entries. 

(d) A "deadreckoning distance" indicating the total distance (e.g., wheel turns or odometer 
difference) since the most recently previous baseline entry for the corresponding MBS 
location estimator for the location track to which the location track entry is assigned. 

For each MBS location track, there are two categories of MBS location track entries that may be 
inserted into a MBS location track: 

(a) "baseline" entries, wherein each such baseline entry includes (depending on the location 
track) a location estimate for the MBS 148 derived from: (i) a most recent previous output 
either from a corresponding MBS baseline location estimator, or (ii) from the baseline 
entries of other location tracks (this latter case being the for the "current" location track); 

(b) "extrapolation" entries, wherein each such entry includes an MBS location estimate that 
has been extrapolated from the (most recent) location track head for the location track (i.e., 
based on the track head whose "latest timestamp" immediately precedes the latest 
timestamp of the extrapolation entry). Each such extrapolation entry is computed by using 
data from a related deadreckoning location change estimate output from the 
deadreckoning MBS location estimator 1544. Each such deadreckoning location change 
estimate includes measurements related to changes or deltas in the location of the MBS 
148. More precisely, for each location track, each extrapolation entry is determined using: 
(i) a baseline entry, and (ii) a set of one or more (i.e., all later occurring) deadreckoning 
location change estimates in increasing "latest timestamp" order. Note that for notational 
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convenience this set of one or more deadreckoning location change estimates will be 
denoted the "deadreckoning location change estimate set" associated with the 
extrapolation entry resulting from this set. 
(c) Note that for each location track head, it is either a baseline entry or an extrapolation entry. 
5 Further, for each extrapolation entry, there is a most recent baseline entry, B, that is earlier 

than the extrapolation entry and it is this B from which the extrapolation entry was 
extrapolated. This earlier baseline entry, B, is hereinafter denoted the "baseline entry 
associated with the extrapolation entry." More generally, for each location track entry, T, 
there is a most recent previous baseline entry, B, associated with T, wherein if T is an 
10 extrapolation entry, then B is as defined above, else if T is a baseline entry itself, then T=B. 

Accordingly, note that for each extrapolation entry that is the head of a location track, there 
is a most recent baseline entry associated with the extrapolation entry. 
Further, there are two categories of location tracks: 

(a) "baseline location tracks," each having baseline entries exclusively from a single 
15 predetermined MBS baseline location estimator; and 

(b) a "current" MBS location track having entries that are computed or determined as "most 
likely" MBS location estimates from entries in the other MBS location tracks. 

MBS Location Estimating Strategy 

In order to be able to properly compare the track heads to determine the most likely MBS location 

20 estimate it is an aspect of a wireless location system according to the novel method and/or system 
disclosed herein that the track heads of all location tracks include MBS location estimates that are for 
substantially the same (latest) timestamp. However, the MBS location information from each MBS 
baseline location estimator is inherently substantially unpredictable and unsynchronized. In fact, the only 
MBS location information that may be considered predicable and controllable is the deadreckoning 

25 location change estimates from the deadreckoning MBS location estimator 1 544 in that these estimates 
may reliably be obtained whenever there is a query from the location controller 1 535 for the most recent 
estimate in the change of the location for the MBS 148. Consequently (referring to Fig. 1 3), 
synchronization records 1790 (having at least a 1790b portion, and in some cases also having a 1790a 
portion) may be provided for updating each location track with a new MBS location estimate as a new 

30 track head. In particular, each synchronization record includes a deadreckoning location change estimate 
to be used in updating all but at most one of the location track heads with a new MBS location estimate by 
using a deadreckoning location change estimate in conjunction with each MBS location estimate from an 
MBS baseline location estimator, the location track heads may be synchronized according to timestamp. 
More precisely, for each MBS location estimate, E, from an MBS baseline location estimator, the MBS 

35 also substantially simultaneously queries the deadreckoning MBS location estimator for a corresponding 
most recent change in the location of the MBS 148. Accordingly, E and the retrieved MBS 
deadreckoning location change estimate, C, have substantially the same latest timestamp". Thus, the 
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location estimate E may be used to create a new baseline track head for the location track having the 
corresponding type for E, and C may be used to create a corresponding extrapolation entry as the head 
of each of the other location tracks. Accordingly, since for each MBS location estimate, E, there is a MBS 
deadreckoning location change estimate, C, having substantially the same "latest timestamp", E and C 

5 will be hereinafter referred as "paired." 

The foregoing description of preferred embodiments of the present invention has been provided for 
the purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to 
the precise forms disclosed herein. Modifications and variations commensurate with the description 
herein will be apparent those skilled in the art and are intended to be within the scope of the present 

10 invention to the extent permitted by the relevant art. The embodiments provided are for enabling others 
skilled in the art to understand the invention, its various embodiments and modifications as are suited for 
uses contemplated. It is intended that the scope of the invention be defined by the following claims and 
their equivalents. 
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APPENDIX A: MBS Function Embodiments 

Mobile Base Station Controller Program 
mobile_base_station_controller() 

{ 

waitJorJnputjofJirstJMBSJocation(ever\\y t r "event" is a record (object) with MBS location data */ 
WHILE (no MBS operator input to exit) DO 

CASE OF (event): /* determine the type of "event" and process it. 7 

MBS LOCATION DATA RECEIVED FROM GPS: 

MBS LOCATION DATA RECEIVED FROM LBS: 

MBS LOCATION DATA RECEIVED FROM ANY OTHER HIGHLY RELIABLE MBS LOCATION 
SOURCES (EXCEPT LOCATION CENTER): 

{ 

MBS_new__est <— getjiew_MBSJocation_using_estimate{eyen{)\ 
r Note, whenever a new MBS location estimate is entered as a baseline estimate into one 
of the location tracks, the other location tracks must be immediately updated with any 
deadreckoning location change estimates so that all location tracks are substantially 
updated at the same time. 7 
deadreck_est <— getjjeadreckoningJocationjzhangejastimate(even{y, 
MBS_curr_est <- DETERMINE JWBS_LOCATION_ESTIMATE(MBS_new_est, 
deadreck_est); 

if (MBS_curr_est.confidence > a predetermined high confidence threshold) then 
resetjdeadreckoning_MBSJocationj3stimatoiieven\)] 

r deadreckoning starts over from here. 7 

I* Send MBS location information to the Location Center. 7 

if ( MBS has not moved since the last MBS location estimate of this type and is not now 
moving) then 

{ 

configure the MBS on-board transceiver (e.g., MBS-MS) to immediately transmit location 
signals to the fixed location BS network as if the MBS were an ordinary location 
device (MS); 

communicate with the Location Center via the fixed location BS infrastructure the 
following: 

(a) a locate me" signal, 
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(b) MBS_curr_est, 

(c) MBS_new_est and 

(d) the timestamp for the present event. 

Additionally, any location signal information between the MBS and the present target 
MS may be transmitted to the Location Center so that this information may also be 
used by the Location Center to provide better estimates of where the. MBS is. 
Further, if the MBS determines that it is immediately adjacent to the target MS and 
also that its own location estimate is highly reliable (e.g., a GPS estimate), then the 
MBS may also communicate this information to the Location Center so that the 
Location Center can: (a) associate any target MS location signature cluster data 
with the fixed base station infrastructure with the location provided by the MBS, and 
(b) insert this associated data into the location signature data base of the Location 
Center as a verified cluster of "random loc sigs"; 
r note, this transmission preferably continues (i.e., repeats) for at least a 
predetermined length of time of sufficient length for the Signal Processing 
Subsystem to collect a sufficient signal characteristic sample size. */ 

} 

else SCHEDULE an event (if none scheduled) to transmit to the Location Center the 
following: (a) MBS_curr_est, and (b) the GPS location of the MBS and the time of the 
GPS location estimate; 

r Now update MBS display with new MBS location; note, MBS operator must 
request MBS locations on the MBS display; if not requested, then the following 
call does not do an update. 7 
updateJ^BSjoperatorjdisplayj^ithJ^BSjast(MBSj^unj8s\)] 

> 

SINCE LAST MBS LOCATION UPDATE 

MBS HAS MOVED A THRESHOLD DISTANCE; { 

deadreck_est <— getjdeadre<XoningJocationjzhangejestimate{even\)] 

I* Obtain from MBS Dead Reckoning Location Estimator a new dead reckoning 
MBS location estimate having an estimate as to the MBS location change 
from the location of the last MBS location provided to the MBS. 7 
MBS_curr_est <— DETERMINE_MBS_LOCATIONJESTIMATE(NULL, deadreck.est); 

r this new MBS estimate will be used in new target MS estimatesV 
update JABS jdisplayjNithjjpdated JABS Jocation(M^ 
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SCHEDULE an event (if none scheduled) to request new GPS location data for MBS; 
SCHEDULE an event (if none scheduled) to request communication with Location Center 

(LC) related to new MBS location data; 
SCHEDULE an event (if none scheduled) to request new LBS location communication 
between the MBS and any LBS's that can detect the MBS; 

r Note, in some embodiments the processing of MBS location data from LBS's may 
be performed automatically by the Location Center, wherein the Location Center 
uses signal characteristic data from the LBS's in determining an estimated 
location of the MBS. 7 
SCHEDULE an event (if none scheduled) to obtain new target MS signal characteristics 
from MS; I* i.e., may get a better target MS location estimate now. */ 

} 

TIMER HAS EXPIRED SINCE LAST RELIABLE TARGET MS LOCATION INFORMATION 
OBTAINED: { 

SCHEDULE an event (if none scheduled) to request location communication with the target MS, 
the event is at a very high priority; 
RESET timer for target MS location communication; r Try to get target MS location 
communication again within a predetermined time. Note, timer may dynamically 
determined according to the perceived velocity of the target MS. */ 

} 

LOCATION COMMUNICATION FROM TARGET MS RECEIVED: { 

MS_raw_signal_data <— getJdSjsignaljcharaderistic_raw_data(eveniy t 

I* Note, a MS_raw_signal_data" is an object having substantially the unfiltered 
signal characteristic values for communications between the MBS and the 
target MS as well as timestamp information. 7 
Construct a message for sending to the Location Center, wherein the message includes at 
least u MS_raw_signal_data" and "MBS_curr_esf so that the Location Center can also 
compute an estimated location for the target MS; 
SCHEDULE an event (if none scheduled) to request communication with Location Center 
(LC) for sending the constructed message; 

r Note, this data does not overwrite any previous data waiting to be sent to the LC. 
7 

MS_signal_data <— get_MSjsignaljcharacteristic_data(ever\\)\ 

/* Note, the MS signal data obtained above is, in one embodiment, "raw" signal 
data. However, in a second embodiment, this data is filtered substantially as in 
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the Location Center by the Signal Processing Subsystem. For simplicity of 
discussion here, it is assumed that each MBS includes at least a scaled down 
version of the Signal Processing Subsystem (see FIG. 11). 7 

MS_new_est <— DETERMINE_MS_MOST__RECENT_ESTIMATE{MBS_cun_e$t, 
MS_curr_est, MS_signal_data); 
r May use forward and reverse TOA, TDOA, signal power, signal strength, 
and signal quality indicators. Note, u MS_curr_est M includes a timestamp of 
when the target MS signals were received. */ 

if (MSjiew__est.confidence > min_MS_confidence ) then 

{ 

ma^_MS_esLasJempora/y(MS_new_est); 

/* Note, it is assumed that this MS location estimate is "temporary" in the sense 
that it will be replaced by a corresponding MS location estimate received from 
the Location Center that is based on the same target MS raw signal data. That 
is, if the Location Center responds with a corresponding target MS location 
estimate, E, while "MS__new__est B is a value in a "moving window" of target MS 
location estimates (as described hereinbelow), then E will replace the value of 
"MS_new_est". Note, the moving window may dynamically vary in size 
according to, for example, a perceived velocity of the target MS and/or the 
MBS. 7 

MS_moving_window <— getJrfSjnoving_window{eveni)] 

r get moving window of location estimates for this target MS. */ 

addJMSjestimateJoJrfSJocationjwindo^MSj]^ MS_moving_window); 
r Since any given single collection of measurements related to locating the 
target MS may be potentially misleading, a "moving window" of location 
estimates are used to form a "composite location estimate" of the target MS. 
This composite location estimate is based on some number of the most recent 
location estimates determined. Such a composite location estimate may be, 
for example, analogous to a moving average or some other weighting of target 
MS location estimates. Thus, for example, for each location estimate (i.e., at 
least one MS location area, a most likely single location, and, a confidence 
estimate) a centroid type calculation may be performed to provide the 
composite location estimate.7 

MS_curr_est <— DETERMINE_MS_LOCATION_ESTIMA TE(MS_moving_window); 
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r DETERMINE new target MS location estimate. Note this may an average 
location or a weighted average location. */ 
mmovejscheduled_events( tt T ARGETJAS^SCHEDULE", event.MSJD); 

r REMOVE ANY OTHER EVENTS SCHEDULED FOR REQUESTING 
LOCATION COMMUNICATION FROM TARGET MS 7 

} 

else r target MS location data received but it is not deemed to be reliable (e.g., too much 
multipath and/or inconsistent measurements, so SCHEDULE an event (if none 
scheduled) to request new location communication with the target MS, the event is 
at a high priority*/ 

addJo_scheduied_events( u TARGEY_WiS_SCHEDUlE\ event.MSJD); 
updateJ^BSj)peratorjdisplayj^ith_MSjest(MSj:un_es\)] 

r The MBS display may use various colors to represent nested location areas 
overlayed on an area map wherein, for example, 3 nested areas may be displayed 
on the map overlay: (a) a largest area having a relatively high probability that the 
target MS is in the area (e.g., >95%); (b) a smaller nested area having a lower 
probability that the target MS is in this area (e.g., >80%); and (c) a smallest area 
having the lowest probability that the target MS is in this area (e.g., >70%). Further, 
a relatively precise specific location is provided in the smallest area as the most 
likely single location of the target MS. Note that in one embodiment, the colors for 
each region may dynamically change to provide an indication as to how high their 
reliability is; e.g., no colored areas shown for reliabilities below, say, 40%; 40-50% is 
purple; 50-60% is blue; 60-70% is green; 70-80% is amber 80-90% is white; and red 
denotes the most likely single location of the target MS. Further note the three 
nested areas may collapse into one or two as the MBS gets closer to the target MS. 
Moreover, note that the collapsing of these different areas may provide operators in 
the MBS with additional visual reassurance that the location of the target MS is 
being determined with better accuracy.*/ 

r Now RESET timer for target MS location communication to try to get target MS 
location communication again within a predetermined time. */ 

reseU/me/t w TARGET_MS_SCHEDULE B , event.MSJD); 

} 

COMMUNICATION OF LOCATION DATA TO MBS FROM LOCATION CENTER: { 

r Note, target MS location data may be received from the Location Center in the seek state, 
contact state and the control state. Such data may be received in response to the MBS 
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sending target MS location signal data to the Location Center (as may be the case in 
the contact and control states), or such data may be received from the Location Center 
regardless of any previously received target MS location sent by the MBS (as may be 
the case in the seek, contact and control states). */ 
if ( (the timestamp of the latest MBS location data sent to the Location Center) <= 
(the timestamp returned by this Location Center communication identifying the MBS 
location data used by the Location Center for generating the MBS location data of the 
present event) ) 

then r use the LC location data since it is more recent than what is currently being used. 
7 

{ 

MBS_new_est <— get_Location_Center_MBS_est(eventy, 
deadreck_est <— getjJeadreckoningJocationj^hange_estimate{e\feni)] 
MBS_curr_est < — DETERMINEJVIBSJ-OCATION_ESTIMATE(MBSjiew_est, 
deadreck_est); 

if (MBS_curr_esLconfidence > a predetermined high confidence threshold) then 

resetjdeadreckoning^MBSJocationjBStimatorieyetA): 
updateJ^BSjoperatorjJisplayj^ith_MBSjest(MBSj^unjss\)) 

} 

if ( (the timestamp of the latest target MS location data sent to the Location Center) <= 
(the timestamp returned by this Location Center communication identifying the MS 
location data used by the Location Center for generating the target MS location 
estimate of the present event)) 
then r use the MS location estimate from the LC since it is more recent than what is 
currently being used. */ 

{ 

MS_new_est <— get Socat'onJCenter_MSjast(eyer\i)\ 

r This information includes error or reliability estimates that may be used in 
subsequent attempts to determine an MS location estimate when there is no 
communication with the LC and no exact (GPS) location can be obtained. That 
is, if the reliability of the target MS's location is deemed highly reliable, then 
subsequent less reliable location estimates should be used only to the degree 
that more highly reliable estimates become less relevant due to the MS moving 
to other locations. V 

MS_moving_window <— getJMSjnovingjwindowieventy, 

r get moving window of location estimates for this target MS. 7 
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if ( (the Location Center target MS estimate utilized the MS location signature data 
supplied by the MBS) then 
if (a corresponding target MS location estimate marked as "temporary" is still 
in the moving window) 

then /* It is assumed that this new target MS location data is still timely 
(note the target MS may be moving); so replace the temporary 
estimate with the Location Center estimate. */ 
replace the temporary target MS location estimate in the moving 
window with "MS_new_est"; 
else r there is no corresponding "temporary" target MS location in the 
moving window; so this MS estimate must be too old; so don't 
use it. */ 

else F the Location Center did not use the MS location data from the MBS even 
though the timestamp of the latest MS location data sent to the Location 
Center is older that the MS location data used by the Location Center to 
generate the present target MS location estimate. Use the new MS 
location data anyway. Note there isn't a corresponding "temporary" target 
MS location in the moving window. */ 
addJ^SjestimateJoJ^SJocationj^indo^MSj^e^esi): 

} 

else T the MS location estimate from the LC is not more recent than the latest MS location 
data sent to the LC from the MBS. */ 

if (a corresponding target MS location estimate marked as "temporary" is still in the 
moving window) 

then I* It is assumed that this new target MS location data is still timely (note 
the target MS may be moving); so replace the temporary estimate with 
the Location Center estimate. */ 
replace the temporary target MS location estimate in the moving window 
with "MS_new_est"; 
else r there is no corresponding "temporary" target MS location in the 

moving window; so this MS estimate must be too old; so don't use it. 
7 

MS_curr_est <— DETERMINE_MS_LOCATION_ESTIMA 7E(MS_moving_window); 

update_MBSjoperatorjdisplayj^ith_MS_est(MS_cun_es\y, 

reset _f/meA( tt LC_COMMUNICATION" ( event.MSJD); 
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NO COMMUNICATION FROM LC: { 

/* i.e., too long a time has elapsed since last communication from LC. 7 
SCHEDULE an event (if none scheduled) to request location data (MBS and/or target MS) 

from the Location Center, the event is at a high priority; 
resefJ/meffLC^COMMUNICATION", event.MSJD); 

} 

REQUEST TO NO LONGER CONTINUE LOCATING THE PRESENT TARGET MS: { 

if (event not from operator) then 

request MBS operator verification; 

else{ 

REMOVE the current target MS from the list of MSs currently being located and/or 
tracked; 

SCHEDULE an event (if none scheduled) to send communication to the Location Center 

that the current target MS is no longer being tracked; 
PURGE MBS of all data related to current target MS except any exact location data for 

the target MS that has not been sent to the Location Center for archival purposes; 

} 

} 

REQUEST FROM LOCATION CENTER TO ADD ANOTHER TARGET MS TO THE LIST OF 
MSs BEING TRACKED; { 

r assuming the Location Center sends MBS location data for a new target MS to locate 
and/or track (e.g., at least a new MS ID and an initial MS location estimate), add this 
new target MS to the list of MSs to track. Note the MBS will typically be or transitioning 
to in the seek state*/ 
if (event not from operator) then 

request MBS operator verification; 

else{ 

INITIALIZE MBS with data received from the Location Center related to the estimated 
location of the new target MS; /* e.g., initialize a new moving window for this new 
target MS; initialize MBS operator interface by graphically indicating where the new 
target MS is estimated to be. */ 

CONFIGURE MBS to respond to any signals received from the new target MS by 
requesting location data from the new target MS; 
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INITIALIZE timer for communication from LC; f A timer may be set per target MS on list. 
7 

} 

} 

REQUEST TO MANUALLY ENTER A LOCATION ESTIMATE FOR MBS (FROM AN MBS 
OPERATOR): { 

/* Note, MBS could be moving or stationary. If stationary, then the estimate for the location of 
the MBS is given high reliability and a small range (e.g., 20 feet). If the MBS is moving, 
then the estimate for the location of the MBS is given high reliability but a wider range 
that may be dependent on the speed of the MBS. In both cases, if the MBS operator 
indicates a low confidence in the estimate, then the range is widened, or the operator 
can manually enter a range*/ 
MBS_new_est <— getjiewJrfBSJocationj3StJmmjoperator{eventy t r The estimate may 

be obtained, for example, using a light pen on a displayed map 7 
if (operator supplies a confidence indication for the input MBS location estimate) then 

MBS_new_est.confidence <— gef_WBS_operafor_conffofence_of_esf/mafe(event); 
else MBS_new_est.confidence <— 1 ; I* This is the highest value for a confidence. 7 
deadreck_est <— getjdeadre(^oningJo<^tionj^angejdstimate{ewen\)] 
MBS_curr_est <- DETERMINE_MBS_LOCATION_ESTIMATE(MBS_new_est, 
deadreck_est ); 

if (MBS_curr_est.confidence > a predetermined high confidence threshold) then 

resetjdeadreckoning_MBSJocationjBstimatoi{eveni)\ 
update_MBSjnperatorjdi$playjwith„MBSjest{MBSj;u 

r Note, one reason an MBS operator might provide a manual MBS input is that the MBS 
might be too inaccurate in its location. Moreover, such inaccuracies in the MBS location 
estimates can cause the target MS to be estimated inaccurately, since target MS signal 
characteristic values may be utilized by the MBS to estimate the location of the target MS 
as an offset from where the MBS is. Thus, if there are target MS estimates in the moving 
window of target MS location estimates that are relatively close to the location 
represented by B MBS„curr.est", then these select few MS location estimates may be 
updated to reflect a more accurate MBS location estimate. 7 

MS_moving_window <— get_MS__moving_window(eventy, 

if (MBS has not moved much since the receipt of some previous target MS location that is 
still being used to location the target MS) 
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then 

{ 

UPDATE those target MS location estimates in the moving window according to the 

new MBS location estimate here; 
MS_curr_est <— DE7ERM/A/E_MS_LOC/17/0/V_£Sr/MAr£(MS_moving_window); 
updateJrfBSj)peratorjdi$playjMith_MS_est(M$jzuFr_^^ 

} 

} 

} I* end case statement */ 

Lower Level MBS Function Descriptions 

I* PROCEDURE: 

DETERMINE JABS J-OCATION_ESTIMATE REMARKS: 

It is assumed that with increasing continuous dead reckoning without additional MBS location 
verification, the potential error in the MBS location increases. 

It is assumed that each MBS location estimate includes: (a) a most likely area estimate surrounding a 
central location and (b) a confidence value of the MBS being in the location estimate. 

The confidence value for each MBS location estimate is a measurement of the likelihood of the MBS 
location estimate being correct. More precisely, a confidence value for a new MBS location estimate is a 
measurement that is adjusted according to the following criteria: 

(a) the confidence value increases with the perceived accuracy of the new MBS location estimate 
(independent of any current MBS location estimate used by the MBS), 

(b) the confidence value decreases as the location discrepancy with the current MBS location 
increases, 

(c) the confidence value for the current MBS location increases when the new location estimate is 
contained in the current location estimate, 

(d) the confidence value for the current MBS location decreases when the new location estimate is 
not contained in the current location estimate, and 

Therefore, the confidence value is an MBS location likelihood measurement which takes into account the 
history of previous MBS location estimates. 

It is assumed that with each MBS location estimate supplied by the Location Center there is a 
default confidence value supplied which the MBS may change. 
7 

DETERMINEJVIBS_LOCATION_ESTIMATE(MBS_new_est, deadreck_est) 
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I* Add the pair, u MBS_new_est B and B deadreck_esf to the location tracks and determine a new current 
MBS location estimate. 

Input: MBS_new_est A new MBS baseline location estimate to use in determining the location 

of the MBS, but not a (deadreckoning) location change estimate 
5 deadreck_est The deadreckoning location change estimate paired with . 

"MBS.new.esf. */ 

{ 

if (MBS_new_est is not NULL) then /* the w deadreck_esf is paired with B MBS_new_esf */ 
{ 

l o if (all MBS location tracks are empty) then 

{ 

insert u MBS_new_esf as the head of the location track of type, "MBS_new_est.type B ; 
insert "MBS_new_esf as the head of the current track; /* so now there is a 
"MBS^cu^est" MBS location estimate to use */ 
15 MBS_curr_est <— gef_cu/r_esf(MBS_new_est.MSJD); /* from current location track */ 

} 

else r there is at least one non-empty location track in addition to the current location track 
being non-empty*/ 

{ 

20 if (MBS.new_.est is of type MANUAL.ENTRY) then 

{ /* MBS operator entered an MBS location estimate for the MBS; so must use 
it*/ 

MBS_curr_est <— add_location_entry(MBS_new_est, deadreck_est); 

} 

25 else r "MBS_new_est" is not of type MANUAL.ENTRY */ 

if (the MBS location track of type, u MBS_new_est.type", is empty) then 
{ r some other location track is non-empty 7 

MBS_curr_est <— add_Jocation_entry(MBS_new_est, deadreck_est); 

} 

30 else /* "MBS_new_est.type" location track is non-empty and "MBS_new_est" 

is not of type MANUAL.ENTRY 7 
{ /* In the next statement determine if "MBS_new_est" is of at least 

minimal useful quality in comparison to any previous estimates of the 

same type; see program def n below 7 
35 continue_to_process_new_est <- FILTER(MBS_new_est); 
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if (continue_to_processjiew_est) then f "MBS_jiew_esf is of sufficient 

quality to continue processing. 7 

{ 

MBS_curr_est <— addJocation_entry(MBS_new_est, deadreck_est); 
5 }r end *MBS_new_j>st" not filtered out 7 

else r "MBSjiew__esf is filtered out; do nothing */; 
y* end else 7 
}/* end else at least one non-empty location track 7 

} 

l o else r MBS_new_est is NULL; thus only a deadreckoning output is to be added to location tracks 7 

{ 

extrapolation_entry <— c/Bafe_an_exfrapo/af/on_enf/y_from(deadreck_est); 

insert JntojeveryJocationJrack(exirBpo\aUon_er)\ry)' t /* including the "current location track" 7 

MBS_curr_est <— gef_cun-_esf(MBS_new_est.MSJD); /* from current location track 7 

15 } 

RETURN(MBS_curr_est); 
} END I* DETERMINE_MBS_LOCATION_ESTIMATE 7 

add_location_entry(MBS_new_est ) deadreck_est); 

20 r This function adds the baseline entry, "MBS_new_esf and its paired deadreckoning location change 
estimate, M deadreck_esr to the location tracks, including the "current location track". Note, however, that 
this function will roll back and rearrange location entries, if necessary, so that the entries are in latest 
timestamp order. 

Returns: MBS_curr_est 7 

25 { 

if (there is a time series of one or more dead reckoning extrapolation entries in the location track of 
type "MBS_new_est.type" wherein the extrapolation entries have a "latest timestamp" more recent 
than the timestamp of "MBS_new_esf ) then 
{ f Note, this condition may occur in a number of ways; e.g., (a) an MBS location estimate 
30 received from the Location Center could be delayed long enough (e.g., 1-4 sec) because of 

transmission and processing time; (b) the estimation records output from the MBS baseline 
location estimators are not guaranteed to be always presented to the location tracks in the 
temporal order they are created. 7 
roll back all (any) entries on all location tracks, including the "current" track, in "latest timestamp" 
35 descending order, until a baseline entry, B, is at the head of a location track wherein B is a 

most recent entry having a "latest timestamp" prior to "MBS_new_esf ; let "stack" be the 
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stack of a location track entries rolled off the location tracks, wherein an entry in the stack 
is either a baseline location entry and a paired deadreckoning location change estimate, 
or, an unpaired deadreckoning location change estimate associated with a NULL for the 
baseline location entry; 

insert "MBS_new_est" at the head of the location track of type B MBS_new_est.type n as a new 
baseline entry; 

insert the extrapolation entry derived from "deadreck_esf in each of the other baseline location 
tracks except the current track; 

/* It is important to note that "deadreck_est B includes the values for the change in the MBS 
location substantially for the time period between the timestamp, T, of u MS_new_esf and the 
timestamp of the closest deadreckoning output just before T. Further note that if there are any 
extrapolation entries that were rolled back above, then there is an extrapolation entry, E, 
previously in the location tracks and wherein E has an earliest timestamp equal to the latest 
timestamp of B above. Thus, all the previous extrapolation entries removed can be put back 
if E is modified as follows: the MBS location change vector of E (denoted herein as E.delta) 
becomes E.delta - pocation change vector of "deadreck.est"]. */ 

MBS_curr_est <— UPDATE_CURR_EST(MBS_new_est, deadreck.est); 

if (the extrapolation entry E exists) then /* i.e., "stack" is not empty */ 

{ 

modify the extrapolation entry E as per the comment above; 
I* now fix things up by putting all the rolled off location entries back, including the 
"current location track" */ 

do until "stack" is empty 
{ 

stackjop <— pop_sfac/c(stack); 

r "stackjop" is either a baseline location entry and a paired deadreckoning 
location change estimate, or, an unpaired deadreckoning location change 
estimate associated with a NULL for the baseline location entry */ 
MBS_nxt_est <— get_ba$eline_entry{$iackJop)] 
deadreck_est <— getjdeadreckoning_entry(s\ackJopy } 
MBS_curr_est <- DETERMINE_MBS_LOCATION_ESTIMATE(MBS_nxt_est, 
deadreck_est); 

} 

} 
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else F there is no deadreckoning extrapolation entries in the location track of type 

B MBS_new_est.type" wherein the extrapolation entries have a latest timestamp" more recent 
than the timestamp of "MBS_nevv_esf. So just insert "MBS_new_est n and u deadreck_esf .7 

{ 

5 insert u MBS_new_est B at the head of the location track of type w MBS_new_est.type" as a new 

baseline entry; 

insert the extrapolation entry derived from u deadreck_est B in each of the other location tracks 
except the current track; 

MBS_curr_est <— UPDATE_CURR_EST(MBS_new_est, deadreck_est); F see prog def n below 
10 */ 
} 

RETURN(MBS_curr_est); 
} F end add Jocation_entry */ 

1 5 FILTER(MBS_new_est) 

F This function determines whether u MBS_new_est" is of sufficient quality to insert into ifs corresponding 
MBS location track. It is assumed that the location track of "MBS_new_est.type" is non-empty. 
Input: MBS_new_est A new MBS location estimate to use in determining the location 

of the MBS. 

20 Returns: FALSE if "MBS_new_esf was processed here (i.e., filtered out), 

TRUE if processing with "MBS_new_est" may be continued . */ 

{ 

continue _to_process_new_est <-TRUE; F assume u MBS_new_esf will be good enough to use as an 

MBS location estimate */ 
25 F see if 8 MBS_new_est" can be filtered out. 7 

if (the confidence in MBS_new_est < a predetermined function of the confidence(s) of previous MBS 
location estimates of type a MBS_new_est.type n ) 

F e.g., the predetermined function here could be any of a number of functions that provide a 
minimum threshold on what constitutes an acceptable confidence value for continued 
30 processing of "MBS_new_esf. The following is an example of one such predetermined 

function: K*(confidence of u MBS_new_est.type" location track head) for some K, 
0<K<= 1 .0, wherein K varies with a relative frequency of estimates of type 
u MBS_new_est.type B not filtered; e.g.., for a given window of previous MBS location 
estimates of this type, K= (number of MBS location estimates of "MBSjiew_est.type n 
35 not filtered)/(the total number of estimates of this type in the window). Note, such filtering 

here may be important for known areas where, for example, GPS signals may be 
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potentially reflected from an object (i.e., multipath), or, the Location Center provides an 
MBS location estimate of very low confidence. For simplicity, the embodiment here 
discards any filtered location estimates. However, in an alternative embodiment, any 
such discarded location estimates may be stored separately so that, for example, if no 
5 additional better MBS location estimates are received, then the filtered or discarded 

location estimates may be reexamined for possible use in providing a better subsequent 
MBS location estimate.*/ 
then continue_to_process_new_est <- FALSE; 
else if (an area for H MBS_new_esf > a predetermined function of the corresponding area(s) of entries in 
1 o the location track of type u MBS_new__est.type") 

/* e.g., the predetermined function here could be any of a number of functions that provide 
a maximum threshold on what constitutes an acceptable area size for continued 
processing of "MBS_new_esf . The following are examples of such predetermined 
functions: (a) the identity function on the area of the head of the location track of type 
1 5 "MBS_new_est.type n ; or, (b) K*(the area of the head of the location track of type 

u MBS_new_est.type"), for some K, K>=1 .0, wherein for a given window of previous MBS 
location estimates of this type, K= (the total number of estimates in the window)/ (number 
of these location estimates not filtered); note, each extrapolation entry increases the area 
of the head; so areas of entries at the head of each location track type grow in area as 
20 extrapolation entries are applied. */ 

then continue_to_process_new_est <- FALSE; 
RETURN(continue_to_process_new_est) 
} 

25 UPDATEjCURRJEST(MBSjiewjest, deadreck_est) 

r This function updates the head of the "current 0 MBS location track whenever H MBS_new_esf is 
perceived as being a more accurate estimate of the location of the MBS. 
Input: MBS_new_est A new MBS location estimate to use in determining the location of the 
MBS 

30 deadreck_est The deadreckoning MBS location change estimate paired with 

u MBS_new_est n . 
Returns a potentially updated u MBS_curr_est" */ 

{ 

if (MBS_new_est is of type MANUAL_ENTRY) then 
35 { /* MBS operator entered an MBS location estimate for the MBS; so must use it */ 
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insert "MBS_new_est n as the head of the "current MBS location track" which is the location track 
indicating the best current approximation of the location of the MBS; 

} 

else /* M MBS_new_est n is not a manual entry */ 
{ 

MBS_curr_est <— geLcwr_esf(MBS_new_est.MS_ID); /* get the head of the "current location 
track"*/ 

adjusted_curr_est <-- applyjcleadm(XoningJo(MB$_curr_es\, deadreck_est); 

./* The above function returns an object of the same type as "MBS_curr_esr, but with 
the most likely MBS point and area locations adjusted by H deadreck_esf. 
Accordingly, this function performs the following computations: 

(a) selects, Ambs, the MBS location area estimate of H MBS_curr_esf (e.g., one 
of the "most likely" nested area(s) provided by "MBS_curr_est n in one 
embodiment of the present invention); 

(b) applies the deadreckoning translation corresponding to "deadreck_est n to 
Ambs to thereby translate it (and expand it to at least account for 
deadreckoning inaccuracies). */ 

if (reasonaWy_c/ose(MBS_new_est, adjusted_curr_est, MBS_curr_est)) 

r In one embodiment, the function B reasonably_close" here determines whether a most 
likely MBS point location (i.e., centroid) of "MBS_new_esf is contained in the MBS 
estimated area of M adjusted_curr_esf 

Note that the reasoning for this constraint is that if fl MBS_curr_est" was accurate, 
then any "most likely MBS point location" of a new MBS baseline estimate that is also 
accurate ought to be in the MBS estimated area of "adjusted_curr_esr. 

In a second embodiment, the function "reasonably_close" determines whether 
the centroid (or most likely MBS point location) of B MBS_new_est n is close enough to 
"MBS.curL.esf so that no MBS movement constraints are (grossly) violated between 
the most likely point locations of "MBS_new_esf and "MBS_curr_esr; i.e., 
constraints on (de)acceleration, abruptness of direction change, velocity change, max 
velocity for the terrain. Note, such constraints are discussed in more detail in the 
section herein describing the "Analytical Reasoned. Accordingly, it is an aspect of 
the present invention to provide similar capabilities to that of the Analytical Reasoner 
as part of the MBS, and in particular, as the functionality of the "MBS LOCATION 
CONSTRAINT CHECKER" illustrated in Fig. 1 1 . It is assumed hereinafter that the 
embodiment of the function, u reasonably_close\ performed here is a combination of 
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both the first and second embodiments, wherein the constraints of both the first and 
second embodiments must be satisfied for the function to return TRUE. */ 

then 

if (the confidence in MBS„new_est >= the confidence in MBS_curr_est) then 

{ 

if (the most likely MBS area of MBS_new_est contains the most likely MBS area of 
"adjusted_curr_esr as computed above) then 

shrink MBS_new_est uniformly about its centroid (i.e., "most likely MBS point 
location") until it is as small as possible and still contain the MBS 
estimated area of u adjusted_curr_esf. 
insertJntoJocationJrack^cuTTenV, MBS_new_est); 

I* The program invoked here inserts a location track entry corresponding to the second 
parameter into the location track identified by the first parameter (e.g., "current"). It 
is important to note that the second parameter for this program may be either of the 
following data structures: a "location track entry", or an "MBS location estimate" and 
the appropriate location track entry or entries will be put on the location track 
corresponding to the first parameter. The insertion is performed so that a latest 
timestamp" order is maintained; i.e., 

(a) any extrapolation entries in the location track, wherein these entries 
have a more recent latest timestamp" than the ("earliest" or only) 
timestamp (depending on the data structure) of the second parameter 
are removed, and 

(b) conceptually at least, the location change estimates output from the 

deadreckoning MBS location estimator that correspond with the 
removed extrapolation entries are then reapplied in timestamp order 
to the head of the target location track. 7 

} 

else r the centroid of "MBS_new_est", is contained in an area of "MBS_curr_esr, but the 
confidence in "MBS_new_est" < confidence in "MBS_curr_est" 7 

{ 

most_likely_est <- determine a "most likely MBS location estimate" using the set S 
= {the MBS location estimate centroid(s) of any MBS location 
track heads contained in the MBS estimated area of 
t, adjusted_curr_est B , plus, the centroid of "MBS_new_esf }; 
r Note, in the above statement, the "most likely MBS location estimate" 
may be determined using a number of different techniques depending on 
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what function(s) is used to embody the meaning of "most likely". In one 
embodiment, such a "most likely" function is a function of the confidence 
values of a predetermined population of measurements (e.g., the 
selected location track heads in this case) from which a "most likely" 
measurement is determined (e.g., computed or selected). For example, 
in one embodiment, a "most likely" function may include selecting a 
measurement having the maximum confidence value from among the 
population of measurements. In a second embodiment, a "most likely" 
function may include a weighting of measurements (e.g., location track 
heads) according to corresponding confidence values of the 
measurements. For example, in the present context (of MBS location 
track heads) the following steps provide an embodiment of a "most likely" 
function: 

(a) determine a centroid of area for each of the selected track heads 
(i.e., the location track heads having a point location estimate 
contained in the MBS estimated area of H adjusted_curr_est"); 

(b) determine the "most likely location MBS position" P as a weighted 
centroid of the centroids from step (a), wherein the weighting of 
each of the centroids from (a) is provided by their corresponding 
confidence values; 

(c) output an area, Ai, as the "most likely MBS location area", wherein 
the centroid of Ai is P and Ai is the largest area within the MBS 
estimated area of "adjusted_curr_est" satisfying this condition; and 

(d) set a confidence value for Ai as the average confidence value of 
"MBS_new_esf , u MBS_cuiT_est" and the selected location track 
head used. */ 

insert JntoJocationJrack( a current n , most_likely_est); 

} 

else r B MBS_new_esf is not reasonably close to "adjusted_curr_esf (i.e., u MBS„curr_est" with 
"deadreck_est" applied to it), so a conflict exists here; e.g., (i) "MBS_new_est" is not a 
manual entry, and (ii) "MBSjiew_est n does not have its centroid contained in the MBS 
estimated area of "adjusted_curr_esf , or, there has been a movement constraint 
violation. Note that it is not advisable to just replace H MBS_curr_est n with 
"new_est_head" because: 

(a) w MBS_new_est" may be the MBS location estimate that is least accurate, 
while the previous entries of the current location track have been accurate; 
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(b) the "MBS_curr_esf may be based on a recent MBS operator manual entry 
which should not be overridden. */ 

{ 

MBS_curr_est <— resolve_conflicts(MBS_new_est, adjusted_curr_est, MBS_curr_est); 

5 } 

} /* end else "MBS_new_esf not a manual entry */ 
if (MBS is a vehicle) and (not off road) then 

/* it is assumed that a vehicular MBS is on-road unless explicitly indicated otherwise by 
MBS operator. */ 

1 o MBS_curr_est <— snap JojbesfJS_$freef(MBS_curr_est); /* snap to best street location 

according to location estimate, velocity, and/or direction of travel. Note, this is a translation 
of a MBS_curr_esf. */ 
RETURN(MBS_curr_est) 
} r END UPDATE(MBS_CURR_EST) */ 

15 

resolve_conflicts(MBSjiew_est, adjusted jcurr_est, MBSjcurrjest) 

r There is a basic conflict here, 

(i) "MBS_new_esr is not a manual entry, and 

(ii) one of the following is true: u MBS_new_est B does not have its centroid contained in the 
20 area "adjusted_curr_esP, or, using a MBS_new_est" implies an MBS movement constraint 

violation. 

Input: MBS_new_est The newest MBS location estimate record. 

adjusted_curr_estThe version of a MBS_curr_esf adjusted by the deadreckoning 
location change estimate paired with "MBS_new_esf. 
25 MBS_curr_est The location track entry that is the head of the "current" 

location track. Note that a MBS_new_est.confidence n > 
"MBS_curr_est.confidence B . 
Output: An updated tt MBS_curr_esf\ */ 

{ 

30 mark that a conflict has arisen between tt MBS_curr_esf and w MBS_new_esf; 

if (the MBS operator desires notification of MBS location estimate conflicts) then 

notify the MBS operator of an MBS location estimate conflict; 
if (the MBS operator has configured the MBS location system to ignore new estimates that are not 
"reasonably close" to adjusted_curr_est) or 
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(MBS_curr_est is based on a manual MBS operator location estimate, and the MBS has moved 
less than a predetermined distance (wheel turns) from where the manual estimate was provided) 
then 

RETURN(adjusted_curr_est); 
5 else I* not required to ignore "MBS_new_esf t and there has been no recent manual estimate input*/ 

{ /* try to use u MBS_new_est n */ 

if ((MBS_new_est.confidence - adjusted_curr_est.confidence) > a large predetermined threshold) 
then 

I* Note, the confidence discrepancy is great enough so that "MBS_new_est B should be 
1 o the most recent baseline estimate on current MBS location track. Note that the 

threshold here may be approximately 0.3, wherein confidences are in the range [0, 
1]*/ 

insertJntoJocationJrack( u currer)\ n , MBS_new_est); 

r insert tt MBS_new_est B into "current" location track (as a baseline entry) in latest 
15 timestamp" order; i.e., remove any extrapolation entries with a more recent 

latest timestamp" in this track, and reapply, in timestamp order, the location 
change estimates output from the deadreckoning MBS location estimator that 
correspond with the removed extrapolation entries removed; */ 
else r 0 MBS_new__est.confidence ,, is not substantially bigger than 
20 "adjusted^unLest-confidence"; so check to see if there are potentially MBS location 

system instabilities */ 
{ r check for instabilities 7 

if [ (there has been more than a determined fraction of conflicts between the "MBS_curr_est n 
and tt MBS_new_esf within a predetermined number of most recent "MBSjiew_esf 
25 instantiations) or 

(the path corresponding to the entries of the "current location track" of the MBS has 
recently violated MBS movement constraints more than a predetermined fraction of 
the number of times there has been new instantiation of u MBS_curr_esf, wherein 
such movement constraints may be (de)acceleration constraints, abrupt change in 
30 direction constraints, constraints relating to too high a velocity for a terrain) or 

(there has been an MBS operator indication of lack of confidence in the recently 

displayed MBS location estimates)] 
then f the MBS location system is likely unstable and/or inaccurate; check to see if this 
condition has been addressed in the recent past. */ 
35 {/* fix instability */ 
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if (fix_instability_counter equal to 0) then f no instabilities have been addressed here 
within the recent past; i.e., "fix_instability_counter" has the following semantics: 
if it is 0, then no instabilities have been addressed here within the recent past; 
else if not 0, then a recent instability has been attempted to be fixed here. 
Note, "fixjnstability.counter" is decremented, if not zero, each time a new 
baseline location entry is inserted into its corresponding baseline location 
track. Thus, this counter provides a "wait and see" strategy to determine if a 
previous performance of the statements below mitigated the (any) MBS 
location system instability. */ 

{ 

most_Jikely_est <- determine a new "most likely MBS location estimate"; [30.1] 
T Note, a number of MBS location estimates may be generated and 
compared here for determining the u mostJikely_esf. For example, 
various weighted centroid MBS location estimates may be determined by 
a clustering of location track head entries in various ways. 

In a first embodiment for determining a value (object) for 
"mostJikely_est", a "most likely" function may be performed, wherein a 
weighting of location track heads according to their corresponding 
confidence values is performed. For example, the following steps provide 
an embodiment of a "most likely" function: 

(a) obtain a set S having: (i) a centroid of area for each of the 
track heads having a corresponding area contained in a 
determined area surrounding the point location of 
"adjusted_curr_est n (e.g., the MBS estimated area of 
"adjusted„curr_esf ), plus (ii) the centroid of "MBS_new_esf; 

(b) determine the "most likely location MBS position" P as a 
weighted centroid of the centroids of the set S from step (a), 
wherein the weighting of each of the centroids from (a) is 
provided by their corresponding confidence values; 

(c) output an area, A, as the "most likely MBS location area" 
wherein A has P as a centroid and A is a "small" area (e.g., a 
convex hull) containing the corresponding the centroids of the 
set S; and 

(d) set a confidence value for A as the average confidence value 
of the centroids of the set S. 
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In a second embodiment, "mostJikely_esr may be determined by 
expanding (e.g., substantially uniformly in all directions) the MBS location 
estimate area of tt MBS_new_est n until the resulting expanded area 
contains at least the most likely point location of u adjusted_curr_est" as its 
most likely MBS location area. */ 
insert JntoJocation_track( tt cumn\", most_likely_est); 
fixJnstability_counter <— a predetermined number, C, corresponding to a 
number of baseline entries to be put on the baseline location tracks until 
MBS location system instabilities are to be addressed again here; /* when 
this counter goes to zero and the MBS location system is unstable, then the 
above statements above will be performed again. Note, this counter must 
be reset to C (or higher) if a manual MBS estimate is entered. */ 

} 

}/* fix instability*/ 

else T The MBS location system has been reasonably stable, and 
B MBS_curr_est.confidence" is not substantially bigger than 
u adjusted_new_est.confidence B . 7 

mostJikely_est <- determine a most likely MBS location estimate; 
/* The determination in the statement above may be similar or substantially the same as 
the computation discussed in relation to statement [30.1] above. However, since 
there is both more stability in this case than in [30.1] and less confidence in 
"MBSjiew^est", certain MBS movement constraints may be more applicable here 
than in [30.1]. 

Accordingly, note that in any embodiment for determining a most_likely_est B 
here, reasonable movement constraints may also be used such as: (a) unless 
indicated otherwise, an MBS vehicle will be assumed to be on a road, (b) a new 
MBS location estimate should not imply that the MBS had to travel faster than, for 
example, 120 mph or change direction too abruptly or change velocity too abruptly 
or traverse a roadless region (e.g., com field or river) at an inappropriate rate of 
speed. 

Thus, once a tentative MBS location estimate (e.g., such as in the steps of the 
first embodiment of [30.1]) for "mostjikely.est" has been determined, such 
constraints may be applied to the tentative estimate for determining whether it 
should be pulled back toward the centroid of the "MBS.curr.est" in order to satisfy 
the movement constraints*/ 

105 



insert JntoJocationJrack( u current n , most_likely_est); /* note, the second parameter for 
this function may be either of the following data structures: a "location track entry", or 
a "MBS location estimate" and the appropriate location track entry or entries will be 
put on the location track corresponding to the first parameter. */ 

} 

} r check for instabilities */ 

MBS_curr_est <— geLcurr_esf(MBSjiew_est.MSJD); /* from current location track */ 
} /* try to use tt MBS_new_est n */ 
RETURN(MBS_curr_est) 
} /* END resolve_conflicts */ 
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